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The L/K-Capture Ratio of Germanium 71+ 


By R. W. P. Drever and A. MousK 
Natural Philosophy Department, University of Glasgow 


[Received November 9, 1956] 


ABSTRACT 


A direct measurement of the L/K-capture ratio of 74Ge has been made 
using a special proportional counter system which avoids the effects of 
X-ray escape. One proportional counter is effectively enclosed by another, 
without any intermediate wall. With a gaseous source of germanium 
hydride and an argon filling at 6 atmospheres, the ratio of the L- and 
K-peaks from the central counter, in anticoincidence with the surrounding 
counter, was determined. The L/K-capture ratio was found te be 
0-1287 'h03- In addition, in’ a separate measurement, a value of 
0-53-+0-03 was obtained for the K-fluorescence yield of gallium. 


§ 1. INTRODUCTION 


THE ratio of the probabilities of decay by the capture of orbital electrons 
from the L-shell and K-shell has been determined experimentally for a 
number of nuclei, but the data available are still insufficient to provide a 
satisfactory check of the relevant theory. The situation has been 
reviewed recently by Radvanyi (1955) and by Robinson and Fink (1955), 
and very few nuclei were found for which experimental results could be 
compared with theory. 

In general there are large discrepancies between calculated and 
measured values of the L/K-capture ratio. Relatively good agreement is 
found only for 37A and 7°Kr, in the latter case after a recalculation of the 
experimental results with a revised value of fluorescence yield. In some 
of the other nuclei the disagreement may be due to uncertainties in decay 
schemes. A particularly large discrepancy in the case of germanium 71 
cannot be explained in this way, however, since the decay is simple and 
well established. A further study of electron capture in this nucleus is 
therefore desirable and this is the subject of the present work. 

Germanium 71 is a pure electron capture source with a half-life of 12-5 
days (Bisi et al. 1955). The internal bremsstrahlung spectrum associated 
with the electron capture was studied first by Saraf et al. (1953) and Saraf 
(1954) who found agreement with the theoretical shape from the end point 
down to 100 kev. Langevin (1954 b) observed a low energy component 
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with longer half-life in his measured spectrum and suggested that the 
decay might be complex. However more recent work by Bisi et al. 
(1955) has shown that this effect was due to the presence of a radioactive 
contaminant, since their spectrum agreed closely with theory down to 
50 kev when well purified sources were used. From the measured end | 
point of the bremsstrahlung spectrum the transition energy is found to 
be 231 kev (Bisi et al. 1955) and the decay is allowed with log f7=4-6. 
The theoretical value for the L/K-capture ratio of 71Ge obtained from 
the calculations of Brysk and Rose (1955) with these data is 0-106, 
Experimentally a higher value of 0-30 has been found by Langevin (1954 b 
and 1955) using a proportional counter containing the source in gaseous 
form. The method employed in the present experiment, although some- 
what different in principle, also depends on the measurement of a gaseous 
source with a proportional counter. It may therefore be useful to give a 
brief review of some of the features common to this type of measurement. 
In a determination of an L/K-capture ratio with a gaseous source 
proportional counter the intensities NV, of the L-peak and NV, of the 
K-peak are measured. The observed K-peak is due to Auger electrons 
and to those x-rays emitted following K-capture that are absorbed in the 
filling gas. The L-peak is due to Auger electrons and x-rays following 
L-capture events as well as to L-radiation emitted, after K-capture, 
simultaneously with Ka x-rays that escape undetected from the counter. 
If F is the K-fluorescence yield of the daughter element, & the fraction of 
Ke x-rays in the K series and P the probability of a K x-ray escaping 
from the counter without being detected, then the L/K-capture ratio R 
is given by 
Ny, 


R= N, 


(1—PF)—PFk. 
The escape of L x-rays from the counter is usually small, and is neglected 
here. 

In practice it is convenient to choose the dimensions of the counter and 
the filling gas so that the probability P of x-ray escape approaches either 
1 or 0. The former method is experimentally the simpler one. A small 
low pressure counter is used, preferably with a filling gas of low atomic 
number. ‘The measurement of the L/K-capture ratio of 71Ge by Langevin 
was performed in this way, propane being used as filling. Unfortunately, 
the wide scatter of measured values of fluorescence yield introduces some 
uncertainty into the result. Langevin took for the fluorescence yield 
of gallium the value 0-45 given by a semiempirical formula by Burhop 
(1952). However more recent values of fluorescence yield suggested by 
Broyles ct al. (1953) and by Burhop (1955) are higher, and results of 
calculations by Rubenstein and Snyder (1955) and measurements by 
Heintze (1955) for argon and krypton show the same tendency. With 
the value of 0-50 of Broyles et al. the L/K-capture ratio found by Langevin 
is brought down from 0-30 to 0-19. 
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The alternative method of measuring an L/K-capture ratio, in which 
P=0, gives a result independent of the fluorescence yield and is therefore 
to be preferred. All radiations resulting from a K-capture or an L-capture 
event are integrated and recorded only in the corresponding K- or L-peak, 
so that the L/K-capture ratio is given directly. This method has been 
employed in the case of *’A by Pontecorvo et al. (1949), and recently by 
Langevin and Radvanyi (1955). Its application to 7!Ge is more difficult. 
since the half-thickness for absorption of the 9-2 kev K x-rays in argon 
at 1 atmosphere is 5-0 cm. For an escape correction to the final result 
of 5°, the diameter of a counter filled with argon to 6 atmospheres would 
have to be 1:0 metre. However with a technique developed in the present 
work the same accuracy can be achieved in a practicable counting system 
of moderate dimensions. 


§2. OUTLINE OF THE METHOD 


The arrangement consists effectively of two proportional counters, 
one inside the other, connected in anticoincidence. A diagram of the 
system is given in fig. 1. The inner circle of wires acts as cathode for the 
central proportional counter, while alternate anode wires in the outer 


Fig. 1 
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Internal construction of the counter. 


circle turn the surrounding layer of gas into a separate counter. With 
a gaseous source, and an argon filling at 6 atmospheres, the spectrum of 
pulses from the central counter in anticoincidence with pulses from the 
ring counter is measured. It may be shown that the ratio of the intensities 
of the L- and K-peaks then observed is closely equal to the L/K-capture 
ratio. 

It should be noted first that with the filling used the half-thickness 
for absorption of K x-rays is 0-81 cm, and the probability of an X-ray 
passing out from the central counter and not being absorbed and detected 
in the ring counter is less than 0-1%. It will be neglected in this pre- 
liminary discussion, as will end effects. 
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All pulses due to L-capture events in the central counter are recorded 
in the L-peak. Incomplete absorption of radiations following a K-capture 
event in the central counter does not lead to the recording of a spurious 
L-count, since the escaping Ka x-ray is detected in the surrounding ring 
counter and closes the gating circuit. Thus the number N,, of counts 
in the L-peak is equal to the number of L-capture events that have taken 
‘place in the central counter. 

The total number of counts in the recorded K-peak is also equal to the 
number of K-capture events that have occurred in the central counter, 
although in this case the peak is made up of two separate components. 
One of these is due to K-capture events that take place, and are also 
detected in the central counter. The other component is due to K x-rays 
which are produced by K-capture in the surrounding layer of gas, pass 
into the central counter, and are absorbed there. Since the layer of gas 
is assumed to be infinitely thick, these x-rays compensate for the X-rays 
that escape from the central counter, provided the gating circuit is biased 
so that pulses corresponding to the L-absorption energy do not close it. 
The two components of the peak differ slightly in energy since the latter 
is caused by K x-rays while the former has the total K-absorption energy. 
In practice they are not resolved and the area of the complex peak gives 
the total number of K-capture events in the central counter. 

The L/K-capture ratio is therefore given directly by the ratio N;,/\V x, 
though some small corrections are required for end effect and inefficiency 
of the ring counter. 


§ 3. CONSTRUCTION AND OPERATION OF THE COUNTING SYSTEM 


The wires making up the ring counter system are mounted between 
‘brass plates which are held apart by five supporting rods. The whole 
cage-like system can be removed from the containing vessel as a unit, 
which makes the fitting of the 20 wires relatively simple. The 10 cathode 
wires of the inner circle, and the 5 in the outer circle, are stainless steel 
0-2 mm in diameter. The 5 counting wires, 0-07 mm in diameter, are 
fitted with guard tubes and are connected together inside the counter, the 
connections also being screened by an annular metal box. The central 
counter is end-corrected and has a sensitive length of 76-0cm. Its 
diameter is 4:5cm while the ring counters are slightly larger. Both 
central and ring counters operate at voltages in the region of 6-8 kv 
when filled to 6 atmospheres. 

Pulses from the central counter, after amplication pass through a 
60 psec delay line to a linear gate circuit and thence to a multichannel 
kicksorter, as shown in fig. 2. From the ring counter, pulses are fed 
to an amplifier, discriminator, and lengthener controlling the gate. 

The counting rate in the ring counter is 15 times greater than that in the 
central counter. However with the electronic system shown, the high 
counting rate in the ring counter does not affect the efficiency of gating 
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Since the circuits controlling the gate are sensitive the whole time, even 

immediately after a pulse. ‘There may be a reduction in the effective 
counting time but this does not alter the ratio of the measured intensities 
of the L- and K-peaks since the relative loss of time is the same in both 
cases. It is necessary that the time constant of the lengthener does not 
change and this condition is checked by a timing unit shown in the 
diagram, which measures directly the effective counting time, by counting 
pulses of constant repetition rate with a gated scaler. 


Fig. 2 
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Schematic diagram of the electronic equipment. 


Fig. 3 
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In view of the unusual design of the counter system its operation was 
checked by measuring the pulse distribution of the K-peaks given by the 


central centre and ring counters with a 7!Ge gaseous source. In the case 
of the central counter the symmetrical arrangement of the many wires 
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forming the cathodes leads to a field distribution which differs little from 
that obtained with a complete cylindrical electrode and the counter 
was found to behave normally. In the ring counters, the cathode wires 
are arranged so that the field distribution around the counting wires is 
distorted as little as possible. Satisfactory performance is obtained, as 
is shown by the 10-2 kev peak in fig. 3, which indicates in addition that 
there is no appreciable variation of gain between individual counting 
units. 

For an accurate measurement of the L/K-capture ratio it is essential 
that the ring counter system has a high efficiency, in that every K x-ray 
absorbed in the gas causes a pulse large enough to trigger the gate circuit. 
The small amplitude of the low energy tail of the peak in fig. 3 gives 
evidence that the counter fulfills this condition sufficiently well. The 
major part of the tail is caused by end effect due to the absence of field 
correcting tubes, and by the natural background. 

It may be remarked that the movement of the ions in the central 
counter induces small pulses of opposite sign on the anode wires of the 
ring counter system. Their amplitude is about 1/40 of that of the corre- 
sponding normal pulses, so they have no effect on the measurements. 


§ 4. EXPERIMENTAL PROCEDURE AND RESULTS 


Germanium 71 was made by (n,y) reaction on germanium oxide at 
A.E.R.E. Harwell, and was stored for a time long enough to allow the 
12-hour activity of ““Ge to decay. Germanium hydride gas was then 
prepared by the action of sodium amalgam on a solution of the germanium 
oxide. The gas was dried by cooling with dry ice, and hydrogen was 
pumped off with the germanium hydride frozen at liquid air temperature. 

The germanium hydride was found to have no effect on the operation of 
the counter at the low partial pressure of less than 0-1 mm of mercury 
used. The absence of residual activity after pumping out the counter 
indicated that there was no chemical reaction with the counter materials. 

Measurements were carried out with argon fillings at various pressures 
up to a maximum of 6 atmospheres, this limit being set by the mechanical 
strength of the counter. Methane at a pressure of 15 cm of mercury was 
used as quenching gas. At each pressure, pulse spectra in the region 
of the 10-2 kev K-peak and the 1:3 kev L-peak were analysed. 
The change in energy range was made by altering the amplifier gain, all 
other conditions being kept constant. The gate control circuit was 
biased at 3 kev throughout. 

Typical spectra of K- and L-peaks obtained with the central counter in 
anticoincidence with the ring counter are given in figs. 4 and 5. In these 
diagrams the background, which was measured separately, has not been 
subtracted. 

The half-width of the K-peak shown is slightly greater than would 
normally be expected at this energy due to the presence of the two 
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components, of K-absorption energy and K x-ray energy, and also partly to 
a technical fault in the assembly of the counter which made it difficult to 
centre the anode wire accurately. However the area of the peak is still 
well defined and this broadening does not affect the measurement. 

The ratio of the counting rates observed at various pressures from 2 to 
6 atmospheres are given in the table. The small corrections required for 
kicksorter dead-time, which was 0-6 msec, have been applied. 


Fig. 4 
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At the lower pressures the condition of negligible x-ray escape through 
the ring counter, assumed above, is not fulfilled and high values of the 
L/K-ratio are obtained. At higher pressures when escape through the 
ring becomes small there is little variation with pressure. However even 
at 6 atmospheres some corrections have still to be applied to find the true 
value of the L/K-capture ratio since the shielding of the central counter 
by the ring counter is not quite complete. 

K x-rays which leave the sensitive volume of the central counter 
through the ends are not detected. An approximate calculation by the 
methods used for shielding problems (see for example Glasstone 1956) 
shows that at 6 atmospheres the fraction P, of K x-rays emitted in the 
central counter that escape in this way is P}=710-3. Further, some 
X-rays are prevented from entering the ring counter by the wires making 
up the cathode. The calculated probability P, for this is P,=7 x 107%. 
The remaining correction, for x-rays passing through the layer of gas 
forming the ring counter without being absorbed, is relatively small at 
6 atmospheres. The width of the ring is then 7-6 half-thicknesses, and 
the probability P, for escape is P,=6 x 10-4. 


Pressure in atmospheres 2 3 4 5 6 


Observed L/K-ratio (°%) 16-5 14-4 13-6 13-5 13-5 


When these corrections are taken into account, the L/K-capture 
ratio R is obtained from the measured ratio V,/N, by the expression 


N,, 

Nx 

A mean of the values of V/V, obtained from three separate measure- 
ments gives a value for the L/K-capture ratio of germanium 71 of 


R 


[1—(P,+P,)F]—(Pi4-P.+P,) Fk. 


The statistical error of the measurement is 1°. However there is a 
possibility of a systematic error in establishing the shape of the low 
energy part of the L-peak, and to a lesser extent of the K-peak. With 
an estimate of this error, the overall accuracy of the final result may be 
stated as (15)%. 

In calculating the L/K-capture ratio above the values taken for the 
K-fluorescence yield of gallium and for the fraction of Kx x-rays in the 
K-series were 0-50 and 0-86 respectively. However the result obtained 
is very insensitive to these factors, a change of 10°/ in either of them 


/0 
leading to a change in the L/K-capture ratio of only 0-5%. 


§ 5. Discussion 


The L/K-capture ratio obtained here is considerably lower than that 
found by Langevin. However, as mentioned above, it is probable that 
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his high result is due in part to the value taken for the fluorescence vield 
of gallium. 

In an attempt to check this explanation, the counter described here 
has been used also to make an independent measurement of the 
flourescence yield. The experiment is described briefly in the Appendix. 
The value of the fluorescence yield obtained was 0:53 with an estimated 
accuracy of 6°, and in addition the fraction of Ka x-rays in the K-series 
was found to be 0-83. With these values the experimental result of 
Langevin comes into closer agreement with the result of the present 
measurement. 

The theoretical calculation of the L/K-capture ratio originally carried 
out by Marshak (1942) involves ratios of the components of L- and 
K-electron wave functions. These are becoming more accurately known 
now and have been recently computed by Rose and Jackson (1949) and 
by Brysk and Rose (1955). The L/K-capture ratio of 7!Ge obtained 
from the results of the latter authors is 0-106, which is appreciably 
lower than the value found in the present measurement. 

It seems difficult to explain the difference in the results by a systematic 
error in the experiment, and it is possible that the discrepancy arises in 
the calculation of the theoretical value. Recent work by Odiot and 
Daudel (1956) shows that when correlations between the positions of the 
electrons are taken into account higher values of the L/K-capture ratio 
result, but at present correction factors are available only for helium, 
beryllium, and argon. 


APPENDIX 
The K-fluorescence Yield of Gallium 


The fluorescence yield was measured, independently of the L/K-capture 
ratio, using the counter described above with the anticoincidence circuit 
disconnected. The counter was filled with methane to a pressure of 
1 atmosphere together with a trace of germanium hydride source. The 
intensity of the K-peak in the central counter, which in this case was due 
to Auger electrons only, was measured. Argon was then added to bring 
the total pressure to 5 atmospheres and the intensity of the K-peak, due 
to both Auger electrons and x-rays was measured again. 

After applying corrections for the absorption of K x-rays in the methane, 
which has a probability of 2-5°,, and for the escape of K x-rays when the 
argon filling is used, which has a probability of 0-89, the fluorescence 
yield of gallium was found to be 0-53. The accuracy is estimated as 6". 

In addition the fraction of K« x-rays in the K-series was determined 
from the intensity of the L-peak observed with the methane filling. 
After subtracting the component due to L-capture, calculated from the 
known L/K-capture ratio, the rate of emission of Ka x-rays is obtained. 
With similar corrections to those mentioned above, a value of 0°83 was 
found for the fraction of Ke x-rays in the K-series of gallium. However 
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the difference between this value and the value of 0-865 used by Langevin 
is not considered significant. 
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Some Observations on the Nature of Fatigue Damagey 
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§ 1. INTRODUCTION 
Ir has been observed in a number of studies on the fatigue behaviour 
of metals that fatigue cracks are initiated along slip bands. This paper 
describes some observations on the nature of these bands at the stage 
when fatigue cracks are forming. 

When the surface steps or grooves of slip are removed by electrolytic 
polishing, faint residual traces may remain or may be revealed by 
subsequent etching. These traces may be revealed in certain materials even 
after unidirectional strain (Jacquet 1945, McLean 1948, Samuels 1954). 
If slip occurred by a shear displacement only, then it should not be 
possible to reveal these bands. However, dislocations along a slip 
plane may make this a region of preferential etching if a reagent sensitive 
to these differences of lattice disorder is used. Certain reagents may 
produce rows of pits along slip bands and these pits have been interpreted 
as sites of discrete dislocations. Apart from the presence of disturbed 
lattice, other local changes may allow the slip traces to be revealed 
(Forsyth and Stubbington 1955). These changes may be very local 
recrystallization or polygonization, or slip band ageing accelerated by 
strain. Another suggested cause for the persistence of slip traces is the 
diffusion of oxygen down slip planes during a fatigue test (Wadsworth 
and Thompson 1954). 

All of these effects are more likely to occur with fatigue stresses and to 
become more prominent with time. They are therefore most likely to 
be visible along slip bands which are in the stage of development into 
cracks, particularly near the ends of the initial cracks. 

These bands are the only signs of permanent fatigue damage which 
can be detected microscopically. A particular case of residual bands 
or zones in aluminium 4°% copper alloy fatigued at room temperature, 
has already been reported (Forsyth and Stubbington 1955). In order 
to reveal these bands specimens were etched in 10°, NaOH solution, 
a relatively drastic etchant which no doubt greatly widened the affected 
fatigue zones. 


§ 2. DamacsE Banps In ALUMINIUM 7:5% Zinc 25% MAGNESIUM ALLOY 

A number of specimens of the above alloy have been fatigued in 
torsion using a pendulum torsional vibrator. It has been observed 
that slip traces can be revealed by subsequent electrolytic polishing 
alone. The material was tested in the following condition : 
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Solution heat treated 450°c for 1 hour, cold water quenched, then aged 
at 150°c for 4 hour unless otherwise indicated. This treatment is not 
long enough to age harden the alloy to peak hardness but brings the 
hardness from 80 d.p.h. (solution treated condition) to 150 d.p.h., whereas 
peak hardness is about 185 d.p.h. 

It was found that the effects which will be described also occurred 
in the solution treated and over-aged conditions but were less marked. 


Fig. 3 


A schematic illustration of a damage band with tubular holes distributed 
along it. 


The zones which were observed are illustrated in fig. 1 (Pl.). They 
appear to be due to more rapid polishing along the slip traces and are 
best observed by oblique lighting. A series of holes are also revealed 
along these zones, as indicated by arrow A, fig. 1 (Pl.) and arrow A, 
fig. 2(PL.). It appears that fatigue cracks develop by breaking through 
these holes and it is very noticeable that they are very evenly spaced. 
The mean spacing measured on a number of specimens being 0-36 pw. 
The distribution of these holes in the slip plane cannot as yet be explained. 
In spite of the random sectioning of the crystal planes which necessarily 
occurs with a polycrystalline specimen, there is a consistent spacing 
of holes from grain to grain. This suggests that the holes are in some 
way aligned with the surface and may exist as tubes normal to the 
specimen surface, the marked decrease in the number of holes visible 
as the surface is removed being due to the different lengths of tube. 
By sectioning specimens in a direction normal to the surface it is sometimes 
possible to observe these tubes as shown in fig. 2 (Pl.) (arrow B) although 
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due to their small size it is only rarely that they can be revealed in length 
by a random section. Figure 3 illustrates in schematic form a possible 
arrangement of these holes. The direction of oblique lighting is denoted 
by X normal to the slip trace and Y parallel to the slip trace. It can 
be seen from this diagram that X direction lighting is necessary to reveal 
clearly the grooves whereas Y direction lighting more clearly reveals the 
holes as there will be no obscuring shadow from the groove in this direction. 

The fact that the slip zones are preferentially dissolved by the electrolyte 
suggests that these regions are more anodic than the general matrix, 
perhaps due to less stable anodic films forming at these places during 
the electropolishing process. It had been found in earlier work that 
zones depleted of solute atoms such as occur at the grain boundaries and 
subgrain boundaries in this alloy are less readily attacked and therefore 
appear as ridges after electropolishing. This is shown in fig. 4 (PI1.). 
Arrow A indicates an artificially made reference groove which is included 
to show that the boundary zones are in fact ridges. This specimen had 
been electropolished, at which stage only the ridges could be observed. 
It was subsequently etched for 20 seconds in Wassermann’s reagent 
(25 cm* 10% H,Cr,0, 10 cm? HCl, 100 cm? NHO,, 2-7 cm? HF (60%) 
500 em? H,O) which then revealed the precipitate. The electropolishing 
alone will only reveal large precipitates which are formed after much 
longer ageing times. It is therefore suggested that the pits or holes which 
appear in the © as polished ’ condition in the fatigued specimens are not 
precipitates which have been dissolved out, but are in fact voids along the 
slip planes. This has now been confirmed by the fact that these holes 
are clearly visible after mechanical polishing of sections of these specimens. 
It might also be suggested that as the slip traces electropolish as grooves 
whereas the depleted zones at the grain boundaries appear as ridges, 
the slip zones are richer in solute atoms than the surrounding matrix, 
and therefore more anodic. However, this is not necessarily the explana- 
tion as the presence of lattice defects would also be expected to make these 
zones more anodic than the less disturbed matrix. 

If the holes along the slip zones are in fact due to the condensation of 
vacancies then it may be that the zones represent regions of fairly dense 
vacancy distribution where the defects are still in free motion and not 
yet attached to the holes. However, this does not seem likely as the slip 
zones in solution heat treated material do not etch as grooves ; they are 
in fact extremely difficult to reveal but there is some evidence that they 
form ridges when electropolished. Experiments with other alloys and 
with pure aluminium have shown that the formation of holes seems to be 
a common occurrence although the grooves along slip striations are 
only clearly observed in aged alloys. 

These slip zones and holes occur in profusion in the surface layers 
of the specimen but extend only a few microns into the surface grains, 
except for a few actively propagating cracks. Sections taken through 
fatigue specimens revealed that the cracks advance by the formation 
of holes along advancing slip zones. This is shown in fig. 5 (PI.). Smal 
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pits of a similar spacing to the holes found in the surface layers can be 
detected on the as polished surface of the fatigue specimens with no 
subsequent polishing. This is shown in fig. 6 (Pl.). This specimen was 
fatigued in the solution treated condition. That these pits represent 
the holes as they appear along the slip band itself is not absolutely certain, 
but it seems probable that this isso. Slip band extrusion occurs profusely 
along these bands, and if the surface is then repolished and refatigued the 
extrusion often occurs as narrow filaments from the slip zones between 
holes. 

A subsequent solution treatment of specimens which had been fatigued 
(450°c for + hr) removed the slip zones but did not heal all of the holes. 
It is not possible to ascertain whether any of the holes were healed as 
repolishing is necessary after the heat treatment, and there is such an 
intense density gradient of holes from the surface that no numerical 
comparison can be made. The fact that re-solution does not necessarily 
remove the holes also suggests that they are not the sites of precipitate. 


§ 3. CONCLUSIONS 

Fatigue stresses produce damage in such materials as aluminium 
7:5°, zine 2-5°,, magnesium in the form of persistent bands which at 
an advanced stage contain rows of holes. 

The bands can be removed by heat treatment, but the holes are not 
healed. 

There is some evidence that these holes are in the form of tubes, 
penetrating to different depths in the surface of the grains. 

Slip band extrusion occurs along these damaged planes and the presence 
of the holes may on occasion cause very fine extrusions in the form of 
filaments. 

The holes can be observed after electrolytic polishing or mechanical 
polishing, but electrolytic polishing also reveals the bands in the form of 
grooves. 

It has been ascertained that grooves may be produced by the presence 
of concentrations of solute atoms, and conversely zones depleted of 
solute atoms are electrolytically polished less rapidly than the matrix 
and appear as ridges. 

Although pure aluminium and also solution heat treated aluminium 
‘5% zine 2-5%, magnesium alloy show holes along the slip band traces, 
marked grooves are not observed. This suggests that the grooves are 
not primarily due to the presence of defects, but to variations in the 
distribution of the solute atoms. 
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ABSTRACT 


Measurements of the velocity and attenuation of plane sound waves 
in a number of organic liquids at pressures up to 10000 Ib/in.? have 
been made. From the velocity, the adiabatic compressibility is derived 
and from this, by comparison with the measured isothermal compres- 
sibility, the ratio of specific heats as a function of pressure is calculated. 

By introducing the measured attenuation into the Kirchhoff formula 
an ‘ultrasonic’ viscosity as a function of pressure is also derived. 


§ 1. INTRODUCTION 


MEASUREMENTS of the (isothermal) compressibility of liquids under 
pressure have in the recent past been obtained in the main in piezo-metric 
experiments and of their viscosity from transpiration experiments. 
Recent developments in ultrasonic techniques have made it possible to 
measure the velocity of sound with precision and so to derive the adiabatic 
compressibility. In conjunction with the aforesaid value of isothermal 
compressibility, it is then possible to derive the ratio of specific heats of 
a fluid. 

From measurements of the absorption of sound waves in the liquid 
it is possible, putting the coefficient into the Stokes—Kirchhoff formula, 
to find how the coefficient of (kinematic) viscosity varies with pressure, 
though it is true that the measured value of the latter coefficient may be 
a function of the frequency employed. The expression for the absorption 
coefficient («) of plane waves of frequency f in a liquid of kinematic 
viscosity v is 872fv/3c*, c being the velocity of sound in the liquid. 

The first measurements of ultrasonic velocity in liquids under pressure 
were made by Swanson (1933, 1934) using a Pierce interferometer at 
200 ke/s in nine organic liquids. Next Biquard (1939) used the optical 
diffraction method devised by himself and Lucas (Lucas and Biquard 
1932) in five liquids up to 650 atmospheres pressure (at 7 and 10 Mc/s). 
His accuracy was of the order of that reached by Swanson, the limit in 
both cases being governed by uncertain temperature fluctuations. 
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More recently a number of measurements of ultrasonic velocity to 
higher pressures have been made by several American investigators. 
As they confined their experiments to water they will not be discussed 
here. 

Only Biquard has reported absorption measurements (in toluene) but 
his optical diffraction technique, though capable of good precision in 
measuring ultrasonic wavelengths, is scarcely amenable to that deter- 
mination of amplitudes in different parts of the ultrasonic beam, necessary 
if absorption coefficients are to be deduced. 


§ 2. APPARATUS 


For the present experiments, both variable-path and fixed-path 
ultrasonic systems were built. It is common when working in liquids 
to send pulses, or short bursts of radiation, from the crystal source and 
to receive the echoes back from a reflector placed in the path of the wave. 
It is comparatively easy then to estimate the size of the echo as the path 
is lengthened at constant pressure or kept constant while the pressure 
is varied and so to measure the amplitude absorption coefficient. It is 


Fig. 1 
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Block diagram of ultrasonic pulse generator and receiver, showing wave- 
forms at three stages of modulation. 


less common to use the pulse technique to measure wavelength and so to 
derive velocity, but Mason (1950) has described a method which is well 
adapted to experiments under high pressure, and a circuit based on this 
principle was used to measure wavelength in the present work. 

A continuously running radio-frequency oscillator (fig. 1) supplies 
short pulses selected by a gating valve and applied to the transducer. 
Echoes from the reflector are received in the usual way but are then mixed 
with a set fraction of the original r.f. As the reflector is moved the phase 
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of the echo alters continuously, and as it moves in and out of phase with 
the original r.f., an echo ‘ beat’ effect is observed on the oscilloscope 
screen. One such ‘ beat ’ or pulse maximum occurs for each half wave- 
length traversed by the reflector. A large number of such wavelengths 
can be counted, the change in path length observed and from the wave- 
length and frequency the velocity determined. 

The technique thus combines the accuracy obtainable with a continuous 
wave method with all the advantages of the pulse method. An added 
advantage is that it can be used for very viscous liquids where high 
attenuation might make normal continuous wave methods impracticable. 

All measurements were made in the Fresnel zone of the transducer, 
i.e. at distances <R?/\, where R is the radius of the transducer and \ the 
ultrasonic wavelength in the liquid, and with a reflector of greater radius 
than R. A check with a moving reflector showed that the attenuation 
was linear with path length up to ?R?/A, after which the effects of 
beam divergence became apparent. Great care was exercised in setting 
the face of the reflector parallel to that of the crystal and in aligning 
both perpendicular to the path of the waves. 

In using the fixed-path apparatus for absorption measurements some 
difficulty was experienced. The obvious method would be to measure 
the change in echo-size as the pressure is varied, but the adoption of 
such an expedient would overlook the fact that the characteristic imped- 
ance of the liquid varies with pressure and so, therefore, does the load 
on the crystal and the intensity reflected at the solid—liquid boundaries. 

The arrangement finally adopted was to mount in the beam path an 
auxiliary reflector of small dimensions compared with the main reflector. 
A proportion of the beam was thus intercepted to provide a reference 
echo. Measurements were made by comparing echo levels from both 
reflectors and in this manner the main sources of error discussed above 
were eliminated. 

The acoustic system was mounted in a steel vessel capable of with- 
standing 1000 atm pressure. Details of this and its associated equipment 
are given elsewhere (Tait 1957). Liquid was fed in from a reservoir 
and, after closing a cock, pressure was applied by screwing up a piston 
into a compression cylinder, connected to the main vessel. The acoustic 
measuring vessel was surrounded by a liquid jacket with a thermostat. 
The thermostat was set in operation overnight which meant that the 
system was given more than 12 hours to attain temperature equilibrium 
before measurements were started. The range of pressure, measured on 
a Budenberg gauge (N.P.L. calibrated), was then covered both up and 
down. 

Some liquids were used in sufficient quantity (about 2 litres) to fill 
the whole system, including the compression cylinder, but in most cases 
the acoustic measurement chamber, which alone contained the liquid 
under test, was cut off from the compression cylinder and connecting 
pipes by a sylphon bellows, which communicated the stress applied to 
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water in the compression cylinder to the experimental liquid. Figure 2 
shows the acoustic chamber with the subsidiary reflector as mounted 
for the absorption measurements. 
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Fixed-path apparatus for absorption measurements; 8S, sylphon; R, main 
reflector ; A, auxiliary reflector; ©, crystal mounting. 


§ 3. RESULTS 

In this paper results are reported in: 

(A) carbon tetrachloride and benzene, typical of non-associated 
liquids showing ‘ anomalous absorption ’ ; 

(B) ether and pentane, liquids with a high compressibility at room 
temperature ; 

(C) carbon dioxide and monofluor-trichlormethane or freon 11, typical 
of fluids which are vapours at n.t.p. ; 

(D) glycerine, having a high shear viscosity. 
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(A) Carbon tetrachloride and pentane were obtained in Analar purity. 
Tables 1 and 2 give the velocity results at 3 Mc/s and three or four 
temperatures and on fig. 3 they are plotted as isotherms. Swanson’s 
points are rather scattered and bear a different slope to ours, but having 
regard to the difference in temperature there is excellent agreement 
between the present ones and Biquard’s. 


Fig. 3 
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Ratio of velocity V, to that at atmospheric pressure V,) in carbon tetrachloride 
and benzene, with water for comparison. 


Table 1. Velocity in Carbon Tetrachloride 


T=15°5"o Ti=25°o T=40°c 

P \b/in.? V m/sec P lb/in.? V m/sec P \b/in.? V m/sec 
0 950 0 908 0 872 
2280 1010 2000 962 3100 1058 
3900 1047 4100 1014 4000 1080 
5900 1090 5900 1054 6000 1024 
8000 1130 8100 1098 8000 1065 
10000 1165 10000 1132 9000 1083 


In comparison with results obtained for water (for which a mean curve 
is given) the two striking features of the curves for benzene and carbon 
tetrachloride are the greater slope and the curvature. Over the measured 
temperature range V/V, increases steadily with temperature. This is 
in accordance with the observed increase in compressibility with temper- 
ature for these liquids. 
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In table 3 is shown the adiabatic compressibilities 6 , which, together 
with the isothermal compressibilities By measured by Dolazalek and 
Scheidel (1920), enable values of y to be calculated. 


Table 2. Velocity in Benzene 


ARES IE) T'=25°C T =31-0°0 T'=40°C 
P \bjin2 V m/sec |P tbjin.2 V m/sec| P Ib/in.2 V m/sec | P lb/jin? V m/sec 
0 1326 0 1302 0 1276 0 1231 
520 1342 450 1318 530 1291 510 1245 


1030 1359 930 1332 1100 1307 940 1258 
1570 1377 1460 1349 1530 1323 1390 1273 
2130 1394 2010 1367 2060 1337 1870 1288 
2760 1412 2580 1384 2620 1355 2340 1303 
3400 1430 3180 1402 3210 1372 2860 1318 
4100 1449 3820 1420 3830 1390 3390 1335 
4820 1468 4510 1439 4490 1408 3960 1351 
5590 1488 5240 1458 5180 1426 4540 1468 
6400 1508 6020 1478 5770 1444 5170 1386 
7280 1529 6780 1498 6500 1464 5870 1403 
8200 1551 7660 1518 7330 1484 6570 1422 
Freezes 8530 1540 7800 1487 7320 1440 
9270 1557 8070 1504 8100 1459 

9000 1525 8860 1478 

9570 1539 9770 1498 


Table 3 
Carbon Tetrachloride at 15-5°c Benzene at 19°c 
Ee lbyia.2 6 el Ot bro 10% va x 108 fe, 1082 
em? dyne-! | cm? dyne-! Y em?/dyne | cm?/dyne ie 

0 74-2 104-4 1-41 67-8 96-4 1-42 
418 72-6 102 1-41 66-4 94-1 1-42 
1392 68:3 96-3 1-4] 63 87:8 1-40 
2781 63-2 88-7 1-40 58-6 81-6 1-40 
4180 58-8 81-5 1-39 55 75-4 1-37 
5560 54-5 75:4 1:38 51-9 70-2 1-35 
6960 51-3 69-2 1.35 48-9 65-1 1-34 
8350 48-6 63-9 1-32 46-5 61-9 1-33 


The values for y are given to the nearest three figures. It will be 
observed that y for both liquids decreases with pressure. For atmospheric 
pressure there are for comparison the values of y: 1:43 for CCl, and 
1-54 for C.H, given by Schulze (1925). The discrepancy is well accounted 
for by the uncertainty in the isothermal compressibility values. Reference 
to the International Critical Tables and other sources shows that the 
values of isothermal compressibility given by different observers do not 
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agree within 5%. The error in the absolute values of y may well be of 
the order of 3°, but the figures given in the table should represent the 
relative change in y to better than 1°%. 

Absorption in these two liquids was measured at 3 Mc/s and several 
temperatures. Table 4 gives the results in terms of the parameter 
a/f? and the kinematic viscosity corresponding to this frequency (vp) 
when «/f? is substituted in the Stokes—Kirchhoff formula. The fourth 
column gives the low-frequency or shear value of the kinematic viscosity 


Table 4 


Carbon Tetrachloride Benzene 


Pan |bin? x10: | aff? x 1015 veces mets Noll? sc101 rece, ye ce 


0 5-20 14-7 0-530 8-50 66-5 0-652 
2 4-56 16-4 0-636 7-70 70-7 0-737 
4 4-04 16-8 0-706 6-76 74-0 0-833 
6 3°64 Pied 0-805 6-00 76:2 0-937 
8 3°30 19-0 0-902 5°30 77-5 1-05 
10 3-10 20-2 1-020 4-80 78:8 1-18 
12 2-85 19-3 1-045 4-30 80-6 1-34 


(The units of kinematic viscosity in the tables are centistokes.) 


Table 5. Velocity in Ether in M/sec 


| 

P \bjin.? 16°c 20°C 30°C 44°C 
0 1023 976 945 862 
560 1047 1002 976 892 
1000 1071 1025 1003 921 
1500 1093 1051 1028 946 
2000 VV; 1073 1050 968 
2500 1138 1096 1073 992 
3000 (ibsyyi 1116 1094. 1013 
3500 1177 1137 1117 TLS 
4000 1198 1154 1137 1054 
4500 1217 1174 1156 1075 
5000 1233 1192 1176 1095 
5500 1253 1209 1192 1113 
6000 1268 L225 1207 1130 
6500 1283 1241 2 1148 
7000 1300 1256 1241 1164 
7500 1313 Ay [pe 1261 1180 
8000 1328 1287 127 1196 
8500 1342 1302 1287 1211 

9000 1358 1316 1300 12257 


+ Extrapolated value. 
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(v,) calculated from the viscosity and specific volume data of Bridgman 
(1931) and Schulze (1925). 

There are two contrary effects here, which one may expect to find in 
all liquids ; the velocity isothermals (fig. 3) curve towards the pressure 
axis whereas the isotherms of shear viscosity curve up, so that if both 
factors are inserted in the Stokes—Kirchhoff formula to derive what 
one might call the ‘classical’ absorption coefficient, a minimum will 
appear (at about 10000 lb/in.2) and thereafter the coefficient will rise. 

(B) Pentane and ether have the highest compressibilities among 
common organic fluids which exist in the liquid state at n.t.p. Tables 
5 and 6 show the velocity results at 3 Me/s. 


Table 6. Velocity in Pentane M/sec 


At 18°c both have the same velocity but that in ether is subject to 
the greater temperature coefficient. 

The isotherms are shown in fig. 4 and the compressibilities and values 
for y in table 7. 

To calculate the isothermal compressibility of ether there are Bridg- 
man’s values, but we could not find any for pentane. Accordingly the 
compression cylinder itself filled with pentane was used as a piezometer, 
i.e. p-v plots were obtained by finding the number of turns of the (cali- 
brated) screw driving the piston necessary to increase the pressure on the 
gauge by a determinate amount. This, of course, assumes that the 
change in volume of the thick steel vessel can be neglected in comparison 
with that of pentane over the range of pressure up to 10 000 Ib/in.. 
This was checked with—less compressible—water in the apparatus, 
when it was found that deviations from published data on the isothermal 
compressibility did not exceed 3°% at the highest pressures. 

As the absorption shown by these two liquids is much less than by those 
in group (A), it was necessary to measure it at a higher frequency, i.e. 
12 Mc/s. Furthermore it was necessary to make the measurements at 
room temperature or lower because of fluctuations in the signal at higher 
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Fig. 4 
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Ratio of velocity V,, to that at atmospheric pressure V, at 15°c in pentane, 
ether and freon 11, with water for comparison. 


Table 7 
Ether at 16°o Pentane at 15°c 

P \bjin.? 67 <10%* br LO y Baxclo® eare 102 y 

0 > HAS 170 1-34 151 190 1-26 

2000 100 130 1-30 127 153 1-23 
4000 87°5 112 1-28 115 136 1-175 
6000 79-5 98 1-24 96:5 108 1-125 

8000 71-5 87 1-22 87 93 1-07 

10000 65:5 78 1-21 80:5 847 1-045} 


+ Extrapolated value. 
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temperatures, attributed to scattering of the radiation by thermal 
inconsistencies such as occur near the critical state. Even at room 
temperature we must report some lack of reproducibility in the absorption 
measurements which reduces the accuracy of results below those of 


group (A). be 
A comparison of ultrasonic and shear viscosities appears in table 8. 
Table 8 
Ether at 16°c Pentane at 15°c 
P< 10-2) tad f* L022 vp CS V, CS p10 vps vcs 
0 0-60 1-66 0-233 1-00 2-20 0-205 
2 0-475 1-88 0-264 0:93 2-70 0-231 
4 0-398 1-92 0-296 0-865 3°20 0-264 
6 0-334 1-98 0-327 0-81 4-88 0-291 
8 0-288 2-03 0-358 0-77 4-53 0-325 
10 0-255 2-05 0-390 0-745 5:26 0-359 
950 
V 
(0 4 
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750 
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Velocity and attenuation in freon 11 at two temperatures. 

(C) Measurements in the vapour state have been made on carbon 
dioxide (critical temperature 31°C, critical pressure 76 atm.) by Parbrook 
and Richardson (1952) and in monofluor-trichlor methane or freon 11 
(critical temperature 198°C, critical pressure 43 atm.) by Miyahara and 
Richardson (1956). They, as liquids, are therefore in the same situation 


as the last two substances (group (B)) in that at room temperature they 
are highly compressible. 
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Table 9 gives the adiabatic and isothermal compressibilities of the 
liquids, the latter being derived for carbon dioxide from Bridgman’s 
measurements and for freon 11 from our own. 


Fig. 6 
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Velocity and attenuation in carbon dioxide at two temperatures. 


In the present research measurements have been made on the liquid 
state at one frequency (1 Me/s) and «/f? is plotted in the two figures 
(fig. 5 for CO, and fig. 6 for CCl,F) which comprise these results. It is 
characteristic of fluids at super-critical pressures that the velocity iso- 
therms are sharply curved. This is especially shown by these results 
in carbon dioxide, which carry the earlier measurements (Parbrook and 
Richardson 1952) to higher pressures. 
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Table 9 
Carbon dioxide at 19°o Freon 11 at 15°c 
Pb ba 10 by 0 y Plbjin* Vex 10%" | 8; < 10" y 
1000 55-5 130 2-3 0 113 167 =| 1-47 
1500 41-6 78 1-9 2000 92-5 130 1-41 
2000 30:8 55 1-75 4000 80-6 110 1-37 
2500 26-8 40 1-50 6000 71-5 95 1-33 
3000 21-7 32 1-35 | 10000 59-5 74 1-24 
Fig. 7 
230 
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Velocity in glycerine at six temperatures. 
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(D) Glycerine is typical of the rather few liquids which exhibit a 
large viscosity both to shear and to ultrasonic waves. The velocity 
isotherms are shown in fig. 7. Table 10 shows, as for the other liquids, 
values of adiabatic and isothermal elasticity. There appear to be no 
extant values for the latter, so we determined the volume: pressure 
relationships in our compression cylinder. So the y values were calcu- 
lated. 

Finally, table 11 exhibits the («//2) coefficients and the calculated 
‘ultrasonic viscosity ’. 

Bridgman has shear viscosities for comparison. It will be observed 
that the ‘ ultrasonic viscosity ’ is of the same order as that derived from 
shear experiments. 


Table 10. Glycerine at 23°c 


P |bjin.? 0 2000 4000 6000 | s000 10 000 
Ba x10? 21-95 | 20-9 20-4 19-9 19-6 19-4 
pe 10M 28-0 25-5 24-0 23-0 22-5 21-0 

y 1-275 1-215 1-18 1-155 1-125 1-085 
Table 11 
PE ibrn* 10- 4 a ffx 10 v, Stokes v, Stokes 
0 13-00 3-42 3-02 
2 13-72 3:90 3°40 
4 14-45 4-38 3°82 
6 16-20 5-05 4-38 
8 15-90 5-70 4-95 
10 16-60 6-72 5-84 


§ 4. Discussion OF RESULTS 


It is the purpose of this work to provide experimental data which might 
form a basis for theoretical discussion rather than to propound theories, 
to which indeed the velocity results do not lend themselves. No tabula- 
tion of absorption coefficients would however be complete without some 
reference to relaxation processes. It is well known that the postulate 
of a mean relaxation time for the molecules varying directly with pressure 
explains well the propagation of ultrasonic waves through gases, remote 
from their critical states, but this idea meets great difficulties in applica- 
tion to liquids because (a) no dispersion of velocity with frequency has 
been found, (b) the absorption per wavelength in most cases changes 
steadily with frequency in most liquids and does not show the peak 
required by relaxation theory. It is true that most of our attenuation 
measurements have been made at one frequency but the steady change 
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of the coefficient with pressure and its nearly constant ratio to shear 
‘viscosity gives no indication of relaxation within the range of frequeney 


and pressure employed. 

Of course, there remains to be explained the very much greater values 
of the deduced ‘ ultrasonic viscosity ’ in comparison to those in shear, 
except in glycerine. To justify a suggestion, on these grounds, that the 
Stokes—Kirchhoff formula is generally invalid and that its applicability 
to the case of glycerine is fortuitous would take us outside the scope of 
this paper, but the invocation of a bulk viscosity several hundred times 
that of ‘the shear viscosity, which some have favoured, does not really 
help the understanding of this phenomenon. 
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ABSTRACT 


Copper, silver and gold have been evaporated onto heated cleavage 
surfaces of alkali halide crystals in an electron diffraction camera. The 
diffraction patterns have been studied continuously as the depositions 
were carried out and the crystalline arrangement in the films has been 
deduced as a function of thickness. When the growth is epitaxial the 
orientations are predominantly parallel, with subsidiary twinning and 
the (112) orientation. The proportion of the latter orientations, the 
crystal size and the degree of disorientation at different levels in the 
films may vary considerably. Single crystal growth by evaporation is 
found in certain cases to cease beyond a certain thickness ; an explanation 
of this is put forward. Fine structure effects have been observed in the 
patterns. These effects have been analysed in some detail, and possible 
explanations are suggested. 


§ 1. IyTRODUCTION 


Ir is now well known that a metal layer condensed from the vapour onto 
a single crystal surface may, under certain conditions, crystallize epi- 
taxially. Study of this phenomenon originated with the now classic 
experiment of Lassen and Bruck (1935) who deposited silver onto cleavage 
surfaces of rocksalt maintained at various temperatures. For substrate 
temperatures of over 100°c the silver was found to give an electron 
diffraction pattern characteristic of a single crystal oriented with its 
cubic axes parallel to those of the rocksalt. Since then the great deal of 
work which has been carried out on the epitaxy of metallic layers has 
been largely dominated by attempts to explain the observed orientations 
in terms of the geometrical theories of Royer (1928, 1954) and Menzer 
(1938) and the more recent ideas of Engel (1953) and of Frank and van 
der Merwe (1948a,b). These seek to explain orientation solely in terms 
of the relative atomic spacings and arrangements in the contact planes of 
the substrate and the deposit. No other line of approach appears to have 
been considered theoretically. The experimental results on epitaxy 
have led Neuhaus (1952) to suggest that the hardness of the substrate may 
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be important, and very recently Gottsche (1956) has shown that the sub- 
strate hardness can be related to the degree of epitaxy in deposits of 
copper, silver, gold and aluminium condensed onto four different alkali 
halide crystals. 

In a study of twinning in oriented deposits of silver on rocksalt (hence- 
forth denoted thus: Ag/NaCl) Kirchner and Cramer (1938) found that 
the crystallographic structure in the films varied with the mean thick- 
ness. In particular, diffraction spots from twinning on {111} planes, 
not observed from films only 8 A thick, were strong in the patterns from 
films of thickness 200A but were again weak in films of 10004 thick- 
ness. Such a result makes it clear that the greatest care must be exer- 
cised in any discussion on the structure of condensed layers. The true 
orientation of the deposit on the substrate can only be inferred from 
deposits a few angstroms thick. Once the substrate is covered with the 
deposit material it will cease to have any direct effect on further deposits 
which can, and in many cases do, crystallize in a different way. The 
mode of crystallization of the continued deposit will, at any thickness, 
depend primarily upon the crystalline nature of the surface, on the 
temperature of that surface, on the rate of evaporation and on the residual 
gas pressure in the apparatus. The structural changes which occur in 
evaporated deposits at different stages during their growth have been the 
subject of very little experimental work. Uyeda (1942) has studied 
Ag/NaCl at 350°c and also the evaporation of silver onto a number of 
other substrates by carrying out the evaporation in the electron diffraction 
camera. This work was not extensive, the results on Ag/NaCl, for example, 
being little more detailed than those of Kirchner and Cramer (1938). 
Recently much work has been published on the structure and topography 
of oriented films as shown by combined electron microscopy and diffraction 
(e.g. Yoshida 1956, Fischer and Richter 1955), but the films used have 
been comparatively thick (10004) and have not been studied during 
their formation. Trillat et al. (1956) have used a novel technique in the 
study of a Au/NaCl, 400°c deposit. Bombardment of such a deposit by 
high energy ions in the electron diffraction camera causes the gold to 
evaporate and thus the deeper layers successively produce diffraction 
until the substrate pattern appears. It was found that the initial 1000 & 
deposit gives a Debye-Scherrer pattern but that when a certain amount 
of gold is removed, an oriented single crystal pattern of gold appears and 
persists until the rocksalt pattern shows through. The catastrophic 
effect which the ion beam has on the gold must, however, preclude any 
detailed deductions on the structure of deposited layers. 

In the work described here, the metals copper, silver and gold have 
been evaporated onto heated alkali halide cleavage surfaces under care- 
fully controlled conditions. The evaporations were carried out in an 
electron diffraction camera and so the crystalline nature of the films 
could be studied continuously from the thinnest detectable deposits of 
mean thickness less than 1 A to the maximum thickness of about 200 4, 
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§ 2, APPARATUS 


The details of the specimen chamber, which incorporates a specimen 
furnace and an evaporator, have been described elsewhere (Kehoe et al. 
1954,1956a). The diffraction patterns were recorded on film, thus allowing 
a large number of exposures to be made during the course of an experiment. 

The rate of evaporation from a bead of molten metal is approximately 
constant and it is therefore possible to estimate the deposit thickness at 
any stage of an experiment from a knowledge of the time taken for the 
evaporation of a known weight of the metal. Thicknesses found in this 
way were in tolerable agreement with estimates based on the thickness at 
which the Kikuchi pattern from the substrate became obscured. The 
latter thickness has been determined for copper by Newman and Pashley 
(1955) with the aid of a radio-active isotope. The rate of deposition was 
always 3-4A per minute. Each experiment lasted about an hour and 
during this period the temperature of the specimen was held constant to 
within a few degrees. 


§ 3. THE STRUCTURE IN EVAPORATED FILMS 


In nearly all the orientations observed in this work the plane of the 
metal deposit parallel to the (001) cleavage surface of the substrate 
contains a [110] direction oriented parallel to the [110] directions of the 
substrate surface. For such orientations it is necessary to name only the 
deposit plane which is parallel to the substrate surface. 

In general, the diffraction pattern of a deposit consists of spots from 
one or more orientations superimposed on Debye-Scherrer rings from that 
portion of the substance which has crystallized at random. The relative 
intensities of the spot and ring patterns give an indication of the pro- 
portions of oriented and disoriented material. When a surface is studied 
in reflection, only the uppermost 104 or so of the specimen contribute to 
the pattern and the proportions estimated refer to this thickness only. 


3.1. The Evaporation of Copper 

The principal substrate was rocksalt and the first results given are of 
the Cu/NaCl experiments. Estimates of the orientation proportions are 
given in table 1 for a series of substrate temperatures. The character of 
the crystallization at some temperatures depends markedly on the 
deposit thickness. At 160°c and below the orientation tends to improve 
with thickness but at the highest temperatures random crystallization 
occurs on initially perfect orientation. 

In a few experiments with each of the three metals an interesting, 
transitory feature of the diffraction patterns was noticed. The pattern 
from the substrates comprised spots and Kikuchi lines and bands. Im- 
mediately the evaporation was started, very faint diffuse haloes appeared 
in addition to the much stronger substrate pattern. The visibility of the 
haloes was poor and in many experiments it was difficult to be sure of 
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their presence. At a mean deposit thickness of about 2 A, the haloes 
could no longer be distinguished, the diffraction rings or spots from the 
deposit having developed and increased the background intensity of the 
pattern. It must be presumed that these haloes resulted from amorphous 
or liquid-like aggregates of deposit atoms on the flat substrate surface. 
The most prominent halo lay for each metal between the (111) and (002) 
Debye-Scherrer ring positions. 


Table 1. Orientation proportions in the copper/rocksalt films 


Mean deposit thickness 
Substrate 
temperature (°C) az ea 200 & 
20 O O O 
90 O D D 
130 C B B 
160 B B A 
195 A A A 
230 A A A 
265 A A A 
315 A A B 
345 A A B 
385 A A B 
Key: A: completely oriented. 
B: mostly oriented. 
C : oriented and random crystal comparable. 
D: mostly random. 
O: completely random. 


The spot pattern from the substrate became obscured at a deposit 
thickness of about 6 A, while the Kikuchi pattern persisted until the thick- 
ness was about double this value. 

Orientation was first observed at 90°C and the growth characteristics 
were almost identical at 130°c. At deposit thicknesses of less than 12 A 
the orientation was solely parallel, i.e. the cubic axes of the copper were 
parallel to those of the rocksalt. The mean crystallite size was about 30 A. 
At 124 mean thickness two further orientations occurred to a small 
extent and continued in the same proportion as the deposit was further 
built up. These were the (112) and (221) orientations, the latter arising 
from twinning on the {111} planes of the parallel crystallites. 

At substrate temperatures of 160°c and over, the initial, and dominant, 
orientation was again parallel, with the twin and (112) orientations 
appearing at 12 A. Figures 1 and 2 show the diffraction patterns from 
deposits of 2A and 20 A respectively on rocksalt at 160°c. The mean 
crystallite size was about 504 at this thickness. The extent of the 
twinning in the deposit increased with thickness and reached a maximum 
at 50 A to 804 and double diffraction between the parallel and twin 
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orientations was then considerable. The tendency to twin decreased as 
the thickness went beyond 1004. While the twinning was at its maximum, 
the strong diffraction spots of the parallel orientation became slightly 
streaked towards the shadow edge. The behaviour after this stage 
differed in the two temperature ranges 160-265°c and 315-385°c. In the 
lower range the streaks in the diffraction patterns became very strong 
and almost masked their parent spots. This is clearly illustrated in fig. 3, 
the pattern from 160A of copper at 160°c. The streaks arise from re- 
fraction of the spots and show the copper surface to be smooth. Kikuchi 
bands from the copper were usually seen at this stage, showing a consider- 
able degree of perfection in the deposit. Faint streaks were observed 
along the (111) directions of patterns, due to the presence of {111} stacking 
faults in the parallel orientation. Furthermore, the spots of the twin 
orientations were streaked along their (111) directions, again an indication 
of stacking faults. This can be seen in fig. 3. 

In the temperature range 315-385°c, the refraction streaks which 
appeared at about 60 A were sharper and did not lengthen as the evapor- 
ation continued. In the main, the patterns retained their cross grating 
nature, indicating a rough surface. At about 1004 thickness random 
crystallization started, the proportion increasing with thickness. 

A few deposition experiments were carried out with potassium chloride, 
bromide and iodide substrates. Below 220°c for potassium chloride and 
below 150°c for the other two salts the structure of the deposited copper 
was, thickness for thickness, exactly like that on rocksalt below 160°c. 
Above these temperatures the similarity was with Cu/NaCl above 300°c 
in that the surfaces were never observed to become smooth and random 
crystallization occurred when the thickness reached 100 A. 

Fine structure effects were observed in some Cu/KI experiments ; they 
are described and analysed in § 4. 

In diffraction patterns from the thinnest layers (less than about 6 A) 
the sharp streaks and spots from the substrate still showed, and provided 
a standard spacing with which to compare the deposit spacing. This 
could be done with an accuracy of about 1°% and to within this figure the 
copper spacing was found in all experiments to have its normal value. 


3.2. The Evaporation of Silver 

Silver has been condensed on to rocksalt and potassium bromide with 
almost identical results. The relative proportions of oriented and dis- 
oriented crystallites have been estimated in table 2. At all temperatures 
there was a tendency towards disorientation as the films were built up. 

At 50°c and 90° the initial orientation was solely parallel, but at 104 
the twin and (112) orientations appeared to a small extent. The crystal 
size at this stage was about 504. There were no further changes in the 
character of the deposits as more silver was deposited. The mode of 
growth of the silver is clearly similar to Cu/NaCl up to 130°o) @ihe 
similarity with copper growth also holds at higher substrate temperatures. 


Z212 
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Ag/NaCl and Ag/KBr at 145°c, 180°c and 215°C correspond to Cu/NaCl in 
the range 160-265°c. A minor difference was that when the twin and 
(112) orientations appeared, the (111), (115) and 

Ag (112) [110] // NaCl (001) [100] 
orientations were also observed. 

Finally silver deposition at temperatures over 215°c was similar to 
Cu/NaCl 315-385°c. At no temperature were Kikuchi bands observed 
from silver deposits. No changes in the silver lattice spacing were 
observed. 


Table 2. Orientation proportions in silver layers 


Bi betrate Mean deposit thickness 


temperature (°C) Ti Wk 100 & 


20 
50 
90 


215 


Key: Ne KBr Symbols as in table 1. 


3.3. The Evaporation of Gold 


Gold was deposited only on to rocksalt. The proportion of orientation, 
estimated in table 3, again tended to decrease with thickness. 

Orientation was first observed at 130°c and up to 230°C was of a type 
not found with copper or silver. The thinnest deposits gave a diffuse 
halo pattern which resolved into rings and spots from the parallel and (111) 
orientations in comparable proportions. At about 104, the twin, (112), 
(115) and other unidentified orientations began to appear. As the 
deposition continued the intensity in the diffraction pattern arising from 
oriented gold decreased until, at 1504, the only spots remaining were 
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those of the (111) orientations on the centre line. Crystallization had 
thus become random except for a slight tendency for the close packed 
(111) planes to form parallel to the substrate surface. 

With the substrate at 260°c, the sequence was the same until the thick- 
ness of the gold reached 40 4, when the parallel orientation instead of 
dying out became more prominent at the expense of the (111) and other 
orientations. This process reached completion at about 804A when only 
the parallel and twin orientations remained. At 1504 the gold surface 
tended to become smooth. 


Table 3. Orientation proportions in gold/rocksalt films 


Re hetrate Mean deposit thickness 
temperature (°c) Ti 10% 100 & 
20 O O O 
85 O O O 
130 O C D 
160 C C D 
190 C C D 
230 D C C 
260 O C C 
295 B A see text 
ole A A see text 
330 A A see text 
385 A A see text 


Symbols as in table 1. 


At 295°c the change in crystallization sequence reached its final form. 
and did not change further at higher temperatures. The initial orien- 
tation was solely parallel with crystallite size about 504A. Refraction 
streaks appeared at 204. A crystal size increase occurred when twinning 
commenced at 404 thickness. Between 804A and 100A thickness, the 
gold surface became very smooth and also Kikuchi bands appeared. At 
about 1204, a tendency of the gold to crystallize at random became 
apparent and at 150 A the proportion of disoriented gold was appreciable. 
The ratio of intensity in the spot pattern to that in the ring pattern 
depended markedly on the incident angle of the electron beam. This is 
seen in figs. 4 and 5, which are patterns from a 150A deposit. At a 
glancing angle of 0-7° the contribution of oriented crystallites to the 
pattern is negligible, whereas at 1-7° the contributions from oriented and 
random gold are comparable. The situation is qualitatively the same 
after a further 25 A deposit but the proportion of spots for a given angle 
of the electron beam is even less. : 

No changes of lattice parameter were detected in the initial gold 


deposits. 
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§ 4, SupsipraRY DirrRacTION EFFECTS 

In $3, the main crystallographic features of the metal deposits have 
been described. A number of secondary effects observed in the diffraction 
patterns and which do not permit of such definite interpretations form 
the subject of this section. 

In patterns from smooth copper surfaces, Cu/NaCl 160-265°c, extra 
streaks were observed midway between the normal spots at the [100] and 
[010] azimuths. These might have been due to adsorbed monolayers of gas. 
Oxygen monolayers on copper cube faces have been studied by Schlier 
and Farnsworth (1954) by the diffraction of slow electrons. They found 
two stable arrangements, the surface density of oxygen in one being twice 
that in the other. In the present work such monolayers would only be 
detectable on a smooth surface. The more densely packed oxygen layer 
would produce diffraction streaks coincident with the copper refraction 
streaks. However, the other type of layer would produce additional 
streaks midway between those from the copper at the [100] and [010] 
azimuths. It seems likely, therefore, that the streaks arose from an 
adsorbed monolayer of oxygen, but it must be noted that a similar arrange- 
ment of any other atoms, e.g. nitrogen, would produce similar continuous 
streaks. They are not consistent with a distortion of the copper atoms 
on the (001) surface having the periodicity of the unit cell. 


Fig. 6 


e 
(a) () 
Patterns from 200 & of copper on potassium iodide at 220°c. 
(a) [100] azimuth. (6) [110] azimuth. 


The fine structure effects mentioned in connection with Cu/KXI took the 
form of fine lines (* spines ’) and very small spots or short dashes (‘ specks’). 
The appearance of spines and specks can be interpreted unambiguously 
in terms of thin filaments of intensity in reciprocal space (Kehoe et al 
1956 b). The relevant features of the patterns at the [100] and [110] 
azimuths of Cu/KI 220°c are shown diagrammatically in figs 6.) he 
spines at the [100] azimuth are consistent with filaments passing through 
all reciprocal lattice points in the following directions : [102], [102], (012] 
[012], [101], [101], [011] and [011]. These bundles of filaments Drorvice the 
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spines at the [100] azimuth and some of the specks which occur are asso- 
ciated with them. The above filaments give neither spines nor specks at 
the [110] azimuth. The remaining specks at the [100] azimuth and those 
at the [110] azimuth are consistent with (111) bundles of filaments 
passing not through the integral reciprocal lattice points but through 
0, 3, 3 and all equivalent points. The expected hexagonal symmetry of 
these bundles would be achieved if extra parallel filaments passed through 
the integral points in the ¢{111) directions. Such filaments occur as a 
result of stacking faults, in any case, and would mask spines in these 
positions. 

A slightly different type of fine structure was observed in patterns from 
Ag/NaCl 340°c and Ag/K Br 340°c. Short, sharp spines appeared at the 
[100] azimuth passing through spots in the [010], [011] and [011] directions. 
No specks were observed. The intensity in reciprocal space was examined 
by making exposures every 2° between the [100] and [110] azimuths. 
The short spines were found not to disappear off azimuth but continuously 
to change direction and finally to merge with the stacking fault streaks at 
the [110]azimuth. They therefore arise from thin {111} sheets of intensity 
surrounding each reciprocal lattice point. 

The spots and refraction streaks in the patterns from smooth silver and 
gold surfaces always split into two components separated by twice the 
distance of the Ewald sphere from the expected positions of the intensity 
in reciprocal space. The effect was not confined to the parallel orientation 
but was observed with spots of the (112) and (111) orientations also. No 
interpretation of this effect can at present be offered. 

In the patterns from some smooth copper and gold deposits there 
occurred streaks, normal to the shadow edge, and bearing no simple 
relation to the metal lattice spacing. Moreover, their general appearance 
was quite unlike that of refraction streaks or spines. It was concluded 
that they were produced by an oriented impurity, in the form of either a 
monolayer or of crystallites giving diffraction only by double scattering 
with the metal. 


§ 5. DISCUSSION 

The results in §3 show clearly that in the most important aspects the 
crystalline arrangements in layers of each of the three metals on each of 
the substrates are similar. The crystallization is random at low substrate 
temperatures, and at higher temperatures the orientation is predominantly 
parallel. The twin and (112) orientations occur almost invariably at a 
mean thickness of 12 4. The quality of the deposited crystal generally 
improves with thickness until a thickness of 60-100 A is reached when the 
general tendency to disorientation manifests itself. 

There is so far insufficient experimental data to permit of a definite 
interpretation of the diffuse haloes which were observed in a number of 
experiments during the initial stages of deposition. It has been shown 
that during the evaporation of copper onto a flat (111) silver surface the 
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deposited atoms can migrate over distances of the order of I (Kehoe ef 
al. 1956 b). It is probable that the migrating units are not only single 
atoms but also liquid-like groups of atoms which crystallize only when 
sufficient energy is lost to the substrate and a suitable site is found. With 
the initial conditions of a clean flat surface and a sufficient density of 
migrating groups, the latter would be expected to contribute to the dif- 
fraction pattern, and it seems likely that this is the origin of the haloes. 
- The crystallites which thus form at this initial stage are widely separated. 
on the substrate surface and are therefore influenced by only the sub- 
strate during crystallization. It is only for such deposits that epitaxy 
should properly be discussed, as Gottsche (1956) has recently pointed out. 
In his experiments he studied the orientation of deposits sufficiently thin 
for the crystallites to be widely separated. For each metal-substrate 
pair two substrate temperatures can be defined: 7'y, below which the 
deposit is completely random, and 7’,,;,, above which orientation is 
complete. The values obtained for 7’, and 7',,,, in the present work are 
given for rocksalt substrates in table 4. For comparison, Gottsche’s 
corresponding results are given in parentheses. The two sets of results 
show the same trends but naturally differ in details owing to different 
evaporation rates and substrate quality. 


Table 4. Limiting Epitaxial Temperatures 


Ty (°c) PE soins (°c) 
Ag/NaCl 20-50 (100) 50-90 (150) 
Cu/NaCl 90-130 (100) 160-195 (150) 
Au/NaCl 130-160 (120) 285-305 (400) 


The departure from solely parallel orientation occurs when the mean 
thickness is about 12 A, i.e. when the surface coverage is 20-40°,. The 
appearance of the twinning and of the (112) orientation marks the end of 
the ideal epitaxy experiment. The proximity of crystallites on the surface 
is such that strains develop, to be relieved only by the formation of sub- 
sidiary orientations. Twinning is brought about by stacking faults on 
{111} type planes and it might be expected that stacking faults in the 
twins themselves would occur, and produce further derived orientations 
which may be termed * secondary twins’. However, there is no evidence 
of any such orientations in the above or any other work. The (112) 
orientation does not seem to have any simple structural relationship with 
the parallel orientation. 

The size of crystallites in layers thinner than 50 4, and in many cases 
thicker, lies between 40 A and 604. Keith (1956) has studied by x-rays 
the crystallite size in copper deposited on to surfaces at —183°c and sub- 
sequently annealed at room temperature. In high vacuum (better than 
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10-7 mm Hg) the original 40 4 crystallites grew to 400A by this process, 
but in the presence of oxygen the upper limit obtained by annealing was 
only 604 owing to an adsorbed oxygen layer on the copper surfaces. 
This seems to be in agreement with the present work in which the 
residual pressure in the apparatus was relatively high (between 10-3 and 
10-4 mm Hg). 

The respective temperature ranges for copper, silver and gold, in which 
smooth surfaces form, are probably those in which the mobility of the 
deposited atoms allows cracks and strains between adjacent crystallites 
to be coherently covered. The subsidiary orientations then weaken as 
uniform parallel crystals of large lateral size are formed. The relative 
perfection of the crystals is also indicated by Kikuchi patterns in the cases 
of copper and gold deposits. The appearance of disoriented crystallization 
in hitherto perfect deposits of copper and gold shows that crystal growth by 
evaporation on to a smooth (001) surface will not continue indefinitely, 
probably a direct consequence of the high surface mobility. While the 
crystal surface is rough there is an abundance of sites for parallel crystalli- 
zation, but once a smooth surface has developed, sites will no longer be 
present to the same extent and consequently migration distances will be 
much greater. If, for any reason, disoriented patches should be present, 
they would provide sites and would grow more quickly than the regions 
of smooth parallel crystal. The slightly raised polycrystalline islands, 
increasing their surface coverage as the deposition is continued, would 
produce exactly the effects found with gold deposits in §3. Trillat et al. 
(1956) observed that both parallel and polycrystalline gold existed in 
layers on rocksalt at 400°c, but their technique did not allow any further 
investigation of the phenomenon. 

The high sensitivity of the structure in the evaporated films to the sub- 
strate temperature suggests that there is no large temperature increase 
at the surface of the deposit. Rises of a few hundred degrees which are 
claimed to occur in certain instances (Murbach and Wilman 1953, Wilman 
1955) must be ruled out here. 

Spines and specks in a diffraction pattern have for the first time been 
fully interpreted in terms of intensity in reciprocal space (§ 4). Elementary 
diffraction theory shows that sharp, uniform filaments of intensity can 
result either from extensive plane monolayers of atoms or from absent 
monolayers in an otherwise perfect crystal. Similarly, sharp sheets of 
intensity which were observed in reciprocal space from Ag/NaCl —340°c 
deposits can arise from extra or absent lines of atoms in the silver. The 
existence of dislocations in the deposits is an obvious possible explanation. 
The diffraction patterns of two dislocation models have been calculated 
(Wilson 1949, Suzuki and Willis 1956) and that of Wilson may be relevant 
to the effects observed from Ag/NaCl. He has shown that screw dislo- 
cations produce round the reciprocal lattice points discs of intensity 
perpendicular to the dislocation axis. Further work would be needed to 
compare the experimental work with the theory in full. 
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ABSTRACT 


An effective range theory of low energy meson—nucleon scattering is 
developed based on the relativistic pseudo-scalar interaction. An 
approximate relativistic equation is derived which, unlike the ladder 
approximation to the Bethe-Salpeter equation, preserves crossing 
symmetry. The equation is solved approximately and all divergences 
are removed by renormalization before going to the low energy limit. The 
resulting equations for the phase shifts depend on only one parameter— 
the coupling constant. With a suitable choice of this constant the p-wave 
phases, in particular that of the 3, 3 phase, are in good agreement with 
experiment. The s-wave phases are of the right order of magnitude, but 
do not show the observed strong dependence on isotopic spin. 


§1. INTRODUCTION 


In a previous paper (Edwards and Matthews 1957, to be referred to as I) 
a conceptually extremely simple derivation was given of the Chew and 
Low (1956) no-recoil effective range formula for meson—nucleon scattering. 
This derivation was based on the idea of an effective potential arising 
from pseudo-vector coupling. According to the effective range theory, 
the experimental points, appropriately plotted, lie on a straight line. 
The theory contains two positive parameters (the coupling constant g? 
and the cut off) and both are used in fitting this line. Thus apart from 
the general feature of leading to an effective range formula, only the sign 
of the slope, or physically the fact that the resonance occurs in the 
J—=3 T=3 state, is specifically predicted by the theory of paper I. 
The coupling constant required agrees with that determined by photo- 
production. 

In this paper we give a relativistic generalization of I, starting from the 
renormalizable y,; coupling. We reproduce effective range formulae, 
but the theory is renormalized before the low energy limit is taken, so 
that only one parameter (the coupling constant) remains. It is shown that 
it is still possible to obtain a very reasonable fit to the experimental data 
on the (33) phase shift with qualitatively correct results for the (11) and 
(13) phase shifts. The coupling constant required remains at the value 
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which agrees with the photoproduction data. This theory also gives 
s-wave phase shifts of the correct magnitude and some isotopic spin 
dependence observed experimentally. This is discussed further below. 

To derive these results we put as a condition on our relativistic theory, 
that it reduces simply to the non-relativistic theory, and the approximate 
solution must reproduce the well known exact solution in the no-recoil 
limit for pair theory. It is important to note that this condition is not 
satisfied by the ladder approximation to the Bethe—Salpeter equation. 
This gives the Tamm—Dancoff approximation in its non-relativistic limit, 
which, as explained in I, excludes virtual meson pair creation, and actually 
replaces a trivial algebraic equation by an integral equation. 

In § 2 we derive a new relativistic equation for the scattering 7'-matrix, 
which does not suffer from these defects. This equation is solved 
approximately in §3 for the low energy p-wave phase shifts and in § 4 
for the s-wave phase shifts. 


§2. RELATIVISTIC THEORY 
The Lagrangian L is given by 
L=48, $09 Pat eee) +Tabat Hiv :8u—m+ifOnba)p, — - (2-1) 
where (), is the operator which defines the interaction, J, is an external 
source, and ¢, includes a quantized field and a classical unquantized field 


dy. Then the equation for the complete one nucleon propagator is 
(suppressing all space variables) 


1420 
| $2+00a($.-157-) |] 9HI=3, 22) 
where x 
S~t=y,0,--t1(m—te). . . 0) «02 (2:38) 
If 4, is the complete meson propagator 
we ) 
Ba ca here (2.4) 


Iterating eqn. (2.2) once, gives 


S19($)S-1=S-—| $0, (fart Soe 5 )—F°0a (Hart 


o¢, 
5 5 
The T-matrix for meson-nucleon scattering is determined by 
AC) 
Yad BS -1 
[eee ees 


Carrying out the differentiation on (2.5) we get 


D3 Faete Vt A pST s+ Ve Age, 
Ohad ae 


13 S414 om Sree e . . . ° ° (2.7) 
l m 
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where Petes (oO SO.) 5 > Se ye aes (2:8) 


Note that V, S, and 7' are matrices both in y-space and in the position 
of the nucleon, and matrix multiplication is implied. The suffices given 
explicitly carry both the charge label and the position of the nucleon. 

This equation is exact. The approximation of taking only the first 
two terms on the right-hand side is precisely equivalent to the ladder 
approximation to the Bethe-Salpeter equation. This is very artificial 
from our present point of view, and the approximation we shall adopt is 
to take the first three terms. This obviously preserves the crossing 
symmetry, and the simplest technique of solution discussed below when 
applied to pair theory can be shown to reproduce the correct equation in 
the no-recoil limit. (It is true, in general, that the final term is pre- 
dominantly a self energy effect in the limit of low energies and consists of 
a multiple of the first three terms which is absorbed in renormalization. 
Thus the condition mentioned above is satisfied. The proposed equation 
is shown graphically in fig. 1. 


Fig. 1 


Be er 
es 


Written in full with nucleons and mesons having incoming four vectors 
R, P and outgoing S$, Q respectively, and meson charge indices 7, s, 
7, (P,Q; R, S)=V,,(P,@; R,$ +] Vee ( (P, P-+X; R, R—X)S(X-+-P) 
x A(R—X)7',, (P+ X, @; R—X, S)d*X 
+ [Var (—@, X—@; R, R—X)S(X—@)A(R—X) 
xT ,,(xX—@, —P; —X-+R, S)d*X, . (2.9) 
where P+Q-+-+-R+ S=0 and 
rs (P, @; R, S)=f?0,S(P+R) 0, tf?0,S(R—@)0,.  . (2.10) 
Spinor indices are not shown. We shall be interested in the value of 7’ 
on the energy shell in the centre of mass system, i.e. 


P=(ppo, p)=[(p? +)”, Pl, | 
eda) el ees (2nd 
R= hor) [(r2-+-m?) yV/2, r], 

Pau 


time)", s) 
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§3. SoLuTION oF THE Equation. THE P-WAVES 


We now specify the interaction to be 

Olt LS ee ee (3.1) 
The non-relativistic equation, of which (2.9) is a generalization can be 
solved exactly, since in momentum space it is algebraic. Another point 
of view is that if the Fredholm method is applied to the non-relativistic 
equation, the resolvent and determinant terminate with only one and 
two terms respectively. Equation (2.9) is much more complicated, but 
we shall obtain an approximate solution by means which are exact in 
the non-relativistic limit. There it is found that if the equation is iterated 
to give 


[Pes SGI. SIT a) a ee eer an ea 
and after angular momentum and isotopic spin projection operators are 
applied, the terms are recombined in the form 

S,) 
SC Ke 
this form is indeed the exact solution. So our procedure is to evaluate 
the relativistic form of (3.2), and having projected from four spinors to 
two spinors by applying (1, op/2m) and (1, —oq/2m) to left and right 


t (3.3) 


Fig. 2 


respectively, and separating into the states of angular and isotopic spin 
by the well known projection operators P,,, P33 and P,, (as in I), the final 
answer is again expressed in the form (3.3). This is equivalent to taking 
the first approximation of the Fredholm solution, appropriately generalized 
to eqn. (2.9). Our final formula has also been used by Sartori and 
Wataghin (1954) in their application of the Cini and Fubini (1954) 
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variational principle to the non-relativistic problem, as has been mentioned 
in I, so our solution below could also have been obtained by a relativistic 
treatment of the Cini-Fubini procedure. In other problems such as the 
two nucleon system, the two approaches would appear to be quite 
different however. 

The approximate solution of (2.9) can therefore be constructed from a 
knowledge of the graphs of fig. 2 and those obtained from them by 


crossing. Their values, retaining only p-wave contributions, are, after 
renormalization, 


(x)= —1(2:)4f?7,7,0 . po . q/4m2a, ] 

(a) =—Sf'7,7,(0?)2mi)o . po .q/4m?, | 

(b)=(c) Se eae) 
| 


= —f*tr,7,(7?/2mi)o . p. o . q/4m?, 
(4)=—3/4[48,,—7,7, l(a?/2mi)p . q/4m?. 
The renormalization of the vortex part, is done so that terms (b) and 
(c) contain no terms in 1/w, and thus the coupling constant is determined 


by the coefficient of the neutron pole at w=0 in Born approximation. 
Adding the crossed graphs, we find 


Sy?) = — (2m) f%(p?/4m? ei) Pa, 


Oe af An on) pee jy Pe | ee ee) weet (S.0) 

where 
Vipers Bale) edie Ah (3.0) 
N,=—13, —7, —10, OW tee Cr wl ose) 


for ~<=11, 13 and 33 respectively. 
As in I, the phase shifts are now given by 


Ra COO 7 t= 1 Co) hy), <o 2t a Fin a pee (3.8) 
Av=HG7/4)tg se ee ee (8.9) 
Wy=(40?r,/mg}), . . . . . . (3.10) 
where 
ON at et ee F841) 
and g is the pseudo-vector coupling constant 
U(r. ee ete ee ee (3.13) 


Thus the method again leads to the Chew and Low (1956) effective range 
plot, and, although there is now only one free parameter, one can again 
fit the experimental data for the 33 phase with 

(9?/477) pn? = 0-08, 
or (f?/477) = 16, ene? © (0.13) 


in agreement with the photoproduction data. . 
As noted by Chew and Low, crossing symmetry requires that 


See je la ee el, 8.14) 
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In our present renormalized pseudo-scalar theory 


whereas for the cut-off pseudo-vector theory, (1), 
x=. 


The two theories thus can be made to coincide for the, experimentally 
well established, resonant 33 state, but give strikingly different predictions 
for the small p-wave phases. The experimental data with regard to these 
is still unclear, but when available it will discriminate between the two 
theories, within the limits, of course, of the reliability of this approach. 

It is worth remarking that, although the cut-off at nucleon mass of 
pseudo-vector theory appears to correspond in a rough way to renormali- 
zation of pseudo-scalar theory, this is only true of the final answer. In 
terms of individual graphs there is no apparent similarity between 
the two methods of calculation. 


$4. THE S-WavEsS 


The S-wave phase shifts may be calculated by precisely the same 
method. Taking only the leading term in the effective range expansion, 
the expressions corresponding to (3.5) are 


S®)— —i(2m)*f2/m, . Se Rice oa eee ee 
SQ ==90r*ftint 2" 0 een Un oe oe 


The corresponding expression for both the S-wave phase shifts is 


TT . 
p. cot =— na —JinArtam ga giles aeat4ss) 
This is very similar to the corresponding expression in the non-relativistic 
treatment of I. Using the value of g determined by the p-wave, it is of 
the observed order of magnitude, but again does not display the strong 
dependence on isotopic spin observed experimentally. 

In I, the possibility was considered of a strong t-dependent S-wave 
coming from the virtual effects of heavy mesons. The experimental 
result also appears to require that there should be a 7-independent term 
from this source, which cancels the ‘ leading ’ term calculated here. This 
appears to us to be the most probable explanation of this effect. 
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SUMMARY 


Calculations have been made of the absolute values of the cross sections 
for processes of the type 
zNA(n, Py)z_1N*, zN*(n, pojyzei N= 3 ZN4(n, n’),N4 and ZN“ (p, p’)zN4, 
which have been induced by nucleons with energies between 10 and 
20 Mev. 

It has been assumed that the target nucleus may be represented as a 
Fermi gas of nucleons, and that the nuclear processes proceed by direct 
interactions between the incident particle and the target nucleons, followed 
by the de-excitation of the nucleus ; the ejection of a particle can occur 
during either process. 

Comparisons between the calculated and experimental data for the 
total cross sections, energy spectra and angular distributions of the 
emitted particles, indicate satisfactory overall agreement except for the 
process ,N+(p, p’),N4. 

Empirical evidence is cited, which indicates differences in the level 
densities of the excited states of even—even, even—odd and odd—odd nuclei. 


§ 1. INTRODUCTION 


FEsHBACH et al. (1953) and subsequent workers (Culler e¢ al. 1956, Woods 
and Saxon 1954) have shown that the magnitudes of the total and 
differential cross sections for reactions between nucleons and complex 
nuclei may be successfully calculated by using a weak coupling theory. 
In these calculations a complex potential of the form V’,(1+-7¢) was used 
to describe the interaction between the incident nucleon and the target 
nucleus. 

Members of this laboratory (Morrison et al. 1955) and others (Lane and 
Wandel 1955, Clementel and Villi 1955) have shown that it is possible to 
predict the behaviour of ¢ as a function of the energy of the incident 
particle, and also its absolute magnitude. In demonstrating this result 
it was assumed that the nucleus could be represented as a Fermi gas of 
nucleons, and that the incident particle interacts individually with the 
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bound nucleons of the target with a cross section which differs only from 
that for free particles through the restrictions of the Pauli Exclusion 
Principle. 

In this paper we have extended this approach to the description of the 
nucleon—nuclear interaction to the calculation of the cross sections for the 
processes listed below ; in all of them previous workers have shown that 
the normal statistical theory of the decay of the compound nucleus is 
inadequate for explaining the experimental observations. 


(1) ,N4(n, p),_,N for neutrons with an energy of 14-5 Mev. 


(2) The excitation functions for the reactions : 
2Mg?4(n, py)y1Na*4, ,,A]?"(n, py)12Mg?’, 148i?*(n, py)13AP*, 
7CP7(n, py)igS*” and sgSr°*(n, Py) s7Rb**. 


(3) The angular distribution of the neutrons emitted with energies 
between 4 and 12 mey, following the interaction of neutrons of 14 Mev 
with bismuth. 


(4) The energy spectrum of the protons emitted at 60° and 150° from 
iron, silver and platinum when bombarded by protons of 18 Mev. 


§ 2. DIscUSSION OF PROCEDURE 


We have assumed that the target nucleus may be represented as a 
Fermi gas of nucleons. Whilst this choice was primarily dictated by the 
need for a mathematical approach which would be both simple and 
consistent, there appears to be some justification for considering it as an 
approximate physical description of the normal and excited states of a 
heavy nucleus. Thus an assembly of non-interacting particles should 
approximate reasonably to the requirements of the shell model of the 
nucleus, in which definite physical states must exist for considerable 
periods of nuclear time. In addition, Lang and Le Couteur (1954) have 
used the model of a Fermi gas to calculate the level densities of excited 
states of nuclei; their results are in good agreement with the experi- 
mental data. 

We have assumed that a nuclear reaction proceeds in three stages, 
(i) an initial collision between the incident nucleon and one of those 
contained in the target nucleus followed by (ii) the formation of an excited 
nucleus which subsequently decays (iii). The nuclear excitation, according 
to our model, is concentrated initially in one or more nucleons moving 
inside the nucleus at energies above the Fermi maximum but too low to 
escape from the nucleus. Degradation of this energy in further collisions 
can occur, and the nucleus can pass through several modes of excitation 
before stage (ili) occurs. The emission of nucleons by process (i) occurs in 
a period of about 10-8 sec, whilst by (iii) it occurs in times which are 
many orders of magnitude longer than this period. Thus process (iii) 
will be essentially independent of (ii), and we have assumed that the 
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statistical arguments of Weisskopf and Ewing (1940) may be applied to 
the nuclear decay, whilst the direct collision process is examined separately. 
Our procedure is discussed in detail in §§ 3 and 4. 


§3. THe Cross Sections ror Emission or NucLEONS BY Direct 
COLLISIONS 


3.1. Total Cross Sections for the Process ,N4(n, p),_,N4 


Culler et al. (1956) have shown that the experimental data on the 
magnitude of the cross section for the interaction of neutrons of 14 Mev 
with nuclei, may be satisfactorily reproduced by solving the Schrédinger 
equation for a complex potential well of the form V—42(1+-0-25%). 
As stated in § 1, imaginary potentials of this magnitude may be predicted 
by a Fermi gas model. Exact reproduction of the experimental data is 
difficult to achieve, since the theoretical results depend partly on the 
assumed shape of the potential at the surface of the nucleus. For 
convenience, therefore, we will use for the interaction cross sections o,, 
of the incident neutrons the measured values of Amaldi et al. (1946) and 
Phillips et al. (1952) at 14 Mev. 

We have assumed that this interaction occurs through a collision with 
either a neutron or proton in the target nucleus. Assuming that the 
neutrons and protons in the nucleus can be represented as two independent 
Fermi gases, we may write the cross section, (do/de), for the production 
of a proton with energy « in the target nucleus as 

° 7 py(dojde) ‘ 

(F N pa [aX (P 3) |npPpt [aX (P 3) |nnPn 

where X(P,)= the cross section for a collision of a neutron of momentum 
P, with a free nucleon, P,=1/(2M(V,+,)], M=mass of nucleon, 
V,=depth of nuclear potential well, H,—energy of the incident neutron 
outside the nucleus, «,=V,)+#,, do/d-e=cross section for the production 
of a proton with energy « in the collision between a single neutron and a 
proton, «=a factor which reduces the magnitude of X(P,) inside the 
nucleus through the operation of the Pauli Exclusion Principle and, 
p=the density of the nucleons inside the nucleus. 

The subscripts n and p refer to neutrons and protons respectively ; _ 
an expression for « will be given later in this section. 

We have defined the depth of the nuclear potential well Vj as the sum 
of the maximum Fermi energy «,; for the neutrons in the target nucleus 
and the binding energy of the last nucleon in that nucleus. The magni- 
tude of «, was calculated by assuming that the neutrons were contained 
within a sphere of radius 1-37 A1/310~18 cm (A is the atomic weight of the 
nucleus). We have assumed that V, for the protons equals that of the 
neutrons, and that the maximum Fermi energy of the protons ¢, is 
given by the difference between V and the binding energy of the last 
proton. The values for the binding energies were obtained from the 
tables of Feather (1953). 


2K 2 
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_ The expression for do/de in eqn. (1) may be written as 


do 3MX(P,) f Pee Py Pel dP, 


de 27Py° [Pree Py cit a9 | [Pit P2| 
25 3M X(P,) 2 2\1/2 =i 2P .+/[2(P +P") | 
~ /(2)P |, ame a P2432 Aon 
y 
—2/(2)P,tanh-1 5 | Re hale) 
Py 2% 
where y= (2M ep), c= / [2M (e+ ep—e,—Q)]. 


In the above equations P, represents the momentum of the struck proton 
(of energy «,). The cross section for scattering in the centre of momentum 
system of the two colliding particles has been assumed to be isotropic. 
Physically, this assumption is not correct, but the error introduced is 
negligible in this calculationt. In accord with experimental data, the 
cross section for the scattering of free neutrons and protons was assumed 
to vary inversely as the energy available in the centre of momentum 
system. In performing the integration the limits quoted above ensure 
that the Pauli Exclusion Principle is not violated. The symbol Q repre- 
sents the normal definition of the term defining the release of energy in a 
nuclear process. 

The term «X(P,) of eqn. (1) may be obtained by performing the 


integration 
y. &+Q da 
KPa) | (7) a A pate oe eee 


The final expression is unwieldy, but a numerically equivalent expression 
derived by an alternative method is given in eqn. (4). 


2 A 
SED RPS 

5) /Te 2 2 

(PALP 2) P P+ Pe 6PyA+ 6MQ) tanbes eV PS Pa! 
1 2 
x 
| +4V(2)P:8P)2—3MQ—P,) tanh — +2/(2)P.P?? 

ul z 
where z= 4/[2M (2e,—€,—Q)]. ot ee ened 


The energy of a proton which has escaped from the nucleus is given by 
H>=e— Vo. The cross section (do/dE#),,, for the direct ejection of a 
proton with kinetic energy between # and H+-dE when a neutron enters 


a nucleus is given by 
do fag do 
dE np aF (7) 9) ee erg Ce aco) 


last rough calculation indicated that the non-isotropy of neutron—proton 
scattering reduces do/de by about 20% in the energy region relevant to the 
present paper. 
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where ¢(H) represents the probability of escape from the nucleus. (see 
Appendix). 

The total cross section o,(dir.), for the direct emission of protons from 
the nucleus as a result of the first collision of a neutron with a proton can 
be obtained by the numerical integration of eqn. (5). 


o,(dir.)= | ; (i), a . . . ° 7 . (6): 


A rough numerical calculation showed that the probability for the 
direct emission of particles from second collisions is negligible, when 
the incident nucleon possesses an energy of less than 20 Mev. The 
total cross section for the emission of neutrons by direct collisions o,(dir.) 
may also be calculated by the methods outlined above. 


3.2. The Differential Cross Sections for Processes of the type 
m (psp )ZN* and ,N4(n,n'),N* 
For convenience, we will consider the process ,N“ (p, p’),N*. A caleu- 
lation of the differential cross section for the scattering of a nucleon by a 
Fermi gas of nucleons has been given by Hayakawa ef al. (1955), and so will 
not be repeated here. The result of their elegant treatment may be 
written as 
eo ATS eet V (p2*4t-2bp.?4-¢c) + p.*+b (7) 
de dw 47.qP}° \ 

where b=2"+-¢°—P?, c=(P?+q?—29")?, po =P pi — aq", pr = Py — 20 
+ P?—z,?2—P,7, %=(1/2¢) (P?+-¢?—P,”) and, q=P—P,. 

The differential cross section (d?a/dedw),°, for the production of a 
proton with energy « and in a direction @ inside the target nucleus may 
then be written as 


aa \2 oe [2( da ) = ( ds ) | 
G 7) [aX (P31) ]ppPptloX (P31) ]pnPn dedw] yy Pp dedw alls ; 


(9) 
The factor of 2 appears in the first term in the brackets since it represents 
the collision of identical particles. 


/(py*+ 26p,7-+-c)4-py?+6 


$4. Tae Cross SECTION FOR THE EMISSION OF A PARTICLE BY THE 
Decay oF A Compound NUCLEUS 
4.1. The Formation of a Compound Nucleus 
Following the initial collision between the incident nucleon and a 
bound one, an excited nucleus is left ; this nucleus may exist in various 
states of excitation. 


+ Equation (3.8) in the work of Hayakawa et al. (1955) is incorrect. Upon 
insertion of the correct expression the equations for b and c may be simplified 
to those given above. In addition we have allowed for the Q-value of a reaction. 
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When nucleon emission by direct collision occurs, the residual nucleus 
may have an energy of excitation which lies between 0 and F. Nuclei in 
this state of excitation cannot contribute significantly to the processes 
which we wish to discuss, and so will be ignored. 

When nucleon emission by direct collisions does not occur the nucleus 
is left with an excitation energy U, given by 


Un aR 


where B is the binding energy of the incident nucleon in the compound 
nucleus. We have defined a cross section o,, for the formation of a 
nucleus with excitation energy U by 


oto [ep dir.)-o, (dirk cee a en 


Thus o,/o; represents the probability that, subsequent to the initial 
collision, a compound nucleus will be formed with an excitation energy 
(£,-+-B) and, momentarily, no particles with energy greater than V>. 


4.2. The Decay of the Compound Nucleus 


In calculating the relative probabilities for the various modes of decay 
of the excited nucleus we have assumed that the statistical theory of Weiss- 
kopf and Ewing (1940) is essentially correct. Thus the cross section 
co, (comp.), for the emission of a particle a, by the decay of the excited 
nucleus may be written as 


o,(comp.) =o (10) 


La 
Cee 
v 
where J’,=partial width for the decay of the compound nucleus by 
emission of particle a, and J’;,=partial width for the decay of the com- 

pound nucleus by emission of a particle 7. 

A convenient formula for the term I” has been given by Le Couteury 
(1950). We have assumed that, at all energies, the particles resulting 
from the decay process are emitted isotropically. 

Recent comparisons (Paul and Clarke 1953, Forbes 1952, Cohen ef al. 
1954, Gugelot 1954, Graves and Rosen 1953) between the experimental 
data on the cross sections for the emission of charged particles and the 
predictions of the statistical theory have suggested that this theory is 
unsatisfactory. We believe that satisfactory comparisons have not been 
achieved through not taking into account the effect of direct collisions, 
and from a lack of suitable data on the level densities of the excited states 
of nuclei. The latter factor is of great importance, and will be con- 
sidered separately in § 5. 


+ For the term representing the inverse cross section which appears in this 
formula we have used the data of Shapiro (1953), with R=1-5 41/3 10-3 em, 
and of Beyster e¢ al. (1955) at 4 ev. This may be done without considerable 
error since cross sections for the formation of a compound nucleus are roughly 
the same for both weak and strong coupling theories. : 
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4.3. The Level Densities of Bacited Nuclei 


Recently, considerable experimental information has been obtained 
concerning the density w, of levels in many nuclei at excitation energies 
of about 8mey. This region of excitation energy is relevant to the 
present calculation, since that of the residual nucleus is in this energy 
range following the emission of a charged particle from the compound 
state. 

Lang and Le Couteur (1954) have made an extensive theoretical study 
of the spacing of nuclear levels D (where wD=1) and have shown that 
application of the thermodynamic properties of a Fermi gas to the 
nucleons of a nucleus yields the following equation for the spacing of 
nuclear levels of zero spin D, in a nucleus of atomic weight A and excita- 
tion energy U. 


AUNi2® 3 
Dy=0-1144(U-+0? exp.—| 2 (3) +35 (110) | aati ty) 


32 
where A= (“)" reins 
i * aes 
The spacing of nuclear levels of spin 7, D;, was assumed to be given by 
Tee rire coe tr ase (TO) 


In tables 1, 2 and 3 we compare the data on measured level spacings 
with the predictions of Lang and Le Couteur. The experimental data on 
the level spacings in light nuclei has been taken from the tables of Endt 
and Kluyver (1954) and from Buechner ef al. (1956), Paris et al. (1955) 
and Paul et al. (1956); that on heavy nuclei has been taken from the 
work of Harvey et a/. (1955). An average quantum number of spin was 
assigned tentatively to each group of levels, the assignment being made 
either from information obtained in the individual experiments, or from 
reasonable assumptions concerning the magnitude of the spin of the initial 
and product nuclei involved in the individual reactions and the energy of 
the incident particles. 

The agreement between the experimental and theoretical data in table 
2 is remarkable, in view of the variation of level spacing by a factor of 
104. It is apparent from tables 1, 2 and 3 that there is a large difference 
between the level density of even—even nuclei and that of other nuclei ; 
the existence of a difference between even—odd and odd—odd nuclei is 
not so certain. Taking an average of the final columns of the tables it 
is found that for nuclei of the same mass and excitation energy the 
average level densities are given by 

Woo Be Woe ae Woe 
12 5) 

These variations in the level densities have been used in the calculations 
on the decay of the compound nucleus. It is possible that, due to the lack 
of resolving power in the existing experimental techniques, some of the 
levels of excited nuclei have not been observed and, therefore, that the 
factors 12 and 5 could be higher. 


(13) 
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§ 5. RESULTS 


5.1. The Cross Section for the process, ,N4(n,p)z_,N, for Neutrons with 
an Energy of 14:5 Mev 

We have been primarily interested in calculating the magnitudes of the 
cross sections for the process ,N+(n,py),_,N“ and comparing these values 
with those obtained by Paul and Clarke (1953) and Forbes (1952), using 
activation techniques. Therefore, in making comparisons with these 
experiments we have assumed that the permitted limits of energy for the 
emission of protons are (#,-+@) and (#,+ q) where q¢ represents the binding 


Table 1. Odd—odd nuclei 


Measured 
U average | Average| (2/+-1)D, Expt. 
Nucleus Reaction in Mev level i} in ev Theory | 
spacing ; 
in ev 
35 
Aya (d,p) (2-545) | 15105 2 Tox 10° 0-4 
7:3 
gol (psy) (7-0-7-6) 4x 104 5/2 2-4 10° 1-0 
4-0 
(d,p) (2:25-5:85) | 7x 104 2 3-5 x 10° 0-3 
rol 
7-0 
(d,p) (6-3-7-7) 5x 104 3 ox 105 1-5 
3:2 
ise (d,p) (2-242) 210° 2 1-0 10° 0-7 | 
Neutron 
plies Scattering 7:2 14+2 9/2 (1-4+0-2)107 | 0-3 
pels i. 6-6 144-2 9/2 (1-4-+40-2)10? | 0-2 
pest4 ty 6:7 4245 7/2 (3-4+0-4)10? | 1-1 
Al er s 5:7 2:2+0:3 5/2 13-2+41-8 0-03 
Ped Monae : 5:8 10+1-0 3/2 40 +-4 0-13 
g7Ho16 : 5:7 1241-3 7/2 96 +10 0-4 
Pd Wana rf 5:9 15+2 1/2 30+4 0-2 
ye 93 6-0 7T+2 7/2 56 +16 0-6 
wa Lac? 6-0 9+1 7/2 72+8 0-9 
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Table 2, Even—odd, odd-even nuclei 


ce ae eee 


Measured 
U average | Average} (2/+-1) D,; | Expt. 
Nucleus Reaction | in Mev level I in ev ‘Theory 
spacing 
in ev 

5 

11 Na?4 (p.p’) (3-5-6-5) | 2-5 105 3/2 1-0 10° 1-0 
9-8 

11Na?8 (p-y) (9-3-10-3) 5x 104 1 LOnc le 1-4 
3°3 

12.Mg75 (d,p) (2-0-4-6) Saale 2 1-5 x 108 0-8 
; 4-1 

Al? (p.p’) (2-9-5-3) 3x 105 2 1-5 x 10° 1:3 
9-0 

il (p,y) (8-5-9-5) Tl0* 3/2 2-8 x 10° 2-8 
4-0 

Py ke (d,p) (2-0-6-0) 4 10° 2 2-0 x 108 1-9 
3°3 

nes (d,p) (2-3-4:3) 4x 10° 2 2-0 x 108 1-4 
9-3 

te (p.y) (9-0-9-6) | 2-5 104 2 T-25<¢10° 2-0 
3-1 

PS ed (d,x) (2-0-4-2) 2x 105 3/2 8 x 10° 0-6 
3-0 

apatt (d.p) —|(1-95-3-95)| 2-5 « 105 3/2 1-0 x 108 0:8 

Neutron 

. ag Mo%: 5,7, 9,101) Scattering 6-9 (5+1-5)10? 1/2 (10 +3)102 0-5 

=oSn * 8-0 (1-5-40-8)102| 1/2 (3-41-6)102| 1-2 

5oSn!!7 s 7:3 (1-5-40-7)102} 1/2 (3-1-4)102| 0-7 

eri! teks * 6-6 (2+1)10? 1/2 (4+2)10? 0-5 

go S121, 128, 125 xm 6-1 (5 +2)102 1/2 (10+4)102 | 0-9 

oo Hf179, 181 A 6-1 (1-40-5)102| 1/2 (2+1)102 | 28 

ss 4-9 18+2 1/2 36--4 0-7 
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Table 3. Even-even nuclei 
Measured 
U average |Average| (2/+1)Dy Expt. 
Nucleus | Reaction | in Mev level ie in ev Theory 
spacing 
in eV 
4-8 > 
(p,p’) (4-1-5:5) 4x 10° 3 2-8 x 108 3 
M. 24 
12Me “ae 
(p,y) (12-0-13-6)) 7104 2 3-5 x 10° 12 
4-3 
one (d,p) (3-0-5-6) ave 5/2 1-8 x 10° 2 
12-0 ; 
Fiske (p,y) (p,p’) |(11-8-12-25), 5x 104 5/2 3 LO 13 
4:8 
(p,p’) (3-8--5:8) S08 2 1-5 x 10° 3 
Si 
10-0 
(p,y) (9-2-10-9)| 1-1 104 3/2 4-5 x 105 21 
12-3 
1gA°® (pn) | (11-8-12-7)| 5108 5/2 3x 104 5-5 
Neutron 
a29Mo% | Scattering 9-2 (3-7 +1-2)10? 5/2 (22 +7)102 8:5 
loss i 8:3 (3-7 +1-2\102 5/2 (22 +'7)10? 4:5 
sist’ a 9-3 (1-2 +0-3)10? 1/2 | (2:-4+0-6)10? 5:3 
Epos. = 8-6 (3 -£1-4)102 1/2 (6+2-8)102 7:7 
Ble bie s 7-6 5-6 0-6 <3/2 2242-4 2:3 
pplehe tay uM 7-4 8+2 <3/2 3248 2:8 


The figures in brackets in the column for U indicate the limits of excitation 
energy over which the assessment of level spacing was made for light nuclei. 


energy of the last proton in the target nucleus; this corresponds to 
assuming that the appropriate radioactive decay does not take place if, 
after the ejection of a proton, the nucleus is left with sufficient excitation 
energy to emit a neutron. The cross sections were obtained by numerical 
integration of eqns. (6) and (10) between these limits. 

In addition we have calculated the total cross sections for the emission 
of protons, and compared these values with those obtained by Allan and 
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Lot Table 4. 
Element | Q Mev o(n,py) mb. a(n,py) +o(n,pn) mb. 
se Paul i 
op op (comp.) ene op op (comp.)| Experi- 
(comp.) | (dir.) + op Clarke lig (comp) | (dir.) +op ment. 
(dir.) (dir.) 
us —9-4 15 by) 20 49+ 25 25 15 40 
Ee 3:7 15 25 40 | 135445 45 75 120 
Na®* —3-4 LS 30 45 34+15 40 90 130 
Mg” —4:7 160 45 205 190+20 160 50 210 
Mg* 2.9 60 60 120 45-418 60 15 135 
ee —1:9 2D 40 65 52+10 | 80-5 75 100 175 140* 
sv —-3:8 ZS 45 260 220+50 230 50 280 
S12? —3-0 45 55 100 100+ 30 45 60 105 
pe —0-7 30 39 65 64+10 | 90410 105 105 210 
s® —()-9 170 50, 220 370+ 50 340 85 425 
S84 —4:3 70 30 100 85445 110 45 155 
Cl? —3-5 10 30 40 338+6 20 55 75 
ie —1:8 10 35 45 80-430 40 75 1s 
aes —3-l 100 30 130 93+ 30 Tihs 35 150 
hfe —1-4 10 30 40 27-+4 30 55 85 
Cr®? —3-0 aS 30 105 78-E 11 100 35 135 
Fe! 40-3 540 BD 595 | 580* 
Fe*t —2-9 60 30 90 97412 |124+12 80 30 110 190* 
site +0-6 470 60 530 530* 
Ni 29:0) is; 30 145 300* 
Net —0-5 30 45 715 182+30 30 45 75 
Cus3 S207 80 65 145 250* 
Cus a8 5 30 35 1944 15 45 60 | <80* 
Zn'4 +0-2] 110 45 | 155 | 386460 360 55 | 415 | 590t 
Zn* —1-9 45 30 75 100+17 95 Bo 130 
Ga®? —O0-1 5) 25 30 24118 30 55 85 
Ge? = (G9 60 30 90 130-+65 135 30 170 
Ge? —3-2 10 20 30 65-430 20 20 40 
Ge® —0:6 6 30 36 137+68 20 40 60 
As? —0-4 ey 30 33 Lee? 15 45 60 
Se77 0-0 5 20 25 45+23 25 25 50 
Sts —4:4 2 9 TL 18+3 33 9 12 
7,90 exalt 20 20 40 250+ 100 3 25 55 
wares * —4-6 0-3 6 6 10+5 1 9 10 
Mo?? ies 2 20 22 108+ 54 3 21 24 
Ruler E037 3 22 p45) 2+1 4, 21 25 
Pdi Sle 6 lye 23 132+ 66 10 20 30 
Pars | 10.2 4 25 29 | 7404520 6 24 30 10+ 
p27 0-0 1 16 17 | 230-4140 2 16 18 
Balss te 0 3 3 6+2 0 3 3 
TLals9 Te 0 9 9 642-5 0 9 9 
T]205 0 0 3 3 oe eb 0 3 3 
Pps | 4.2 0 0-3 0-3) 141 0 0-3 0-3 


* Allan. + Armstrong and Rosen 


The symbol (n,pn) indicates the emission of a proton followed by a neutron at a later 
period. The measurements of Allan and of Armstrong and Rosen cover the cross sections 


a(n,py), o(n,pn), o(n,pp) and o(n,np). 
by Armstrong and Rosen (private communications), using nuclear 
emulsions. The results for these calculations are shown in table 4. 

It is seen from the table that our calculations follow the trend of the 
experimental results. The large fluctuations in the calculated cross 
sections are produced by variations in both the Q-values and the level 
densities of excited nuclei. It should be noted that, in addition to the 


Cal 
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fluctuations resulting from the even-even, even—odd and odd-odd 
characteristics of nuclei, deviations from the formula of Lang and Le 
Couteur would be expected to arise from the detailed structure of an 
individual nucleus. These fluctuations will affect the magnitude of 
o,(comp.). 

It can be seen that calculated cross sections for the heaviest nuclei are 
due almost entirely to the emission of protons by direct collisions. In 
particular, the cross section predicted for lead by the decay of the com- 
pound nucleus is too small by a factor 104 when compared with the experi- 
mental value. 


5.2. The Excitation Functions for the Reactions ,.Mg?*(n, py),,Na**. 
Al’ (n, py)i2Mg*’, 1Si°"*(n, py) sAP®, 7CP%(n, py)ieS*” and 
ggdr**(11, Py)37Rb*S 


The excitation functions for ,,.Mg?4, ,,8i?°, ,,Cl®’ and ,,.Sr°8° have been 
measured by Cohen and White (1956), and that for ,,Al?’ has been 
measured by workers at Los Alamos (Hughes and Harvey 1955) by using 
activation techniques. We have calculated the cross section for the 
emission of a proton with energy lying between the limits (#,+Q) and 
(H,+¢q) by direct collisions and by the decay of a compound nucleus. 
Since the cross section for interaction of a neutron with the target nucleus 
does not vary appreciably with energy in the range 10-15 Mev (Taylor 
et al. 1955) we have used the values of co, measured at 14 Mev in calcula- 
ting the excitation functions of ,,.Mg?*, ,,Si?8, ,,Cl?? and ,,Sr88. The 
value of o; used in the calculation of the excitation function for ,,Al2* 
was 0:8 barns. 

The excitation functions are shown in figs. 1 and 2. 


Fig. 1 


- millibarns 
a 


(np) 


cs 3 4 5 6 7 8 


E Mev 
Excitation function for the reaction i3Al?"(n,py),.Mg?’. 
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The agreement between the calculated and experimental values of the 
cross sections is satisfactory except for Cl°’ ; we are unable to explain why 


the cross section for this nucleus should rise with increasing neutron 
energy. 


a(n p)- millibarns 


E Mev 


Excitation functions for the reactions ,,Mg”4(n, py),,Na**;  ,,Si?8(n, py)3Al’§ ; 
17C}?7(n, Py)igS?? and ggSr®8(n, py)37Rb®. 


o)(comp.)+-0,(dir.) reins 
5.3. The Angular Distribution of the Neutrons Emitted with Energies 


between 4 and 12 MeV, following the Interaction of Neutrons of 
14 MeV with Bismuth 


The experimental data with which we have made a comparison is that 
of Rosen and Stewart (1955). 
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The differential cross section for the production of a neutron, at different 
angles, with an energy « by direct collisions has been calculated from 
eqn. (8) with suitable alteration of subscripts. For simplicity of calcula- 
tion the parameters used were ¢,=«;=30 Mev and Vo= 36 Mev. It noe 
been assumed (Appendix) that the differential cross section (d?o/dEdw),, 
for emission of a neutron of energy / at an angle 6 to the direction of the 


incident neutron is given by 


hie NG d2a \ Ea (14) 
(ara) = (sey) 2°? - oti ar 


nn 
where \(<)—=the mean free path of a neutron of energy in nuclear matter, 
and R=the nuclear radius. 


Fig. 3 


Ore Teo 
Ome OlneO 


aN 
O 


o (6) -millibarns / steradian 


ce) 


20 40 60 80 100 120 140 160 180 
@-scattering angle 


Oo 


Angular distribution of the neutrons emitted from bismuth with energies 
between 4 and 12 mev. 


The angular distribution of the neutrons with energies between 4 and 
{2 Mev was obtained by the numerical integration of eqn. 14. The 
cross section for the emission of neutrons with an energy greater than 
4 Mev by the decay of a compound nucleus was found to be negligible. 

In fig. 3 a comparison is made between the calculated and experimental 


data. 


5.4. The Hnergy Spectrum of the Protons emitted at 60° and 150° from Iron, 
Silver and Platinum when Bombarded by 18 MeV Protons 

The magnitude of the term o, appearing in eqn. (8) is not satisfactorily 

known for protons with an energy of 18 Mev. We have employed values 
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Fig. 4 
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Energy spectra of the protons inelastically scattered at 60° and 150°, from iron, 


silver and platinum. 


Fig. 5 


Ratio 


E Mev 
Ratio of the cross sections-at 60° and 150° for protons inelastically scattered 
from iron, silver and platinum. 
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extrapolated from the tables of Shapiro (1953). For simplicity of calcula- 
tion the Fermi gas parameters used were ¢,=«,;=30 Mev and V)=38 Mev. 
Rough allowance was made for the effects of multiple reflection by the 
Coulomb barrier. 

The differential cross section for the emission of a proton of energy # 
by the decay of a compound ‘nucleus has also been calculated. The 
sun of this cross section and that for emission of protons produced by 
direct collisions is shown for iron, silver and platinum in fig. 4. For 
platinum o, (comp.) is virtually zero. 

It can be seen that the comparison becomes progressively worse with 
increasing atomic weight ; this might be due in part to uncertainties in 
the magnitude of the penetrabilities of the coulomb barrier, but it should 
be pointed out that, because of this barrier, surface interactions could be 
more important for reactions with protons than for those involving only 
neutrons. The ratio of the cross sections for emission of protons at 
60° and 150° is less sensitive to the penetrability; these ratios are dis- 
played in fig. 5. 


§ 6. Discussion 


It has been found that the nuclear reactions induced in heavy nuclei, 
by nucleons with energies greater than 100 Mev, may be successfully 
described by a mechanical model which assumes the production of a 
nucleon cascade followed by the decay of an excited residual nucleus 
(Bernadini et al. 1952, Goldberger 1948, Morrison et al. 1953). The work 
described in this paper employs a basically similar method, and yields 
surprisingly satisfactory results at an energy, for the primary nucleon, at 
which a more formal theoretical treatment might be expected to be 
necessary. 

It would be of interest to examine the region of validity of this model. 
In particular an inspection of the results quoted in table 4 would 
indicate that the angular distributions are strongly dependent on the 
nature of the target nucleus for the protons emitted in (n,p) reactions, 
since they will be determined by that relative strengths of the direct 
and compound processes ; whilst for nuclei with atomic weight greater 
than 100 the angular distributions of high energy neutrons from (n,n) 
reactions should always be of the form shown in fig. 3. 
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een NAD 1X. 


The Probability of Escape from a Nucleus, for Particles Produced in 
Direct Collisions 


Assuming that the initial collision can occur throughout the nuclear 
volume, it can be shown that the average distance a particle must travel 
to reach the surface is about 0-75R, in a nucleus of radius R. 

Thus the probability that a nucleon, of energy e«, may reach the surface 
of the nucleus and escape is given by the expression 


where A(<)=mean free path of a nucleon of energy « in the nucleus, and 
P=the probability that the particle crosses the nuclear boundary. We 
have assumed that P=1 for neutrons, =o,/c, for protons, where a, 
and o, represent the cross sections for the interaction of protons and 
neutrons respectively with nuclei. Values of o, were taken from the 
tables of Shapiro (1953). Since the values of o,, change very slowly with 
energy in the region relevant to the present calculations, we have used 
the data obtained by Beyster et al. (1955) for neutrons with an energy of 
4 Mev. 

Since P is less than unity and the mean free path, A(«), of the proton in 
nuclear matter is comparable with the dimensions of the nucleus, protons 
may emerge after more than one traversal of the nucleus following reflec- 
tion at the nuclear boundary. Thus we may write the total probability 
for the escape of the protons as 

GUID) IE Dil PEP ae aa a 
eee 
l-—z 
where x=(1—P) exp— [1-33R/A(e) ]. 

The term 1-33R represents the average distance the proton traverses 
in crossing the nucleus. The term x has negligible effect at high energies 
but can increase the emission of protons by about 40% at energies of 
about 5 Mev for nuclei of A~50. 
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The Electronic Structure of Palladium—Uranium Alloys 


By J. A. CaTTERALL 
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[Received November 14, 1956] 


SUMMARY 


The measured effects of uranium upon the paramagnetic susceptibility, 
electrical resistivity, thermoelectric power, and lattice parameter of 
palladium are discussed, and an interpretation of the results in terms of a 
solvent and solute band is suggested. 


§ 1. INTRODUCTION 
In a recent paper dealing with the constitution of uranium—palladium 
alloys (Catterall e¢ al. 1956) a considerable solid solubility of uranium in 
palladium is reported. Measurement of the effect of dissolved uranium 
upon certain of the physical properties of palladium suggested a means of 


Fig. 1 


N(E) 


Density of states of palladium. 


estimating the valency of uranium in the alloys, and the variation with 
composition of the lattice parameter, paramagnetic susceptibility, 
electrical resistivity, and thermoelectric power of the solid solution alloys 
are dealt with here. A comparison is also made, where possible, between 
the effects of uranium and silver on these physical properties, and an 
interpretation in terms of a solvent and solute band is suggested. 


Z2L2 
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In the solid state the electronic structure of palladium is considered to 
resemble that of nickel, but with 4d and 5s levels instead of 3d and 4s. 
In this model a long, low parabolic band contains the electrons of mainly 
5s-type, while a high, narrow inverted band, also parabolic, contains the 
electrons of mainly 4d-type. The relative position of these two bands is 
shown in fig. 1. From a consideration of the paramagnetic susceptibility 
of palladium, and that of its alloys with the noble metals, hydrogen and 
nickel, Wohlfarth (1948) concludes that there are 0-6 holes per atom 
remaining in the 4d band, and this figure is confirmed by Hoare et al. (1953) 
from the determination of the susceptibility and electronic specific heat 
of palladium-silver alloys. The latter workers also deduce a form for the 
shape of the head of the band. The relatively high susceptibility, 
electronic specific heat and electrical resistance of palladium are attributed 
to the unfilled d-band, while the s-band is considered to be responsible for 
conduction and cohesion. 

The band structure of uranium is not yet known, but in the free atom 
the outer electrons exist either in the states (7s)?(6d)1(5f)? (Seaborg 
1949) or in the states (7s)?(6d)* (Dawson 1952). 


§ 2. EXPERIMENTAL METHODS 


In this section the methods are described for the determination of the 
lattice parameter, electrical resistance and thermoelectric power. The 
paramagnetic susceptibilities discussed later are the values obtained by 
Leach (1955). 

The materials and methods used in the preparation of the alloys were 
the same as those of the constitutional investigation. The lattice 
parameters were determined at room temperature on quenched alloys by 
the powder method using a 9 cm Unicam camera and CoKz radiation. The 
filings were prepared from ingots homogenized at 970°c, and they were 
stress-relieved at the same temperature. The accuracy of the measure- 
ments, except for one point, is about +0:0004 A. 

The difficulty of preparing the alloys in the form of wire necessitated a 
somewhat crude method for the determination of the thermoelectric 
power, and no great accuracy is claimed for these measurements, the 
error being of the order of +-0-5 microvolts/degree. A copper tube was 
soldered at right angles to each end of an annealed specimen of approxi- 
mate dimensions 3 x 0-5 x 0-5 em, and steam was passed through one tube 
and water through the other. A platinum/platinum—rhodium thermo- 
couple was attached to each end of the specimen in order to measure the 
temperatures, and the platinum wires were then used to measure the 
e.m.f. of the alloy against platinum. The whole assembly was enclosed 
in a heavily lagged wooden box. The equipment was calibrated against 
nickel, which also served to check the signs of the e.m.fs. 

The electrical resistances of annealed specimens in the form of strip 
approximately 8 x 0-5 0-1 cm were measured at 20°c in the Electricity 


Division, National Physical Laboratory, and they are considered to be 
accurate to within +0-5%,. 
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The chemical analysis of all the specimens was performed in the 
Metallurgy Division, N.P.L. All compositions are reported in atomic per 
cent. 

§ 3. EXPERIMENTAL RESULTS 


3.1. The Lattice—Parameter Composition Curve 


The variation of lattice-parameter with composition is shown in fig. 2. 
Solution of the larger uranium atom expands the palladium lattice, whilst 


Fig. 2 


Lattice 
Parameter 


3 :°9500 


3°9000 


5 10 Is 20 


Atomic per cent Uranium 


The lattice parameter—composition curve. 


a marked change of slope in the direction of an additional expansion occurs 
at a composition between 10 and 10-5%, uranium. The curve has been 
determined up to the solubility limit at 970°C of 22-3°% uranium. 


3.2. The Paramagnetic—Susceptibility Composition Curve 
Leach’s values at room temperature are shown in fig. 3, together with 
those of Hoare et al. for palladium-silver. In this figure, and also in 
figs. 4 and 5, the axis representing the percentage of uranium has been 
drawn six times the scale of the axis representing the percentage of silver. 
In this way if uranium were to dissolve with a valency of six and the 
changes in physical properties were attributable primarily to the influence 
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of the valency electrons of the added element the curves corresponding 
to the uranium and silver alloys of figs. 3, 4 and 5 should be superimposed. 
It is evident that at compositions up to 5% of uranium, uranium is six 
times as effective as silver in reducing the paramagnetism of palladium. 
At higher concentrations it becomes less effective, and the paramagnetism 
rises in alloys containing more than 9°% uranium. 
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Paramagnetic susceptibility—composition curve at 20°c. 


3.3. The Hlectrical—Resistivity Composition Curve 


The variation of resistivity with composition is shown in fig. 4. 
A maximum occurs at 10%, uranium and a minimum at 11-39%. Beyond 
this point the curve rises extremely steeply. The resistivity of palladium— 
silver alloys at 20°c determined by Kemp e¢ al. (1956) is also shown in 
fig. 4. The shape of the uranium curve in the region 8-12°% of solute 
depends upon only three alloys, and in order to avoid the possibility of a 
mistake the compositions of these three alloys (which contained 8-5 
10-5 and 11-3°% uranium) were checked ; and their resistances rodeter. 
mined. In addition although the constitutional investigation has shown 
the absence of order—disorder transformations at temperatures down to 
700°C, the alloy containing 10-5°% uranium was re-annealed for 24 hours 
at 500°c. The resistance was unaltered after this treatment, and it ma 
be concluded that such transformations have no influence on the sha Z 
of the resistivity curve. : 
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3.4. The Thermoelectric-Power Composition Curve 
The variation of the thermoelectric power against platinum with 
composition is shown in fig. 5. The values for palladium-silver alloys 
obtained by Giebel (1911), which are qualitatively confirmed by Taylor 
and Coles (1956), are also shown in this figure. 
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§ 4. Discussion 

In calculating the heats of formation of binary alloys Varley (1954) 
has proposed a two-band model for the distribution of the energies of the 
electrons in a disordered solid solution. In this model the energies of 
the electrons in cells surrounding the solute and solvent atoms resemble 
those of the pure solute and solvent metals respectively. The energy of 
the system is then lowered by electron transfer between the cells if the 
electrons in one type of cell have higher energies than electrons in the 
other type. Varley has used this model to interpret the magnetic 
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properties of copper—nickel alloys, in which he assumes that electrons 
are transferred from the copper cells to the nickel cells. 

A similar model may be used to account for the behaviour of palladium- 
uranium alloys. The distribution of the energy levels in palladium metal 
is shown in fig. 1. The distribution of the levels in uranium metal is not 
yet known, but the x-ray absorption spectra obtained by Cauchois (1952) 
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indicate a high density of d- and f-states at the Fermi surface. It is not 
unreasonable therefore to assume that the energy band in uranium 
metal is derived from the 7s, 6d and 5f levels of the free atom, with the 
d- and f-levels forming a high narrow band overlapping a low band of 
mainly 7s charactery. 

From an examination of the susceptibility curves it is evident that at 
low concentrations uranium dissolves in palladium with a valency of six, 
since the paramagnetic susceptibility is reduced at six times the rate for 
silver. This shows that initial additions of uranium to palladium cause a 
large electron transfer from the uranium cells to the palladium cells, the 
electrons occupying the empty 4d-levels of palladium and consequently 
decreasing the paramagnetism. A continuation of this process with alloys 
of progressively increasing uranium content, however, will not explain 
the subsequent rise in paramagnetism. It is suggested, therefore, that the 
relative position of the solvent and solute bands shown in fig. 6 accounts 
for the observed effects. The occupied states shown as shaded areas in 
this figure correspond to the conditions existing in very dilute alloys. 
Further additions of uranium cause further electron transfer and also 
raise the Fermi surface. The system will eventually arrive at a stage 
where two opposing effects are taking place. These effects are (a) the 
occupation of the remaining empty d-levels in the palladium cells, which 
will cause the paramagnetism to fall, and (6) the occupation of the empty 
d- and f-levels in the uranium cells which will cause the paramagnetism 
to rise. A minimum in the curve of susceptibility against composition 
will therefore be expected. 

The two-band model may also be used to interpret the resistivity against 
composition curve shown in fig. 4. The high electrical resistance of 
palladium has been shown (Mott and Jones 1936) to arise from the 
scattering of the conduction electrons into empty d-levels, which occurs 
in addition to the contribution from s — s scattering. The probability 
of scattering is proportional to the density of states at the Fermi surface. 
Initial additions of uranium to palladium cause the resistance of the alloys 
to rise as the foreign atoms introduce perturbations. This effect is 
opposed with increasing uranium content by a tendency for the resistance 
to fall as the empty d-levels of the palladium are filled, and a maximum 
similar to that found in palladium-silver alloys appears. The resistance 
then rises again as the alloys become still richer in the solute because of 
the presence of the empty d- and f-levels of the uranium in the neighbour- 
hood of the Fermi surface. 

The change of slope of the curve of lattice-parameter against 
composition which occurs at 10 atomic’, uranium suggests a weakening 
of the cohesive forces existing in the alloy at the point where the palladium 
d-band is completely occupied. Until fairly recently the cohesion in the . 
end members of the First Long Period was regarded as arising from the 
conduction electrons, but Friedel (1952) deduced a large Van der Waals 
i gl 
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force between the 3d-shells in copper, amounting to one-third of the 
binding energy. Mott (1953) has suggested that in nickel the Van der 
Waals contribution may exceed that from the conduction electrons and 
that its presence may account for the low number (0-6) of conduction 
electrons in both nickel and palladium. Friedel (1954) does not now 
consider the description of the d-shell interaction as a Van der Waals 
force to be reliable. However the change of slope in the palladium— 
uranium alloys appears to be connected with a variation in the rate of 
decrease of some form of d-shell interaction between the palladium atoms. 

If the relative band positions shown in fig. 6 are correct they provide an 
explanation for the extensive solubility of uranium in palladium. The 
Fermi energy must rise relatively slowly per substituted solute atom in 
order that the structure should remain stable, and for this to occur a high 
N(E) curve is necessary. For dilute alloys this condition is provided by 
the 4d-levels in the palladium cells, and for richer alloys by the d- and 
f-levels in the uranium cells. 
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The y-Rays from the 8:06 MeV Level in “N} 


By C. Broups, L. L. Green, J. J. Stncu and J. GC. Witimorr 
Nuclear Physics Research Laboratory, University of Liverpool 


[Received November 9, 1956] 


§ 1. INTRODUCTION 


THE y-rays from the 8-06 Mey level in !4N which is formed when 550 kev 
protons are captured by °C have been examined by several workers and 
are of considerable interest. The resonance is an s-wave resonance of 
width 32-5 kev which is 12°, of the Wigner s-wave limit, indicating that 
the level may be a single particle level (Clegg and Wilkinson 1953). 
The main E1 transitions from this level have been discussed by Lane and 
Radicatti (1954) and by Wilkinson (1956). The level is known to be 
1~, 7=1 and the E1 transition to the 2-31 Mev 0+ 7'=1 level is isotopic 
spin forbidden. Attempts to observe this transition by Clegg and 
Wilkinson (1953) resulted in an upper limit of 0-7°% of the 8-06 Mev 
transition. Wilkinson (1956) has recently pointed out that the 8-06 Mev 
level may be one for which the isotopic spin impurity is particularly high. 
It is the purpose of this paper to report a measurement of the intensity 
of this unfavoured E1 transition and also other weak transitions. Similar 
measurements have recently been carried out at Saclay by Lehman et al. 
(1956) who obtained very similar experimental results. 


§ 2. EXPERIMENTAL RESULTS 


A thin target of C was bombarded with 550 kev protons and the 
y-rays observed at 0° and 90° to the proton beam. The y-rays were 
observed in a 5in. diameter by 6in. thick Nal crystal mounted on a 
5 in. diameter E.M.I. photomultiplier and the pulses recorded by a pulse 
height analyser. The irradiation of the crystal was restricted to the central 
2in. of the face by a lead collimator 4in. thick. The methods for 


y-ray energy in Mev 8:06 | 5-73 | 4-91 | 4- 
3: 


1 | 3-38 | 3-13 | 2-30 | 1-61 
Relative Intensity 100 36.) £8: {1 2:5 | 2 


1 1 
1 “() | 19:0 | 12-8 


obtaining y-ray intensities from complex spectra have been described in 
earlier papers. The spectra obtained at this resonance are shown in 
figs. 1 and 2. The y-ray energies and relative intensities observed at this 
resonance are shown in the table. The relevant levels in '4N and the 
transitions to which these y-rays are ascribed are shown in fig. 3. The 
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intensities are in excellent agreement with this scheme except for an 
excess of the 2-30 Mev y-ray which is due to a contribution from the 
RC(p, y) reaction. 

The region above 4 Mev is shown in detail in fig. 2. The peaks A, B 
and C are the three peaks due to the 8-06 Mev ground state transition. The 
peak D is of energy 5-73 mev and is of the correct energy for the isotopic 
spin forbidden transition to the 2°31 Mey level. Unfortunately there is also 
a level at 5-69 Mev which could also give rise to a y-ray of very nearly this 
energy. The decay of this state has been reported by Bent et al. (1955), and 


Fig. 3 
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by Thomas and Lauritsen (1953). Both groups report that the 5-69 Mev 
level branches to the ground state and the 2-31 Mev state in the ratio of 
30 to 70, the latter transition giving a y-ray of 3-38 Mev. Peak number 1 
of fig. 1 is the total capture peak of a y-ray of this energy and it has an 
intensity of 2:5°% of the ground state transition. The corresponding 
intensity of the 5-69 Mev y-ray, using the above branching ratio, should 
thus be 1-1°% of the ground state transition, whereas peak D of fig. 2 
corresponds to a y-ray of 3:6% of the intensity of the ground state 
transition, Hence if the previously reported branching ratio is correct, 
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two-thirds of the 5:73 Mev peak is due to the isotopic spin forbidden 
transition to the 2-31 Mev level. The broad width of this resonance 
enables us to check this interpretation. The 5-75 Mev y-ray is a primary 
y-ray whose energy should change with bombarding energy, while the 
5-69 Mev y-ray should be of constant energy. The region of the spectrum 
containing peak D of fig. 2 was observed under large energy expansion 
with the proton bombarding energy at 30 kv below the 550 kev resonant 
energy and again at 30 kv above. In the first case the y-rays should be 
of energies 5-69 and 5-72 Mev and in the latter case 5-69 Mev and 5-78 Mev. 
The two spectra are shown in figs. 4 and 5. The shift of the mean energy 
of the single escape peak leaves the 4-91 Mev y-ray more clearly resolved. 
The pulse height curves of y-rays of 5-69 Mev energy have been carefully 
determined. We have constructed the expected shapes for a composite 
peak of y-rays of energy 5-69 and 5-72 Mev and in the intensity ratios 
to be expected from the measured intensity of the 3-39 Mev y-ray and the 
branching ratio of Bent ef al. (1955). This enables us to determine from 
fig. 4 the exact shape of the tail from the 8-06 Mev y-ray. Ifa tail of the 
same shape is subtracted from the pulse height distribution shown in 
fig. 5 the peak remaining is both broader and at a slightly higher energy, 
thus confirming that the 5-73 Mev peak does contain a high proportion 
of a primary y-ray. 


§ 3. DiscussIon 


The intensity of the isotopic spin unfavoured transition to the 2-31 Mev 
state is approximately four times the limit set by Clegg and Wilkinson 
on this transition and corresponds to a contamination of the 8-06 Mev 
T=1 state with 8% of a 7’=0 state in intensity ; in the notation of 
Radicatti (1954) «,? (0)~0-08. This assumes that the 2-31 Mev state is 
a pure isotopic spin state; MacDonald (1956) has calculated that the 
impurity of the 4N ground state should not be greater than 3x 10-3 
and the 2-31 Mev state should not be very different from this. The 
8°% impurity is rather a high contamination for a state of this excitation 
and a possible reason for this has already been pointed out by Wilkinson. 
The width of the 8-06 Mev level indicates that this state probably has a 
unique parent in the ground state of %C. The interaction between the 
lp and 2s orbits is weak, as has been pointed out by Lane (1954). Thus 
thomicummavates: 0-0 = Sl Vests 07, V=0F 17) 70, of which 
the 8-06 Mev level is the 1~ 7-1 state, should not be far separated in energy. 
The 1~- 7'=0 state should be very favourably situated for contaminating 
the 8-06 Mev 1- 7’=1 state. The large impurity measured here lends 
some further support for Wilkinson’s suggestion. 

The 3-13 Mev y-ray corresponds to a transition from the resonance to 
the 4:91 mev state. This is known to be 0~ or 1~, of which the 0~ is to be 
strongly favoured for reasons given by Ajzenberg and Lauritsen (1955). 
The matrix element | M [? for this transition is 0-3 Weisskopf M.1. units, 
indicating that it is probably M.1. in agreement with the 0~ assignment. 
As no change of configuration is allowed for an M,1. transition, this means 
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that the 4-91 mev level is probably the 0~, 77=0 member of the quartet 
of states based on the ground state of #%C with an odd 2s particle. On 
this assumption we have calculated the transition probability for the 
3:13 Mev y-ray, following Lane and Radicatti (1954), in jj-coupling, 
including the contributions from the 2s particle and the p,/, particle. In 
terms of their parameter A the results are as follows : 

Aexpt. Ajj 

2:36 2-66 

This M.1. transition probability is thus consistent with assumption that 
the 4:91 mev levelis the 0-, 7/0 member of the quartet mentioned above. 
It may be noted that the 8-70 Mev state in 14N is also formed by s-wave 
proton capture in !8C and has a large width. This state is known to be 
0-, 7=1 and is probably a third member of this quartet and indicates a 
spacing between similar levels of different isotopic spin of this quartet of 
about 4 Mev which is of the order of magnitude required to produce the 
measured isotopic spin impurity. 

The other level fed from the 8-06 Mev level is the 5-69 Mev which is 
fed with about the same intensity as the 4-91 Mev state. The spin and 
parity of this state are not known but the intensity of this y-ray is correct 
for an M.1. transition which would be the case if this level were the fourth 
member of the quartet and 1~, 7=0. It may be noted that the decay of 
the 5:69 Mev level is to the ground and first excited state and the intensities 
of these transitions are just those to be expected for a 1~, 7’=0 state to the 
ground and first excited state with the 7'=0, 1~ state contaminated with 
8% of 7’ =1 in intensity. 
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The Effect of Added Titanium and Aluminium on the Magnetic 
Behaviour of a Ferric Oxide} 


By G. Haier 
Department of Physics, Imperial College of Science and Technology, London 


[Received November 24, 1956] 


ABSTRACT 


An examination is made of the magnetic behaviour of hematite, 
« Fe,O;, in the region of the antiferromagnetic transition temperature 
(—15°c) when titanium and aluminium are added in small amounts as 
an impurity. Contrary to an earlier observation by Morin, the transition 
is still observed in the synthesized material for relatively large impurity 
contents (~10%). The magnetic properties and behaviour of some 
naturally occurring hematites are compared with the synthesized 
materials, and are shown to be very structure sensitive. Thus the fact 
that this antiferromagnetic transition is not observed in natural hematite- 
bearing sandstones is thus probably not due to the presence of impurities 
but to some effect of crystal structure. 


§ 1.1 


In 1952-53 a series of measurements were made by Clegg et al. (1954) 
of the remanent magnetization of Triassic sediments in England. Red 
rocks, characteristic of this geological period, were collected from widely 
separated sites and were found to have consistent polarizations, being 
magnetized in an approximately north-east to south-west direction 
with dips significantly less than that of the present earth’s field, and with 
about half the specimens showing reverse polarization. The important 
palaeomagnetic implications of these results, both to continental move- 
ment and to the history of the earth’s magnetic field, made it desirable 
to investigate the origin of the natural magnetism of these rocks. 

An examination of the thermal properties of these rocks by Leng 
(1955) has revealed that the magnetic constituent is hematite, « Fe,O3, 
but that it exists in two different forms. One is a detrital form consisting 
of black specks of material which x-ray powder photography revealed 
to be hematite, but which gave a Curie point in the region of 550 to 
580°c, significantly lower than the Curie temperature of pure hematite, 
which is 685°c. The second is a red form appearing as a coating or surface 
deposit on the silicaceous grains, or as a cement between these grains. 
The latter was found to be much less magnetic than the former, but with 
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a much higher coercive force. It has been shown by Leng that it is this. 
latter form which gives the sandstones their natural remanent 
magnetization. The Curie point of this second form is observed to be 
closely that of pure hematite. 

Leng found no detectable difference between the magnetic properties 
of the normal and reversely magnetized specimens, thus strengthening 
the evidence for, but by no means finally proving the reality of repeated 
reversals of the earth’s magnetic field. A general discussion of the 
argument. has recently been given by Blackett (1956). 


§ 1.2 


For a further study of the possible origin of magnetic reversals in such 
rocks, a further knowledge of the magnetic properties of hematite is 
required, not only at normal temperatures, but throughout the possible 
range of terrestrial temperature. 

Morin (1950) has observed the effects of cooling on the susceptibility 
of powder specimens of « Fe,O, and, in particular, of specimens of 
« Fe,O, with small quantities of titanium added. Morin found that the 
susceptibility of « Fe,O,, and « Fe,0, with small percentages of titanium 
added (up to 0-2% Ti) exhibited a marked transition at —15°c when 
cooled : the susceptibility decreasing rapidly at this temperature, to be 
about one-ninth of the room temperature susceptibility. This transition, 
however, was not observed when titanium was present to the extent of 
{oi 

Similar observations were made by Pauthenet (1952, 1953) on a 
single crystal of natural rhombohedral « Fe,0, from Elba; he also 
observed the transition occurring in the remanent magnetization, that is, 
when measured in zero field. 

This phenomenon has been described theoretically by Neel (1952), 
wherein he attributes the transition to one of two parasitic ferromagnetisms 
which are present in hematite together with the fundamental 
antiferromagnetism. One of the parasitic ferromagnetisms he assumes 
to be closely tied to the direction of antiferromagnetism, which at room 
temperature and above is in the basal plane at right angles to the trigonal 
axis. On cooling below —15°c, the direction of antiferromagnetism 
turns into the trigonal axis taking a component of ferromagnetism with 
it. This change produces the decrease in remanent magnetization 
observed at this temperature. 


§ 2.1 


In view of these results, it seemed likely that by observing the remanent 
magnetization of a natural or saturated British Triassic sandstone 
(a typical specimen contains some 5% of hematite (Blackett, p. 55)), 
as a function of temperature, a similar transition at or about —15°c 
should be found. Figures 1 (a) and 1 (b) show the remanent magnetiza- 
tion of rock specimens as measured by an astatic magnetometer, which 
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(6) Saturation magnetization of Triassic sandstone while being cooled. 
(Saturating field 6000 oersted.) 
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were cooled to —80°c. It will be seen that no transition is observed, 
either from the natural or the saturation magnetization. 

Spectroscopic analysis has since shown that for detrital particles at 
least, about 5°, of titanium is usually to be found in the hematite in the 
‘Triassic rocks. This, according to Morin’s results, would account for the 
transition not being observed in these rocks. 
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ae ‘ 
Variation of the saturation remanence of hematite with temperature 
(measured in zero magnetic field). 


§ 2.2 
By way of comparison, fig. 2 shows the behaviour of the remanent 
magnetization of a specimen of chemically pure hematite when cooled 
to —80°c. It is interesting to note that the reheating curve does not 
follow the cooling curve, but by replacing the specimen in the same 
magnetizing field as initially, the portion F to A is completed and the 
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cycle can then be retraced. This cyclic hysteresis effect, involving both 
temperature change and remagnetization, will be further described 
in a later paper. 


SS 

The present paper describes experiments to elucidate more fully the 
manner in which titanium, and also aluminium impurities, influenced 
the Morin transition at —15°c. : 

I am indebted to Mr. D. E. Cox of the Chemistry Department of 
Imperial College, for making up a series of samples in which varying 
amounts of titanium were added to pure « Fe,0O,. The method of 
preparation was as follows: powders of high purity (at least 99-9°/), 
a Fe,O; and TiO,, were weighed, mixed and ground in an agate mortar. 
Compositions in which titanium (expressed as a percentage of the total 
metal atoms in the mixture) was 0-1, 0-3, 0:5, 0-7, 0-85 and 1-0 were made 
up. The mixtures were pressed into pellets under a pressure of about 
18 000 p.s.i., and fired for about 16 hours at 1100°c in an atmosphere of 
oxygen. The fired pellets (grey in colour) were ground to a fine powder 
in an agate mortar. Small discs, ~21 mm diameter and 4 mm thick, 
were then made, using 0-5 gm of the magnetic material and 2-0 gm of 
plaster of Paris. These discs, magnetized diametrally in a field of just 
over 6000 oersted, were then examined individually in a cooling chamber 
situated in zero field beneath an astatic magnetometer. The variations 
of the remanent moment of the specimen could then be followed during 
the cooling, measurements of the magnetization being made at 10° 
intervals from +15°c down to —95°c. Because of the hysteresis effect 
observed during such a temperature cycle, measurements of the 
magnetization were always made during the cooling. 

A measure of the temperature at which this transition occurs can be 
more easily obtained by taking from the cooling curve, the temperature 
T,, at which the magnetization has fallen to a value given by 
Trrp=2Irinit.tr,r,)> Where Jp, is the magnetization at the 
temperature 7',: this being the temperature corresponding to the 
magnetic state in the low temperature region D of fig. 2, the decrease 
in magnetization being assumed to be complete at 7';. It is the variation 
of the ‘ characteristic temperature ’ 7',, with titanium content, which is 
examined here. 

A selection of cooling-curves for specimens of this series is shown in 
fig. 3: the heating portions of the curves being omitted for clarity. 
These curves may be compared with curves obtained by cooling a natural, 
and a magnetically saturated, Wirral sandstone (figs. 1 (a) and 1 (b)). 

It was found that as the titanium content was increased from 0% to 
0-98%, the characteristic temperature varied in a continuous manner 
from —14-0°c to —51-0°c. Even for a titanium content of 192% the 
transition was still observed. A curve showing the variation of 7’, 
with percentage Ti content is given in fig. 4. From this curve it appears 
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that 7’, reaches a lower limit of about —52-0°c as the percentage of 
titanium is increased. This non-disappearance of the transition with 
increased titanium content conflicts with the work of Morin. 

An attempt was made to increase still further the titanium content in 
% Fe,O; by the method described, but on examination of the resulting 
powder by means of x-ray powder photography, lines other than those of 
% FeO, were found. As these only occurred when more than 5°, 
titanium was present, it must be assumed that for these samples the 
“Fe,O; lattice had not accommodated all the Ti atoms present. 

Although other compounds were found in the samples containing 
5% and 10%, Ti, a specimen was made up from the 10% sample and 
examined in the low temperature region. Here again the transition was 
still observed, but with a 7',-value of —53°c. 


§ 2.4 

From these observations it would seem that titanium is not the cause 
of the disappearance or suppression of the transition in « Fe,O, at low 
temperatures. Consequently, it would appear that the suppression of 
the transition in the natural rock must be due to some factor other than 
the titanium impurity. It was, therefore, useful to investigate if the 
second most common impurity associated with the magnetic fraction 
of the natural rock, aluminium, was the cause of the suppression, or if 
it be due to some other entirely different mechanism which has not so far 
been considered. 


§3.1 


The first of these possibilities has been examined in some detail. 
A further series of samples was made by Mr. D. E. Cox, in which varying 
amounts of aluminium were introduced as an impurity into « Fe,Og. 
The amount of Al present being expressed as a percentage of the total 
metal ions in the sample and the method of preparation was identical 
with that for the Ti series: the Al ions being introduced in the form of 
Al,O;. Samples were made with the following percentages of Al: 
0-099, 1-005, 1-99, 4-986, 9-996. As with the Ti series, the samples were 
examined by means of x-ray powder photographs for other phases present. 
The general result indicated a higher solubility of Al in « Fe,O, than 
with the Ti: only a few weak lines appeared in the 5% and 10°, samples, 
indicating the presence of small amounts of a second phase. 

As with the Ti series, a value of 7’, was determined for each sample ; 
and fig. 4, the curve of 7’, plotted against the percentage of Al content 
in « Fe,O3, shows that the decrease in 7’, with increase in Al is linear 
within the limits of Al concentration considered here: the transition 
again being found with high percentages of impurity. A comparison 
of the curves of fig. 4 would also seem to indicate a greater solubility of 
Al,O, in « Fe,O, than of TiO, in « Fe,O;: in the latter one gets a 
“ saturation ’ effect for impurity percentages greater than 1%. Comparing 
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7,-values of the two series, those for the aluminium series are much 
Hieher than those for the titanium series for Coens impurity 
concentrations : 0-17°/, Ti being as effective as 10% Al in lowering the 
T p-value. 

Teromn these two experiments it can be concluded that foreign ions of 
Ti and Al, when introduced into « Fe,O, by the method described, do 
not cause a suppression of the tnecee raat below room temperature. 


§ 3.2 

There is the possibility, particularly in the case of the titanium impurity, 
that Ti ions introduced in amounts in excess of 1° are not taken up by _ 
the « Fe,O, lattice, but remain in the form of a titanium oxide as a 
simple mixture with the « Fe,0;. If this were the case, the x-ray 
picture of the 2% Ti in « Fe,O, sample would be expected to show lines 
other than those of « Fe,O;; this is not found to be so. It is only with 
the 4:8°%, Ti in « Fe,O;, sample that extra lines begin to appear, though 
very faintly : these lines being the stronger lines of the compounds ferric 
titanate (Fe,O0,.TiO,), and cubic iron-titanium oxide (FeO . Ti,O,). 
From this it may be concluded that the amount of Ti wholly accommodated 
by the « Fe,O, lattice, by the method of preparation used here, is certainly 
greater than 1% and probably very nearly 4:8°%. 

These conclusions on the miscibility of the foreign ions Ti and Al in 
# Fe,O, seem to be in fairly good agreement with the general results. 
found by other workers (see review by Nicholls (1955)). In general it 
has been found that Fe**+* ions in hematite may be replaced by larger 
amounts of Al ions than by Ti ions; replacement by Ti ions seems to be 
~10%, whilst the corresponding figure for Al ions is variously quoted 
as ~10% (Nicholls 1955) to ~36°% (Yoshimura 1935). By the method 
of preparation of the samples described here, nearly 5°% of either Ti 
or Al ions could be introduced into the « Fe,0; before compounds began. 
to be formed. In no case, however, was the transition suppressed, 
even though higher percentages were used here than were used by Morin.. 


§ 3.3 

As a consequence of these results it is evident that two outstanding 
problems arise. In the first place, the results do not agree with Morin’s. 
observations on the suppression of the transition, and secondly, they 
fail to explain the suppression of the transition in ae case of the Triassic: 
sandstones. 

It is significant to note here that whilst the introduction of foreign. 
ions to hematite produces a lowering of the transition temperature, 
on no occasion, as a result of any treatment of the hematite, has the- 
transition been observed to occur above —13°c. 


§ 4.1 


Whilst we have not observed the transition in any natural sandstones. 
we have observed it in a sample of massive hematite, just as Pauthenet 
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observed it in natural rhombohedral hematite. Common massive 
hematite, known geologically by the name ‘ red-hematite ’, is brownish 
red in colour, generally possessing a fibrous structure which is cut 
transversely by a curvilinear banding, and with a surface which is often 
botryoidal or reniform. The botryoidal masses with radial-fibrous 
structure, common for instance in Cumberland, break up into conical 
forms known there as pencil ore. 

On the other hand, as will be shown later, specimens of specular 
hematite studied do not show an appreciable transition effect. These 
crystals of specular hematite are brilliant steely plates with modified 
edges, and are usually opaque, but very thin plates when viewed in 
transmitted light appear a blood-red colour. 


§ 4.2 

The direction of magnetization, with respect to the fibre-structure 
in the case of massive hematite, has been found to have an appreciable 
effect on the shape of the curves in the temperature cycle. Figure 5 (a) 
shows the temperature hysteresis of a specimen magnetized in about 
1000 oersted in a direction parallel to the fibre-lines, and fig. 5 (b) shows 
the complementary curves for a similar specimen magnetized in a direction 
perpendicular to the fibre-lines. The most significant feature common 
to both these curves is the relatively low temperature at which the 
transition effects are observed. Following the procedure in § 2.3, the 
effective transition temperatures in these curves are of the order of 
—60°c or —70°c, which is somewhat lower than any previously observed 
value. The low temperature at which the transition takes place is 
further emphasized in the reheating portions of these curves, where the 
magnetization attains a maximum at a temperature as low as —20°c, 
after which thermal demagnetization causes a decrease in intensity as 
the specimen regains room temperature. 

When a sample of massive hematite was ground into a powder, made 
up into a specimen in the way previously described, and subjected to the 
usual temperature cycle, the main features of the cooling curves of 
figs. 5(a) and 5 (b) were again observed ; the reheating curve, however, 
was somewhat different in that the amount of recovery was much greater 
than in either of figs. 5 (a) or 5(b), the final intensity being about 70%, 
of the intensity before cooling. 


§ 4.3 
Specular hematite was also examined in the low temperature in both 
bulk and powdered forms. With the bulk specimens magnetized in 
about 1000 oersted, the cooling produced a change of less than 7°, of the 
initial intensity. In the coarsely-powdered state, the specimen 
magnetized in a field ~6300 oersted showed a decrease in intensity of 
only about 10% at —75°c. So we conclude that specular hematite, 
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both in the bulk and ground-up form, shows almost no sign of the Morin 
transition. 
§ 4.4 
The results of the preceding sections indicate that whilst the two 


naturally occurring forms of hematite examined, that is, the massive and 
specular types, both have a rhombohedral structure, their magnetic 


Fig. 5 (a) 
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Specimen of massive hematite magnetized in a direction parallel 
to the fibre-lines. 


behaviour is quite different in the low temperature region. Specular 
hematite shows almost no transition effect, while the massive hematite 
shows a marked effect, which is, in fact, not very different from that 
found by Pauthenet for a single rhombohedral crystal from Elba. The 
Triassic red sandstones behaved in a similar way to the specular hematite 
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in that little or no transition effect was observed : however, the cause of 


the suppression of the transition effect may not be the same in both 
cases. 


In general, one must conclude from these results that the Morin 


transition must be very sensitive to the conditions of growth of the 
crystals. 


Fig. 5 (b) 
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Specimen of massive hematite magnetized in a direction perpendicular 
to the fibre-lines. 


§ 5.1 


To elucidate this fact, the order of magnitude of the size of the 
crystallites of massive and specular hematite was determined by X-ray 
photographs. It was found that the massive form of hematite consists 
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of very fine crystallites, probably less than 1p in size, whilst single 
crystallites of specular hematite have been found up to about 1mm in 
size, that is, a thousand times as large. It has been suggested that 
massive hematite occurs in regions where growth is partially inhibited 
by the pressure of surrounding material, thus restricting the size of 
crystals formed ; whereas specular hematite forms in rock cavities, etc., 
where growth is relatively free. This, however, is probably an over- 
simplified picture of the situation. The natural growth of crystals, and 
crystal habit, is a subject about which little is definitely known. 


§ 5.2 

The possibility that the observed difference m magnetic properties 
was simply due to a difference in crystallite size was examined. A sample 
of specular hematite was wet-ground in an agate mortar until the particle 
size was less than 1. A specimen of this material was then examined 
in the region of the transition, but no appreciable transition effect was 
observed. From this we may conclude that whilst a difference in crystal- 
lite size does exist between these two materials, it does not appear to be 
this difference alone which so markedly affects the transition. 


§ 5.3 


Spectroscopic analysis of the specular and massive hematite samples. 
yielded little new information. Both samples were found to contain 
0-7% Si and 0-1% Al. Whilst the massive hematite contained less than 
0-01°%, of either vanadium or titanium, the specular hematite contained 
0-06% V and 1-0% Ti. No other elements were detected in either sample. 
In view of our findings that quite large impurities of titanium do not 
suppress the transition materials in synthetic materials, it is unlikely 
that the difference between the magnetic properties of the two forms of 
natural hematite studied are due to chemical impurities. 


S08 


Throughout the course of these experiments it has become increasingly 
evident that the magnetic properties of hematite are extremely variable, 
being dependent not only upon the method of preparation of the sample. 
but also upon the past magnetic and thermal history of the sample. 
A further study was therefore made of the properties of various samples. 
of synthetic hematite powder. The three curves of figs. 6 and 7 show how 
the remanent magnetization of various samples decreases as the 
temperature is lowered. Curve 1 is for sample S,, which is a relatively 
poorly crystalline, commercial preparation of hematite of a brown-red 
colour ; curve 2 is for sample 8,, a high purity hematite rather purple 
in colour ; whilst curve 3 is for sample S,, a dark grey powder formed by 
heating S, in oxygen for 16 hours at 1100°c: this treatment normally 
increases the homogeneity as well as the crystal size of such a powder. 
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Whilst all the samples have been shown by x-ray powder photography 
to be purex Fe,0,, they exhibit widely different magnetic properties. 


§ 6.2 


For magnetizing fields ~17 000 oersted, the intensities of magnetization 
of 8, and S, are ~0-012 e.m.u./g and 0-135 e.m.u./g respectively : 


Fig. 6 
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Variation with temperature of remanent magnetization of three different 
samples of hematite (a Fe,Os). 


a difference of more than ten times. In spite of this large difference in 
saturation magnetization, the remanence coercivities at room temperature 
for the two samples are nearly the same, viz. ~4100 oersted and 
~4650 oersted for the above saturating field. For specimen 8; the 
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corresponding values for saturation magnetization and remanence 
coercivity are 0-15 e.m.u./g and 2400 oersted. Thus, in the heat 
treatment of S, to form S3, the coercivity has decreased by a factor of 2, 
whilst the saturation magnetization has remained practically unaltered. 
It would seem that heating samples S, at high temperature for any 
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Normalized cooling curves for three different samples of hematite (« Fe,0,). 


length of time tends to decrease the coercivity of the material without 
significantly altering the saturation remanence. Apart from the large 
difference of saturation magnetization at room temperature, the striking 
feature of the curves of fig. 6 is the large difference in the magnitude of 
the decrease in intensity as the specimens are cooled: The commercial 
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hematite, curve 1, fig. , shows a decrease to about one-half while the 
heat-treated pure specimen, curve 3, shows a decrease to one-thirtieth 
of the value at room temperature. 


Si lel 

The experiments decribed in the present paper bring out several 
important facts concerning the magnetic properties and behaviour of 
hematite. Firstly, naturally occurring hematites can be divided into 
two groups : those which show the Morrin transition on cooling, and those 
which do not. Of the samples examined in the present work, the greater 
proportion fall in the latter group—these being the crystalline specular 
hematite, and both the detrital and cement forms of hematite found in 
the natural sandstones. The only sample we have studied falling in the 
first group is the massive hematite ; however, Pauthenet showed that 
a single rhombohedral crystal from Elba shows the same effect. 

Whilst there is a difference in the crystallite size of our two samples, 
all the samples show the normal hematite lattice spacing when examined 
by X-ray powder methods. 

From the spectroscopic analyses of the specular and massive hematite, 
the only noticeable ion present in some and not in the other is the 1-:0°% 
of titanium in the specular material. Whilst this is one of the samples 
not showing the transition at —15°c, it is in marked opposition to the 
results of the chemically synthesized series of Ti in « Fe,O;. There is 
the possibility that the Ti ions in the natural material are more closely 
linked with the Fe ions, that is, in a more highly ordered state, than in 
the synthesized material, but the method of synthesis makes this unlikely. 
However, this does not explain the failure to reproduce Morin’s observa- 
tions, particularly when the samples examined here contained higher 
percentages of titanium than those used by Morin. 

Throughout this work the variability of the properties of hematite 
and the dependence of the samples on past history has been most 
noticeable, thus making any generalized statement of the expected 
properties of any hematite sample almost impossible. Hematite thus 
lends itself to phenomenological description, but not to a precise quanti- 
tative estimate of its properties. 
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ABSTRACT 


The entropy, the specific heat and the absolute temperature have been 
determined as a function of the magnetic temperature for a single crystal 
of neodymium ethylsulphate in the temperature range between 0-015 
and 1°K by adiabatic demagnetization experiments. The specific heat 
values were confirmed by determination of the magnetic specific heat 
between 1°k and 2°K by the paramagnetic relaxation method. The 
results are compared with those expected from paramagnetic resonance 
data and from Daniels’ extension of Van Vleck’s theory of dipolar inter- 
action effects. Good agreement is found down to temperatures of 
about 0-05°K. 


§1. INTRODUCTION 


In this paper an account is given of demagnetization experiments 
with a single crystal of neodymium ethylsulphate and the establishment 
of the absolute temperature scale down to 0-015°K. Results are also given 
of specific heat measurements by the paramagnetic relaxation method 
between 2°K and 1°K. This is part of a programme of investigation by 
different workers in this laboratory of the properties of several ethyl- 
sulphates at low temperatures. Unpublished measurements of the 
thermal and magnetic properties have so far been carried out on cerium, 
dysprosium and ytterbium salts above and below 1°k and on lanthanum, 
praseodymium, terbium and erbium salts above 1°K. 

The rare earth ethylsulphates are interesting for several reasons : 
first, since all these salts form isomorphous hexagonal crystals, a systematic 
investigation of the magnetic and thermal properties of the complete 
series from lanthanum to ytterbium is possible. This is not the case with 
the rare earth double nitrates, where only the salts of the elements from 
lanthanum to gadolinium crystallize. Secondly, extensive use has been 
made in this laboratory of several of the ethylsulphates as cooling agents 
for the experiments on nuclear alignment (Cacho ef al. 1955, Bishop et al. 
1957). For these experiments it is important to know the absolute 
temperature scale of the salts. Finally, it is of interest to compare the 
experimental results with those expected from paramagnetic resonance 
experiments (Cooke 1955) and from the theory of magnetic interactions of 
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Van Vleck extended by Daniels (1953). This comparison is done in the 
present paper for the case of neodymium ethylsulphate. It is shown in 
particular that the agreement with Daniels’ calculations of the suscepti- 
bility is satisfactory down to about 0-04°x and that the measured specific 
heat from demagnetization and relaxation experiments also agrees well 
with the predicted one down to about 0-05°K. 


§2. GENERAL PROPERTIES OF NEODYMIUM ETHYLSULPHATE 


Neodymium ethylsulphate has been the object of several investigations. 
Its crystal structure, together with that of several other ethylsulphates, 
has been determined by Ketelaar (1937). The chemical formula is 
Nd (C,H,SO,),.9H,O. The salt crystallizes in hexagonal prisms and has 
a trigonal axis of symmetry. All the paramagnetic ions are equivalent, 
and there is only one neodymium ion per unit cell. The natural 
concentration of 147Nd and 145Nd, each having a nuclear spin, is respec- 
tively 12:2°% and 8-3%, the remaining 79-5°, being made up by even 
isotopes. The ground state of Nd** is 4/4). and at temperatures below 
say 30°K, only the lowest doublet of this state is populated, the next one 
having an energy of 130 cm~! (Elliott and Stevens 1953). The experi- 
ments of Van den Handel and Hupse (1942) showed that at helium 
temperatures the magnetization of the Nd** ion in the ethylsulphate both 
parallel and perpendicular to the crystalline axis is represented exactly 
by the Brillouin function. This substance behaves therefore in this 
respect like an ideal paramagnetic. The results of paramagnetic resonance 
measurements with the field parallel and perpendicular to the crystal 
axis (Bleaney et al. 1954) showed that the spin Hamiltonian % of the 
neodymium ion can be written 

H = PS ,A +98 (S,H,+8,H,)+ KH vac : : (1) 
where A yuc is the Hamiltonian for the nuclear splitting of the ground 
state 

FO ng HAS DSL Si) ee ee ee) 
Here g,=3-54 and g,=2-07. S=} and I=$ for both odd isotopes. 
A and B have the values, expressed in cm~!, 1434 —0-0380, 143B=0-0199, 
1454 —().0236, 145B—0-0124, 

Roberts et al. (1953) measured above 1°K the paramagnetic susceptibility 
and the ‘ magnetic ’ specific heat by the paramagnetic relaxation method 
(Casimir and du Pré 1938). Their measurements were carried out on 
spherical powdered samples of different isotopic compositions. They 
were thus able to separate the specific heat into its nuclear and elec- 
tronic constituents. They found good agreement with the values 
calculated from paramagnetic resonance data for the nuclear specific 
heat while the electronic term was wholly accounted for by magnetic 
dipole interaction showing that exchange effects were negligible. For the 
sample of natural abundance the specific heat they obtained was given 
by C7*/R= 1:09 10-*. 
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Daniels (1953) extended Van Vleck’s treatment of the magnetic dipole 
interaction to the anisotropic crystals. He expressed the susceptibility 
and the entropy as a power series in 1/7’ and calculated the coefficients of 
the series as far as the term in 1/7'4 for ions with magnetic dipole inter- 
action in presence or absence of nuclear hyperfine structure (h.f.s.) 
interaction. He showed in particular that (1) the contributions to the 
susceptibility y and the entropy S from h.f.s. and dipolar interactions are 
additive as far as the term in 1/7", (2) that for the susceptibility the effect 
of the h.f.s. does not appear in the term in 1/7”, i.e. does not affect the 
Curie-Weiss 4 and that (3) 4 for a powder is not zero in general. Daniels 
gives detailed expressions for 4 in the case of spherical specimens for 
several salts, among them the ethylsulphates and a general formula for 
the higher terms. For a spherical specimen one has per gramme molecule 


of substance : 


m\ A EK 
Xsphere ~ 777 (1+ 7 == Wie +... :) Sr ing Bes (3) 


and for an ellipsoid the susceptibility is 


ees ES aentole (4) 
Xellipsoid — PzS TIE CAP) a= (e715 at a eS 


Ne 
5=5(F —s) 3 


A is the Curie constant of the substance, V the volume of one gramme 
molecule and « is the demagnetization coefficient of the ellipsoid. It is 
customary to express the magnetic properties below 1°K by the magnetic 
temperature 7'* defined as the ratio A/y for a spherical specimen (Kurti 
and Simon 1938). One has therefore in the case of an ellipsoid 


A 


Xellipsoid 
At sufficiently high temperatures the susceptibility is 


itll alk oth wee 
Xeupsoia = 7 | Tl ~ T—(A+8)° 


In the case of neodymium ethylsulphate, Daniels calculated for the 


where 


T= 


Be Sie sweat amma? Te ahh re) 


(6) 


parallel direction 4,=1310-*°K and for the perpendicular direction 
4,=—2%X10-3°x. In Appendix 1 of the present paper, the term FH is 
calculated for the case of crystals with axial symmetry, in particular for 


neodymium ethylsulphate. 


§3. THE SPECIMEN 


Neodymium ethylsulphate was prepared from Johnson Matthey 
‘Specpure’ Nd,O, and ‘Analar’ materials. An almost transparent 
single crystal was grown, which was then shaped to an ellipsoid weighing 
about 2-6 g. and having a length of 1-6 cm and a diameter of 1-3 cm, with 


2N2 
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the large axis parallel to the crystalline one. There were probably 
some water inclusions in the crystal, as the density was only 1-80 compared 
to the expected value of 1-84 (Ketelaar 1937). 


$4, THe APPARATUS FOR DEMAGNETIZATION EXPERIMENTS 


The specimen was suspended on eight cotton threads in a double walled 
brass demagnetization cryostat with its crystalline axis either parallel or 
perpendicular to the cryostat axis. The outer vacuum jacket sur- 
rounding the cryostat was made of pyrex and supported a mutual 
inductance coil with compensated secondaries to measure the suscepti- 
bility of the specimen. This unit was immersed in liquid helium and 
surrounded by a metal dewar of the type described by Henry and Dolocek 
(1950), using liquid air as external cooling agent. 1-5 litres of liquid 
helium were sufficient to cool the apparatus and to maintain the bath at 
4-2°x and the cryostat at temperatures between 4-2 and 0-85°K for 17 
hours. It was not found possible to measure the paramagnetic suscepti- 
bility of the specimen with a conventional audio-frequency a.c. bridge 
because the large amount of metal around the cryostat gave rise to a.c. 
losses and drifts in the mutual inductance of the coil, and the set-up had 
to be somewhat modified. Through an automatic contacter, pulses of 
current were sent through the primaries of the coil and a five decade 
standard compensating mutual inductance. The signal from the 
secondary circuit was amplified by a system of galvanometers and photo- 
cells (MacDonald 1947) and displayed by a luminous spot on a graduated 
screen. When the compensation was not exact, the spot moved across 
the screen at every change of the primary current by an amount pro- 
portional to the error in balance. The period of switching the primary 
current on and off was tuned to that of the amplifier, about 1-5 sec, 
so that the whole system was in resonance, thus reducing the irregular 
effects of thermal e.m.f.’s. The measuring field in the primary of the 
measuring coil was 23 gauss during the calibration and 2-3 gauss when the 
temperature of the specimen was below about 0-5°K. The sensitivity of 
the amplification system was sufficient to allow the detection of signal 
corresponding to a change of mutual inductance of 1-5 10-8henrys 
or a difference of susceptibility of 5-8 x 10-6 e.m.u./g. The magnetic field 
for the adiabatic demagnetization was provided by an iron-free solenoid 
of the type described by Daniels (1950). The maximum field was 
24 kilogauss corresponding to a power input of 200 kw to the solenoid. 


§5. GENERAL DESCRIPTION OF THE METHOD 


After the cryostat was cooled down to the temperature of liquid 
helium, the mutual inductance M of the measuring coil surrounding the 
specimen was determined at different temperatures between 4:2 and 
1-2°K. This calibration gave the relation between 7'* and the suscepti- 
bility. A series of isentropic demagnetizations was then carried out 
from known temperatures between 0-85 and 2-0°K and magnetic fields up 
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to 24 kilogauss. These experiments gave the entropies as a function of 

the magnetic temperature 7*. The measurement of the specific heat 

C*=dQ/dT* against 7’* was made by heating the specimen with y-rays. 

(e.g. Kurti and Simon 1935). The relation between the absolute and 

magnetic temperature was then obtained using the well-known expression. 

cr 

(ie Ty 

ERIEViGh ey en ee (7) 

and from this relation, the entropy and the specific heat were calculated 
as a function of the absolute temperature. 


§6. THE CALIBRATION OF THE MAGNETIC THERMOMETER 


Over the whole series of experiments four calibration curves, each with 
about 20 points, were taken with the crystal axis of the specimen parallel 
to the magnetic field and one with the axis perpendicular to the magnetic 
field. The demagnetization correction was calculated to be 52-2 x 10-3 
and 6,——0-7+10-*° in the parallel and in the perpendicular direction 
respectively. The temperature of the helium in the cryostat was 
determined from the 1955 helium vapour pressure scale, account being 
taken of the hydrostatic pressure effect and the thermomolecular pressure 
correction. In this temperature range, (M—M,), which is proportional to 
the susceptibility of the salt, was found to follow a Curie-Weiss law as. 
predicted by eqn. (6) with 4,—(14+7) x 10-%° in the parallel direction 
and 4,—(0+7) x 10-°° in the perpendicular direction, in agreement with 
the values predicted by Daniels. M, is the mutual inductance of the 
coil in absence of the specimen. The experimental splitting factors were 
found to be g,=3-51+0-05 and g,=2-05-+0-04, in agreement with the 
paramagnetic resonance data. 


§7. THE DETERMINATION OF THE ENTROPY AGarnst 7'* 


The isentropic demagnetizations from known fields H and temperatures. 
7 , were carried out in the usual way. The temperature 7’, was measured 
with the magnetic thermometer and its constancy during the magnetiz- 
ation was checked by measuring the vapour pressure of the helium in the 
cryostat with a McLeod gauge via a special manometer line. After the 
demagnetization, readings of the mutual inductance were started as 
soon as possible (15 seconds after the demagnetization) and continued 
every 10 seconds for about 2 minutes. The mutual inductance as plotted 
against time gave a straight line which was extrapolated back to the 
time of demagnetization. From eqn. (5), 7* was obtained. From the 
values of H and 7’, the entropy of a gramme molecule of an ideal para- 
magnetic salt S;,,,;, could be obtained using the tables of Hull and Hull 
(1941). To obtain the actual entropy, S. of the salt, S;,.,,; had however 
to be corrected for the effects of the lattice as well as of the dipole and 
h.f.s. interaction, these being to a good approximation respectively 
FT?/3 and D—(G/2T?) between about 0-25 and 4°x. Here F and G@ are 
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constants to be determined by experiment. D is the entropy of the 
degenerate system of the nuclear spins. As the total concentration of the 
odd isotopes is 20-5%, one has D=0-205 x log (2/-++1)=0-426. A term 
¢(H/T), which is calculated in Appendix 2, has also to be added to take 
account of the additional magnetic field in the crystal created by the 
magnetic dipoles. It follows that : 


S Sideal le fo ‘42 . 8 
Ro GRao 8) S72 ees one oes 
‘The specific heat in zero field is then 
OS orp tee x 
poll + m- (9) 


With the crystal axis parallel to the field, 10 isentropic demagnetizations 
from low fields and different temperatures were carried out in the tempera- 
ture region between 2° and 0-25° to determine / and@. Using the method 


He 


He in kilogcuss/ degree 


0.01 0.02 0.1 0.2 0.3 


The entropy S/R as a function of the magnetic temperature. A Crystal axis 
parallel to the field H. @ Crystal axis perpendicular to the field H. 
Columns on the right-hand side: H/7', for the crystal axis (A) parallel 
to the field (B) perpendicular to the field. 


of least squares, and assuming Curie—Weiss’s law 7=7* +14 10-8 to be 
valid down to 0-25° (see Appendix 1) these constants were found to be 
F=(2-10-4)x 10-3; G=(1-16+0-04-4 x 10-3. 


Figure 1 shows the experimental points of the entropy against the magnetic 
temperature for demagnetizations of the crystal with the axis parallel 


N eodymiwm Ethylsulphate at Low Temperatures 527 


and perpendicular to the field and the first three columns of table 1 
give smoothed values. It can be seen that the susceptibility in the 
parallel direction reaches a constant value for an entropy S/R=0-45 
corresponding to 7',*=0-0132. At temperatures below about 0-1°K, 
relaxation effects of the spin system, already described in detail by 
De Klerk (1951) and Beun et al. (1955) were observed. 


$8. THe DETERMINATION OF THE SPECIFIC HEAT AND OF THE ABSOLUTE 
TEMPERATURE SCALE 


The determination of the specific heat C* below 1°K was done with the 
crystal axis parallel to the field. The experiment gave O* in arbitrary 


Table 1. Neodymium ethylsulphate 


S * * os Cexp ov Cc Caan cale 
rie | os a aa R SS o> ga eres eke ar 
1-115 | 0-478 0-490 |0-50x10-7|} 0-492 |0:50x10-?} 1:63x10-% |0-46x 10-3 
1-110 | 0-285 0-300 1-25 0-300 1-25 3:15 1-23 
| 1-105 | 0-211 0-228 2-13 0-225 2-13 4:35 2-11 
1-100 | 0-169 0-188 3:03 0-183 3:03 5:46 3-03 
1-090 | 0-123 0-139 4-95 0-136 4-95 7-20 4-95 
1-070 | 0:0845 | 0-102 8-50 0-0965 8-50 9-60 8:50 
1:050 | 0-0685 | 0-085 10-5 0:0795 LOS 11-2 11-0 
1-025 | 0:0540 | 0-070 14-2 0-0645 14-0 13-0 14-2 
| 1-000 | 0:0450 | 0-0625 17-2 0-0550 16-7 14-5 17-0 
| 0-950 | 0-0330 | 0-0500 26°3 ():0426 25-0 17-2 20-5 
| 0-900 | 0-0267 | 0-0433 35:5 0-0370 31-5 18-5 23-0 
0-850 | 0-0228 | 0-0380 48-0) 0-0318 41-0 20:5 26-0 
0-800 | 0-0200 | 0-0341 63-0 0:0283 46-5 21-9 — 
0-750 | 0-O181 | 0-0320 79-0 0-0260 51-0 22-9 — 
| 0-700 | 0-0165 | 0-0305 96-0 0):0226 52°5 24.7 — 
0-650 | 0-0154 | 0-0298 | 113 0-0208 53-0 26-0 —- 
0-600 | 0-0148 | 0-0289 | 125 00-0195 54-0 26:8 — 
0-550 | 0-0143 oo 134 0-0185 55-0 — = 
0-500 | 0-0137 — 150 00-0167 —— — —- 
0-450 | 0-0133 - 165 0:0160 — — a= 
Remarks 


1. At ‘infinite’ temperatures, the entropy of the magnetic effects in zero 


field is S/R=1:119 
2. In column 7, the energy H/R is given, which has to be extracted if the salt 


is to be cooled from 1°K to the temperature 7’. 


units as a function of 7',*. From eqn. (7) the relation between 7’ and 
T,* was obtained in arbitrary units of 7. By fitting this curve at 0-5°K 
to the Curie-Weiss law 77 =7',*+ 14x 10-3, 7’ was obtained in degrees as 
a function of 7',* between 0-015 and 0-5°. Smoothed values are given in 
table 1 which shows that the Curie-Weiss law extends down to about 
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T=0-15°K. In this region C* is therefore equal to C=dQ/dT. Again, 
using eqn. (7), C* could be obtained in absolute units. 

Above 0-24°x, C obeyed eqn. (9) with G@=(1:14+0-04)x 107% and 
F=(2+1)x10-* in good agreement with the results obtained in $7. 
This shows that the use of the Curie-Weiss law makes the results con- 
sistent. The experimental points of C*/R are given in fig. 2. Combining 


Fig. 2 


02 04 06 08 10 
T * 
The magnetic specific heat O*/R as a function of 7*. 


the different experimental data obtained with the crystal axis parallel 
and perpendicular to the magnetic field, the relationship between S/R 
C*/R, C/R, T, T\* and T7,* was obtained. The smoothed results oe 
given in table 1 in the first six columns. The enthalpy derived from 
the measurements is given in column 7. Table 2 gives the absolute 
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temperatures reached after adiabatic demagnetization as a function 
OneEL Gs 

An independent check of the specific heat values below 1°K was made 
by carrying out measurements of the magnetic specific heat between 1°K 
and 2°k with the paramagnetic relaxation method (Casimir and du Pré 
1938) using the same crystal as before. The experiments were made at a 
frequency of 900 c.p.s. and in magnetic fields up to 380 gauss on an 
apparatus kindly lent by Dr. A. H. Cooke. The results can be seen in 
table 3. j 


Table 2. Absolute temperature reached after demagnetization as a 
function of the ratio H/T’ , in kilogauss/degree 


(1) Crystal axis parallel to the field H. 
(2) Crystal axis perpendicular to the field H. 


ol di) OY 1 2 3 4 8 10 15 20 25 


Or 


T(1) | 0-900 | 0-304 | 0-144 | 0-089 | 0-062 | 0-037 | 0-029 | 0-024 |0-0180 |0-0165 |0-015* 


7T(2) | 0-900 | 0-428 | 0-248 | 0-160 | 0-118 | 0-063 | 0-052 | 0-040 |0-0265 |0-0210 |0-018* 


* By extrapolation. 


Table 3. Magnetic specific heat of neodymium ethylsulphate above 1°K 


T CTR 
2-27°K (1-13, +0-015) x 10-3 
2-01 1:13, 
1:59 1-14, 
1-32 1:13, 
1-06 Lies 
| 


§9. DiscussION 


The measured specific heat of neodymium ethylsulphate may be 
compared with the value obtained from the paramagnetic resonance data 
of Bleaney et al. (1954). The magnetic specific heat at temperatures 
above about 0-25°K, which is the sum of the h.f.s. and the dipole inter- 
action contributions, is calculated from a series expansion of the partition 
function (Bleaney 1950) to be CT?/R=1-11,x10-%, in good agreement 
with all our experimental results, which therefore confirm Roberts's 
findings and conclusions. 

At lower temperatures, the nuclear specific heat has to be calculated 


from the exact partition function Z= > exp (—W,/kT’). The energies of 
j 
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the levels W,; are to be found in the paper by Bleaney et al. We then 
assume that the relation for the dipole interaction (Daniels loc. cit.) 
is valid down to 0-05°K. The total magnetic specific heat is given in 
table 1, column 8. The agreement with the experimental values is good 
down to about 0-05°K. 


Fig. 3 


Tee K 
0.12 


0.10 


0.08 


0.06 Parallel to the axis 


0.04 x 


——— _—s_ experimental 


calculated 


0.02 


ie) 0.02 0.04 0.06 0.08 0.10 0.12 


The magnetic temperature as a function of the absolute temperature in the 
parallel and in the perpendicular direction to the axis. 


In fig. 3, the experimental relations between magnetic and absolute 
temperatures in both directions are compared with the calculated ones 
(see Appendix 1). Over the whole temperature range where the power 
series converge, the agreement is within the experimental error, +6°% 
Thus for the parallel direction the comparison extends down © 0-04°K. 

In conclusion it may be said that neodymium ethylsulphate is 
suitable substance for demagnetization work and for testing the theories 
for the thermal and magnetic properties. It would be interesting to 
demagnetize samples of this salt with greatly enriched even neodymium 
isotopes as one would then be able to reach lower temperatures and stud 
the cooperative phenomena due to pure magnetic dipole intereeiten 
Lower temperatures could also be reached by demagnetizing crystals of 
neodymium ethylsulphate of natural isotopic composition, magnetically 
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diluted with lanthanum ethylsulphate. From the above calculations and 
experiments, it is easy to estimate fairly well the lowest temperature 
attaimable with a certain concentration of neodymium ions in the crystal 
and to calculate the approximate relation between 7’ and 7*. Thus one 
can expect to reach a temperature of about 8 x 10-°°K for a concentration 
of 20%, Nd®* and a temperature of 6 x 10-°°K for a concentration of 10°, 
Nd** starting from an H/T ,=25 kilogauss/degree. 
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aN aca Sad ch fl Bl BS. Coa 
The Evaluation of E in the Susceptibility Expansion Serves 
As shown by Daniels 
#H=EH,+#} 


where EL, is the effect of the pure dipole interaction and #1 is the effect of 
the h.f.s. 
Calculation of E 


The general expression for this calculation, too long to be reproduced 
here, is derived from the formulae 3, 8 and 10 in Daniels’ paper. For 
a crystal with the symmetry axis parallel to the z-axis of the reference 
system, one obtains : 


parallel direction : 


‘ ae 24 
Tae Lobe gt Oe 
Eou= 16k2 2 TP ri LF ipo tae 


perpendicular direction : 


4R4 1 372 2 l 4__ 192% j22?-++ 924 
n, = 28 |i y(4-3)\) - 33+ 
01 1 6k2 Ge NT ag ri3° a5 iy 


gu'g°P* 9 zt eee 
+ 6a 3| 27,37, 


ij ij 


where r,, is the distance between the ions 7 and j, and - is its component 
along the z-axis. The lattice sums 2M MZ DEM : . have been 
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evaluated by Daniels for the case of the ethylsulphates. These are : 
) 


y— =245-53 
j (5) 
2 
> ue — 144-32 
gas xa A6 
zt 


a9 


a2 1 \2 
( Ses: :) — 349-69 
Ga ae 


J 
where a is a certain dimension in the unit cell of the lattice. Thus for 
the ethylsulphates one obtains : 
Eyy=—a-® x 10-48 [4-42 X gy4+3°54 X Gn'gs"| 
Eov=—@ xX 10-48[5-72 xg," 1-16 X 9179171. 
Substituting the values g,=3-54, g,=2-07, a=13-992 A for neodymium 
ethylsulphate, one obtains : 
Eyy=—1-16 x 10-4 
Eyy=—2°64 x 10-°. 


Calculation of E+ 


In the case of an ion with hyperfine structure, it has been shown (Cooke 
et al. 1956) that 


2B 
By=—MI+]) = 
(424 BY 
ee Sere esas 
| Hyi=—I(I+1). 36 
For 148Nd, Eiyi=—7-0,+10-4 #,i=—1-64+10% 


while for 14°Nd, #,!=—2-74+10-4 H,1=—6-45+10-4 
1 is thus comparable to Ey. 
Adding H, to H! and taking into account the concentration of the 


isotopes in the natural neodymium one finally obtains for the susceptibility 
of neodymium ethylsulphate : 


mY P72 1°32 10% 2:24 10-2 
i enter he | re Ee AUS x I | 


S= Gia 7 7 an) /\> 
ae La 7. ff Te 
1a Ape eu s03 2S < 10-3. 2 S0e 10a 
Xia ies TT Le (eae 7. eee a 


APC Pe aNeD LX 82 
Evaluation of the term ¢(H/T) 


The treatment is the same for ions with or without h.f.s., since to a good 
approximation the h.f.s. does not contribute to a change of the entropy 
in a magnetic field. For the sake of simplicity, we are also not taking 
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account of the splitting of the energy levels by the dipole interactions in 


zero field eee it 

where M’ is the magnetic moment per cm? and « is a constant to be 
determined. We make the assumption that «M’<H, so that it can be 
treated as a perturbation. As a first approximation, we can put 

M'=N'(4g8) tanhx 
where N’ is the number of ions per em? and x=g8H/2kT assuming S=}. 
We omit here the indices ‘ parallel’ and ‘ perpendicular’ as this calcula- 
tion applies to any direction with respect to the crystalline axis. 
The partition function of an assembly of NV ions is then 


aN 2QR2 N 
hear (cosh e+ ATE tanh a sinh x) ‘ 
From the well known relation S=(0/07T)kT In Z one derives 


eNg*B? 
The sum of the first three terms is the entropy S/R,,.,; mentioned 
before, and the fourth term is due to the field of the magnetic dipoles. 
For small fields, this equation reduces to the one derived from Daniel’s 
series expansion. One obtains by identification: 4J=«N’g?8?/2k. One 
has finally S/R=S,,..,/R+¢(H/T) where 6(H/T)=(x4/T) tanh« sech? x. 
For the case of Neodymium ethylsulphate the effect of ¢ on the entropy 
S/Ris very small. It can easily be seen that 6<(4/7)(In 2—S,,,,,/R). 


a =In 2+-In cosh x—=z tanh x— x tanh x sech? x. 
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ABSTRACT 


Neodymium 147 has been aligned in a crystal of neodymium ethyl 
sulphate. The directional distributions and states of polarization of the 
90 and 530 kev y-rays have been measured. It is concluded that both 
transitions are mixtures of M1 and E2 radiation (6 E2/M1=—017 
+0-15 for the 90 kev and --0-75-L0-25 for the 530 kev y-rays) and that the 
530 kev state has spin 7/2+. From the 530 kev results a value for the 
nuclear magnetic moment of !47Nd of 0:44-0:06 nuclear magnetons is 
found. ; 


§ 1. INTRODUCTION 


THe nuclide neodymium 147 has been the subject of several investigations, 
including nuclear alignment experiments by Ambler et al. (1955, 1956). 
They obtained the alignment by cooling crystals of neodymium—cerium 
magnesium nitrate (Nd 1% Ce 99% ),.Mg;(NOs3)1.. 24 H,O containing some 
147Nd, the cerium being included for the purposes of magnetic cooling, 
and they measured the spatial distribution of y-radiation emitted in the 
radioactive decay. The neodymium ion has a large and anisotropic 
hyperfine structure coupling in this salt and the lowest magnetic state has 
I,=+4, hence at low temperatures the nuclear spins precess approximately 
in the plane perpendicular to the symmetry axis and they are said to 
be ‘ aligned’ in this plane (see e.g. Bleaney 1951). With this type of 
orientation the lower nuclear magnetic levels are compressed compared 
with those of a system where the nuclei are aligned along the axis, i.e. 
where /,—-+J state lies lowest, and for a given hyperfine structure 
interaction and temperature larger directional effects in the y-radiation 
will be observed for axial alignment than for planar ‘ alignment’. For 
these reasons we decided to align '47 Nd in neodymium ethyl sulphate, in 
which axial alignment is produced, and measurements were made of 
both the state of polarization and the directional distribution of the 


y-radiation in order to determine the dipole-quadrupole admixture in the 
y-ray transitions. 


{+ Communicated by the Authors. 
t Now at the Ecole Normale Supérieure, Paris. 
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§ 2> DISINTEGRATION SCHEME 


Kleven-day !4*Nd decays with the emission of two main B-groups which 
lead to excited states of 90 and 530 kev in the daughter nuclide 147Pm 
(see fig. 1). The comparative half-lives (log ft—7-4 and 7-0) of these two 
groups indentify them as first forbidden with a change in spin of 0 or 1 
and a change in parity. The absence of a £-group to the ground state 
suggests that such a transition would involve a spin change of at least 2. 
The half-life of the 90 kev state (2-44 10~® sec) led Graham and Bell 
(1953) to classify the transition of M1. Since both transitions lead to the 
ground state they must involve the same parity change and hence the 
530 y-ray must be M1, E2 or a mixture of these ; higher multipoles are 
excluded by lifetime measurements (<10~7 sec). The internal conversion 
data are in accord with these assignments. Many weaker y-rays are also 
present in the decay of !4’Nd but since their multipolarity assignments are 
uncertain, attention has been concentrated on the prominent y-rays in 
the nuclear orientation experiments. 


Fig. 1 
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The decay scheme of 147Nd. 


The daughter nucleus 147Pm decays to 147Sm the spin and moment of 
which identify it as 7/2~ (Murakawa 1954). From the first forbidden 
character of this £-transition an assignment of 5/2+ to the ground state 
of 147Pm may be made and the absence of a f-transition to this state 
suggests that the ground state of 147Nd is 9/2~. This leads to an assign- 
ment of 7/2+ to the 90 kev state and either 7/2+ or 9/2* to the 530 key state, 
as shown in fig. 1, and these have been assumed in the interpretation of 
the nuclear orientation results. 


§ 3. EXPERIMENTAL 
Bleaney, Scovil and Trenam (1954) have shown that the energies of the 
states of the Nd*+ ion in neodymium ethyl sulphate (diluted with lan- 
thanum ethyl sulphate) can be described by a Spin-Hamiltonian of the 
form . ; 
H = PH S,+9:8(4 S.thS,)+48 1,4 BS are ab) 
where gu=35 y= 352-072 S12 


536 G. R. Bishop et al. on the 
and for the stable isotope 1*Nd 


A=0-0547°K  B=0-0286°%K I[=7/2. 


! 


The magnetic moment of '48Nd is given by Bleaney (1955) as 1-03+-0-06 
nuclear magnetons. 

Since A is greater than B axial alignment of the neodymium nuclei 
will set in at temperatures of the order of A. Experiments have shown 
that the adiabatic demagnetization of concentrated neodymium ethyl 
sulphate from 25 kilogauss at 1°K lowers the temperature to about 
0-018°K and that therefore this is a favourable crystal for aligning neo- 
dymium isotopes. 

A number of crystals of neodymium ethyl sulphate Nd(C,H,;SO 4), . 91,0 
was grown from a solution which contained a trace of radioactive 14’Nd, 
prepared from ‘ Specpure ’ neodymium oxide irradiated in the Harwell 
pile. Five of these, weighing about 2 grams and containing 30 micro- 
curies of radioactivity were mounted in the low temperature cryostat, 
with their crystal axes horizontal. This is convenient for measurement 
of the y-ray angular distribution since it allows both counters to be 
mounted in the horizontal plane. The sample was demagnetized from 
0:9°K using fields up to 25 kilogauss, produced by a solenoid, applied perpen- 
dicularly to the crystal axes. The angular distributions and states of 
polarization of the two strong y-rays were measured at various tempera- 
tures between 0-018°K and 0-9°K. At the same time the magnetic tempera- 
ture 7'’*, defined by 7’*=A/y where x is the susceptibility of the specimen 
and A the Curie constant, was measured by the usual ballistic method 
after calibration in the helium range. A correction on account of the 
non-spherical shape of the specimen was made using the relation 
T®=T*+A, where A was 0-007°k. The absolute temperature was 
determined using the 7’'—7'® correlation of Meyer (1957). 

For measurements of the angular distribution of the y-radiation a pair 
of scintillation counters fitted with lin. cube Nal(T1) phosphors was set 
outside the cryostat around the specimen, one along the axis and the other 
on the equator at a distance of about 18 cm. The phosphors subtended 
at the sample angles of +.4° which are so small that no correction for the 
angular resolution of the detectors was required. The counters were 
connected through conventional amplifiers to a pair of discriminators 
whose settings were adjusted to straddle the photoelectric peak in the 
pulse height distribution and thus provide a coarse pulse height analysis. 
For the 90 kev y-ray the contribution to the counting rate due to the high 
energy radiation was estimated and allowed for using the counting rate 
observed on the high discriminator setting. The validity of this correc- 
tion was checked by using lead absorbers to eliminate the 90 kev radiation. 

The anisotropy « of the y-radiation is defined by 


ex lol (Oa) 2) 5 Se) oan) 


where /(@) is the intensity, normalized to unity at 0-9°x, of radiation 
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Fig. 2 
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emitted at an angle @ to the axis. Figure 2 shows the variation of « 
with 1/7 for the 90 and 530 kev y-rays. 

The y-ray polarimeter described by Bishop and Perez (1955) was set 
up near the cryostat so that the scatterer received radiation emitted in the 
equatorial plane where the polarization effect is a maximum. Measure- 
ments were made of N,, and N,, which are the counting rates for radiation 
scattered parallel to and perpendicular to the alignment axis. The 
effects of the 90 and 530 kev y-rays were separated by pulse height 
selection of the recoil electron pulse in the scattering phosphor. The 
variation with 1/7 of p’=N,/N,, is shown in fig. 3. The degree of polariza- 
tion p is defined by 

ewes 1/2) —I (2/2, 0) 3 
PT (x/2, 7/2) —1(nf2, 0) * i 


where (6, d) is the intensity of radiation at an angle @ to the axis with its 
electric vector inclined at an angle ¢ to the plane containing the axis and 
the direction of propagation. It is related to the measured p’ by 
the formula 

1—p’ R+1 


where R# is the analysing ratio of the polarimeter. 


§ 4. RESULTS 


For low degrees of alignment (A/7’<1) the intensity J(#,¢) can be 
expressed in the form 


I(0, 6)=1+B,(T) [A2P, (cos @)+A,’ cos 26P,? (cos 0)]. . (4) 


(see e.g. Cacho et al. 1955). B,(7') describes the orientation as a function 
of temperature. A, and A,’ are parameters determined by the dis- 
integration scheme and involve the spins of the levels and 54, the 
quadrupole-dipole amplitude mixing ratio. P, (cos @) and P,? (cos @) 
are Legendre and associated Legendre polynomials. 

The advantage of considering low degrees of alignment is that terms in 
P, (cos 6) and higher do not occur, and eqns. (1), (2) and (4) then give the 
following simple relation between polarization and anisotropy at any 
arbitrary temperature: p/e=2A,'/A,. This ratio is a function of 5, 
and tends to unity for a pure dipole y-ray (5=0). Simultaneous measure- 
ments of « and p’ therefore enable 5 to be determined. 

In the table the observed values of « and p’ are given for 1/7’=25. 
The analysing ratio R for the polarimeter was determined experimentally 
as 5-0-+-0-02 and 2-1+-0-2 for radiation of 360 and 935 key energy respec- 
tively. By interpolation a value of 3-5+1-0 for 530 kev radiation was 
found. As it is doubtful how R varies for energies below 360 kev it was 
assumed to lie between 2 and 10 for the 90 kev radiation. Column 4 gives 
the polarization p. The results show that a spin of 9/2+ for the 530 kev 
is not possible since this requires the transition to be pure E2, which would 
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give a positive sign for p whereas a negative sign was observed. The spin 
must therefore be 7/2 which fixes A,. In column 5 is given 5 determined 
from these values and from 6 and « the value of B,(7') for this temperature 
was calculated. 

In the evaluation of B,(7’) as a function of A/T exact calculation of the 
ionic level positions and the contributions to them by different nuclear 
states was made. From our experimental value of B,(7') at 1/7'=25 the 
hyperfine structure constant A for 147Nd was determined. Since AI/u 
must be a constant for all neodymium isotopes the magnetic moment of 
'47Nd may be determined and is shown in column 6. The solid curves in 
fig. 2 and 3 have been derived theoretically using these values for 5 and p. 


Summary of results for 530 and 90 kev y-rays. The observations on « 
and p’ were made at 1/7’'=25 deg-1. 


530 kev | +0-200+0-005 | 1:164+0-01 | —0-136 | +0-75+0-25 | 0-444 


90 kev | —0-093 0-005 | 1-054+0-005 | —0-055 | —0-17-+0-15 | 0-274 


§ 5. Discussion 


It will be noticed that the moment deduced from 90 kev results is smaller 
than that from the 530 kev results. This may be attributed to a perturba- 
tion of this state through coupling with extranuclear fields during its long 
lifetime (2-44 10~® sec). Effectively this reduces 6,(7') for this state 
relative to that of the 530 kev state by a factor A where A=0-5-0-2. 

Our results may be compared with those of Ambler et al. (1955) who 
found a 40% anisotropy in the 530 kev radiation but no significant 
anisotropy for the 90 kev y-ray. This lends support to our hypothesis 
that the 90 kev state is perturbed. In the light of the present experiments 
the absence of anisotropy cannot be caused by the H2/M1 admixture 
which they have suggested as the most likely explanation (Ambler e¢ ai. 
1956). In the results for the 530 kev y-ray there is also a discrepancy 
between the two sets of experiments. This can be seen by calculating 
the theoretical dependence of < on 1/7 in the double nitrate using the 
values for » and 6 obtained from the ethyl sulphate results ; one finds that 
the experimental values of « are up to a factor three times smaller at high 
temperatures giving rise to too low a value for the nuclear moment 
(0:22-0:05 mm). Although the lifetime of the 530 kev state has not been 
measured, for an M1/E2 transition of this energy it would be expected to 
be in the region of 10-11 second which is too short for a sizeable perturba- 
tion to occur, A more probable explanation for the discrepancy lies 
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therefore in some ionic effect, e.g. dipole-dipole interaction. The existence 
of this discrepancy does not however influence the conclusion about the 
mixing ratios for the y-ray transitions, since these are based upon simul- 
taneous measurements of « and p’ and do not depend upon B,(7). 

Since the possibility of reduction of the anisotropy of the 530 kev 
radiation cannot be excluded 0:44+0-06 nuclear magnetons must be 
regarded as a lower limit for the magnetic moment of 14’Nd. 


§ 6. CONCLUSIONS 


The 530 kev level of 147Pm has a spin 7/2+ and decays to the ground 
state with the emission of a mixture of M1 and E2 radiations in the 
proportions given by 6(E2/M1)=+0-75+0-25. The y-ray from the 
90 kev state (7/2+) is principally M1 with 6(EK2/M1)=—0-17-.0-05. 
Assuming that the spin of 147Nd is 9/2 the magnetic moment is 0-44-++ 
0-06 nuclear magnetons. The 90 kev level appears to be perturbed in 
both neodymium cerium magnesium nitrate and in neodymium ethyl 
sulphate. The attenuated angular distribution in the latter salt has 
A=0-5-40-2. 
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ABSTRACT 


Measurements of the thermal conductivity of single crystals of zinc 
and cadmium have been taken between 2 and 20°xK. The results show that 
a plot of W7' against 7? (where W is the thermal resistance) instead of 
being a straight line, has linear sections at high and at low temperatures, 
the slope at high temperatures being about three times that at low tempera- 
tures. This behaviour is shown to be connected with the very marked 
changes which occur in the values of the Debye 0. 


§ 1. INTRODUCTION 


THE theory of the thermal conductivity of a metal using a quasi-free 
electron model gives the sum of two terms as an expression for the total 
electronic thermal resistance W. These are W), the thermal resistance due 
to the scattering of the electrons by the impurities (this will be dependent 
on the particular sample of the metal being measured) and W ,, the thermal 
resistance due to the scattering of the electrons by the lattice vibrations 
(this should be a constant for any given metal). To a first approximation 
these two resistances may be considered as EEE and the total 
resistance is therefore. 
ree OSs © a ee Pe es 2 (1) 
For a discussion of the theory the reader is referred to Wilson (1953). 
The resistance, W,, should be of the form £/7', and W, should be of the 
form «7?. It is this latter term, the so-called ‘ideal’ resistance, with 
which this paper is concerned. 
From (1) we can write 
Wen DO EET: 3: "Mas, maseed (Ae ed tone (2) 
The theory suggests that (2) should be valid up to a temperature of the 
order of 6/10, although in practise it is found to hold up to ~6/5 (@ is the 
Debye characteristic temperature). The usual way of checking (2) is to 
plot WT against T?. This should give a straight line with a slope of «. 
Although the theoretical model is very idealized and might only be 
expected to hold for a metal such as sodium, it is found that for a very 
wide range of metals with complicated electron distributions eqn. (2) is a 
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good representation (e.g. see Rosenberg, 1955). Nevertheless earlier 
work indicated that in some cases (2) was not very accurate, 1.e. « was 
not constant at low temperatures. Hulm’s (1950) results on mercury 
show this very clearly. 

_The present work deals with a similar situation in measurements on 
cadmium and zinc and for these metals, at least, a satisfactory explanation 
can be presented in terms of the variation of the Debye @ at low tempera- 


tures. 
Fig. 1 


0 100073 2000 3000 


A plot of W7' against 7% for the cadmium crystal showing the breakaway from 
the straight line at low temperatures. The inset shows the low tempera- 
ture end of the graph on a larger scale and illustrates the fact that in this 
range the W7' curve becomes linear once more. Also shown in the inset 
is the variation of the Debye @ with temperature. 


§2. EXPERIMENTS AND SPECIMENS 


The thermal conductivity was measured between 2 and 20°k. The 
general technique using gas thermometers to measure the temperature 
difference has already been published (Rosenberg 1955). All the 
specimens were single crystals. Several cadmium specimens of different 
orientations were measured, all made from the same batch of Hilger 
VPS 8670, with a stated purity of 99-995°%. Except for a slight 
anisotropic effect, they all gave similar results and hence only one set 
will be presented in detail in this paper. One zinc crystal was measured. 


This was made from 99-997°, purity metal supplied by the Imperial 
Smelting Corporation, Avonmouth. 
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§ 3. RESULTS 


The results are plotted as curves of W7' against 7%. Those for cadmium 
(Cd5) are given in fig. 1 which shows that over the higher temperature 
range the WT’ curve is linear, but that it tends to decrease in slope at 


Fig. 2 


0 1000 a 2000 3000 


A plot of WT against T° for the zinc crystal. The inset shows the low tempera- 
ture section on a larger scale, together with the variation of the Debye @ 
over that range. 


lower temperatures. The inset to fig. 1 however, shows that up to 
T2—80 the WT curve is once again linear but its slope is only one quarter 
of that of the linear section which occurs at the higher temperatures. This 
means that we have two possible values of « for the high and low tempera- 
ture ends of the WT7' curve, which we can designate «, and «, respectively. 
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The values of these for the four cadmium crystals and the zinc crystal are 
given in the table together with the ratio ay/a,. The orientation of the 
hexagonal axis of the crystals to the crystal axis is also given. 


Specimen Orientation ay, x 10° ay X 10° oy Oey, 


The WT curve for the zinc crystal is shown in fig. 2. Similar behaviour 
to that in cadmium is observed. The graph is linear up to 7?=150 and 
it then curves upwards till its slope is nearly trebled when it again becomes 
a straight line. Due to the fact that the values of « are much smaller 
than those for cadmium (a consequence of the fact that @ is nearly twice 
as great for zinc as for cadmium) there is considerably more scatter in the 
results for zinc, particularly at the lowest temperatures. It will be noted 
that the total variation of WT over the entire «, region is only about 2%. 


§ 4. DIscUSSION 
The theory already referred to gives an expression for « in the following 
form 
GIN28 
“o= K, & oc eet eo . Ss (3) 
where JN is the number of electrons per atom, K ,, is the limiting electronic 
thermal conductivity at high temperatures and G is a constant with a 
value of about 70 (the precise value depends on the method used for 
solving the Boltzmann integral equation, see the reviews by Olsen and 
Rosenberg 1953 and by Klemens 1956). We are not here concerned with 
the numerical agreement between « as calculated from (3) and the experi- 
mental value except to observe that (3) presupposes a constant value for 
the Debye 6. Experiments however show that at low temperatures, 6 as 
determined from specific heat measurements, does not remain constant. 
The general pattern is illustrated in fig. 3 which shows the variation of 
6/0, for cadmium and zine (4 is the value of 6 at 0°K). There is a constant 
region at the lowest temperatures but as the temperature is increased, 0 
decreases quite rapidly, passes through a minimum, and then gradually 
increases and flattens off to a constant value. This drop in @ which is 
discussed in detail by Blackman (1956), is an inherent property of a 
crystal, once it is considered as a network of discrete particles, instead of 
as the elastic continuum of Debye. It follows from the fact that, on the 
particle model, the distribution function for the low frequency lattice 
vibrations rises at a faster rate than the (frequency)? which is used in the 
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Debye theory. This will therefore give a higher specific heat than the 
Debye theory, i.e. the effective @ will be smaller. If eqn. (3) is valid, 
therefore, it is not surprising that the slope of the WT curve changes and 
that there is more than one apparent value for «. It is therefore suggested 
that the increase in slope at the end of the «,, region is a reflection of the 
rather sharp decrease in @ which occurs at about the same temperature. 
The values of @ for cadmium (Smith and Wolcott 1956) and zine (Smith 
1955) are plotted on the W7' curves of figs. 1 and 2. If one takes as a 
criterion that a departure from linearity can be detected when the slope 


Fig. 3 


0 0.1 0.2 
T/g ; 
The dependence of the reduced Debye 6 on the reduced temperature for zinc 
(Smith 1955) and cadmium (Smith and Wolcott 1956). 6)—305°K 
for zine and 188°K for cadmium. 


of the WT curve (i.e. «) has changed by 10°, then if (3) is correct, this 
corresponds to a change in 6 of about 2-5°% +. Such a change from 6, has 
occurred by 3-5°K for cadmium and by 5°xK for zinc. These temperatures 
agree quite well with those at which the W7' curves cease to be linear— 
about 4°K for cadmium and 5°K for zinc. 

If one assumes that the low temperature, constant value of 6, 4), should 
be used for «,, and that «, corresponds to the region where @ has ceased 
to fall rapidly, i.e. around its minimum value, Omin, then from (3) we have 


eee on (4) 


+ Since K.. is proportional to 62, « is proportional to 6-4. [am grateful to 
Dr. P. G. Klemens for pointing this out. 
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From the work already quoted (6,/min)* for both cadmium and zine is 
about 5. This value is similar to that for the ratio a/a;, given 
in the table and is probably as good as one could expect because «; 
becomes constant over a range where @ is still changing although less 
rapidly than at lower temperatures. 

The above discussion shows quite clearly that the change of slope of the 
WT curves is bound up with the variation in the Debye 0. Nevertheless 
it should be emphasized that one must not expect complete agreement 
between « and @. In the first place it is still uncertain as to whether 
one should use the specific heat 6 in (3). Blackman (1951) has suggested 
that one should use a value, 6,, which is associated with the longi- 
tudinal lattice vibrations only, since the Bloch theory only assumes 
interaction with the longitudinal vibrations. Ziman (1954) also assumes 
this in his development of the theory. The experiments of Kemp eé¢ al. 
(1956) on the other hand indicate that the electrons do, in fact interact 
with lattice vibrations of all polarizations. However, whichever @ 
is used, it is probable that it would follow the same trend as the value 
derived from specific heat measurements. 

There is, however, a more important aspect which is sometimes over- 
looked. The 0, as it appears in the conductivity theory, is not an 
expression of the specific heat as such. It arises from the calculation of the 
interaction of the electrons with the lattice vibrations and the averaging 
of this interaction over all lattice frequencies. The way in which this is 
done so as to yield a useful solution is to introduce the Debye 6, which is 
also obtained from an average over all lattice frequencies. The method 
of averaging used, however, is different in the two cases and it is unlikely 
that the results of these two types of averaging should give exactly the 
same solutions. Therefore whilst a change in one average will almost 
certainly correspond to a change in the other, it is unlikely that there will 
be exact agreement. Hence, whilst these changes in « which have been 
described seem to correspond quite well with the observed variation of 
specific heat 0, it is not surprising that exact agreement is not found for the 
values of a ;/;. 

Although this paper has been concerned with cadmium and zine where 
the change in @ is quite large, it seems probable that the increase in « 
should occur to a greater or lesser extent in all metals, depending on how 
great a variation there is in 6. For many metals this variation is quite 
small, but there are a few other cases where an effect should be observed. 
For mercury (Smith and Wolcott 1956) 6 changes from 77° at 1°x to 
50° at 3°K. This gives a value of (4)/@min)* of 5-7 which indicates that « 
should change considerably in this region. This is undoubtedly the cause 
of the non-linearity of the W7' verses T? curves given by Hulm (1950). 
Since the constant value of « in this case would be below 1°K he did not 
find the initial linear part of the curve (corresponding to «,) and he begins 
with the curved section. His straight section, however, «, starts at about 
2-5°K which is just where 6 stops changing rapidly. Another case is that 
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of gallium where (4)/@min)* is 10 (Wolcott 1955), and the thermal conduc- 
tivity results indicate that «/a,, is about 2 (see Rosenberg, 1955). 
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ABSTRACT 


Specimens of medium carbon steel subjected to slow and dynamic 
compression have been examined by optical and electron microscopy 
In slow deformation, coarse slip in ferrite occurs in very narrow regions 
on corrugated surfaces; this behaviour appears to be related to the 
presence of an ageing precipitate. Coarse slip is absent in rapidly 
strained steel ; the metal deforms by fine slip which occasionally forms 
deformation bands. Deformation twins were only observed on one 
specimen of which the microstructure was abnormal. 


§ 1. INTRODUCTION 


OBSERVATIONS on the mechanism of deformation during slow and impact 
compression of mild steel have recently been described by Campbell and 
Duby (1956). It was found that the coarse, wavy slip lines which 
characterize plastic deformation of the « iron at normal strain rates are 
absent after dynamic straining, and that impact seems to cause less work 
hardening than slow compression. Some further metallographic observa- 
tions on a medium carbon steel, made by both optical and electron 
microscopy, are here reported. 


§ 2. EXPERIMENTAL METHODS AND RESULTS 
The steel used was a medium carbon steel (M2185) of the following 
composition : 
C Mn Si 8 BE 
per cent 0-32 0-62 0 235 0.037 0-025 


The steel was supplied in the form of 2} in. diameter bars. In order to 
avoid any segregated impurities, specimens 3 in. in diameter were cut 
from the outer layers of the rod. These were stress-relieved for one hour 
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at 600°C, and then machined and ground to their final dimensions : 
0-499 in. diameter, length 0-5 in. They were finally annealed at 850°c for 
30 minutes in vacuo and slowly cooled. For metallographic examination 
flats were filed on the sides of specimens ; these were ground, electro- 
polished and etched, before test. The method of deformation was 
identical with that previously described by Campbell and Duby. 

The first series of specimens to be examined consisted of two cylinders 
which had been strained to about 4°%, one in slow compression and one 
in impact (at an average plastic strain rate of 500 sec), and a third 
which had been twice impacted to a strain of 8%. 

The microscopic examination was made optically with a light profile, 
and carbon replicas of the surface for electron microscopy were prepared 
by a method previously used by Puttick (1957). 

The light profile showed that the coarse slip steps developed on the 
slowly compressed cylinder varied widely in height, up to a maximum of 
about 0-5 w as in fig. 1. (The vertical magnification is about twice the 
lateral magnification.) As might be expected, little sign of deformation 
could be seen in pearlitic areas, but large shears were observed at the 
boundary between pearlite and ferrite, indicating that such areas can be 
regarded at these strains as elastic inclusions in a matrix of plastic ferrite. 

The oblique illumination required for the light profile produces a 
schlieren effect which shows up, rather more clearly than normal illumina- 
tion, fine slip lines in the rapidly deformed specimen (though they cannot 
be separately resolved). It is interesting to note that under these 
conditions, deformation bands are occasionally formed (fig. 2). If Mott’s 
(1953) suggestion is correct, that deformation bands are only formed by 
fine slip, the fact that bands are not observed on slowly strained « iron 
suggests that there is little or no fine slip between the coarse lines 
characteristic of slip in body-centred cubic metals. 

The ferrite had in places a mottled appearance, such as often indicates 
the presence of a sub-microscopic precipitate. A section was therefore 
cut from one of the specimens, electrolytically polished, and anodically 
etched in 10% chromic acid. This preparation has been found by 
Puttick (1957) to reveal fine cementite particles clearly in steel, and the 
‘rough-cast ’ texture of many of the ferrite grains as in fig. 3 confirms that 
some kind of ageing has occurred in the « iron, probably during the 
furnace cooling. This section revealed a second feature worth noting, 
exemplified by the bright vertical line in fig. 3. As the illumination is 
oblique, this indicates a surface step, which on a deformed surface would 
naturally be attributed to shear of the ferrite between two adjacent 
pearlite areas. However, the metal had not been further deformed after 
polishing, and it is therefore concluded that the electrolytic polishing or 
etching has revealed some details of the previous deformation, perhaps 
due to strain-ageing segregation of carbon atoms to dislocations. 

The condition of the ferrite and its mechanism of deformation were 
further elucidated by electron microscopy. Figure 4 shows a replica of 
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a slowly deformed ferrite grain. The points of interest are : 

(i) The surface texture, confirming the existence of small etch-resistant 
particles. 

(ii) The coarse, forked slip lines. 

(iii) Slip apparently beginning on a second system, transverse to the 
large slip lines ; the lines are faint but seem to have a definite and fairly 
regular spacing. 

The nature of the small particles observed in steel has recently 
been the subject of some controversy ; it has been asserted that they 
may not be true precipitates, but regions of high carbon concentration 
round dislocations. Such segregation could certainly lead to unusual 
chemical effects, as has been suggested above, but we doubt whether it 
is the true explanation of the surface texture of ferrite in the present case. 
For one thing, the particles were almost certainly present before deforma- 
tion, yet their density is in many places 1011—1012/cm? (for instance fig. 7). 
This is very much higher than current estimates of numbers of disloca- 
tions in annealed metals. Also, the density of the particles is not markedly 
higher at grain boundaries, to which carbon atoms must certainly migrate 
even more readily than to single dislocations. On the whole, the surface 
appearance is most simply explained by homogeneous precipitation of a 
second phase, probably cementite, in the ferrite. After all, if such 
precipitation does not occur the ferrite formed in mild steel during 
cooling after austenitizing must be heavily supersaturated at room 
temperature. 

The surface texture of ferrite grains makes observation of fine slip in 
the impacted specimens difficult, except at boundaries between grains, 
or between ferrite and pearlite, where the requirements of continuity may 
impose severe local distortion. This is, of course, the explanation of the 
fact that deformation of rapidly strained metals apparently occurs at 
grain boundaries rather than within the grains ; fine slip cannot be 
resolved, and only the discontinuities at grain boundaries are visible. 
It is quite a different process from ‘ grain boundary slip ’ which accommo- 
dates deformation at high temperatures and low rates of strain. This 
point is emphasized by fig. 5 in which slip steps can be seen following the 
boundaries of a severely deformed grain. 

Some interesting characteristics of coarse slip in « iron can be discerned 
in fig. 6 (also of the slowly compressed specimen). This micrograph is 
actually one member of a stereoscopic pair, measurement of which shows 
that the total shear at the large slip step there visible is about 0-4 ie 
This slip seems to have taken place in a very narrow region ; there is no 
suggestion of a slip band resolvable into an aggregate of lamellae as in 
face-centred cubic metals such as aluminium. The step does, however, 
exhibit a number of striations which may reasonably be taken to indicate 
the slip direction. It has been suggested (Maddin and Chen 1954) that 
‘wavy slip’ in iron on apparently irrational planes is really the result 
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of intimate cross slip on elements of low index planes (such a process 
must indeed be the explanation of the difference between slip in face- 
centred cubic and body-centred cubic metals) but such cross slip has never 
actually been demonstrated. Here, and in similar photographs such as 
fig. 7, the scale of the cross slip seems to be related to the size of the 
precipitate particles. This seems a reasonable supposition, since if 
dislocations are held up by the precipitate particles as in Mott and 
Nabarro’s (1948) model of slip in age-hardened alloys, dislocations in 
screw orientation would climb out of the plane of their source to surmount 
the barriers. It need not, of course, be supposed that a single dislocation 
source is responsible for the large amount of slip in figs. 6 and 7; since 
slip lines can branch and coalesce, the steps may well represent the total 
slip from a number of sources. 

It is interesting to compare slip in this material with that in « iron 
which is apparently free from precipitate. Figure 8 is an electron micro- 
graph of deformed free ferrite in a mainly pearlitic carbon steel, the 
composition, heat treatment, and deformation of which has been described 
elsewhere by Puttick (1957). The specimen was strained about 3% in 
tension. The slip lines in this case are rounded and diffuse, in marked 
contrast to the sharp ‘ shear cliffs’ in figs. 6 and 7. Although the two 
steels are of somewhat different composition, it does seem possible that 
the coarse slip in figs. 6 and 7 is associated with the presence of fine 
cementite particles7. 

An attempt was made to prepare specimens of the 0-32°%, carbon 
steel which would be completely free of precipitate. ‘Two cylinders were 
austenitized in argon at 850°C, transformed at approximately 710°c 
in a lead bath for two hours, and air cooled. Examination of these 
specimens showed that the attempt was not entirely successful, but 
revealed some other features of interest. The microstructure of one 
specimen was similar to those previously examined, except that the 
pearlite was coarser and the precipitate was less apparent ; but the other 
specimen had apparently suffered some surface decarburization, possibly 
the result of inadequate flow of argon through the furnace tube. In this 
specimen the grain size at the surface was much larger than in the 
other, and the amount of pearlite very small. Some precipitate within 
the ferrite was visible with the optical microscope, not randomly distri- 
buted as in the previous specimens, but in the form of a network. 
Electron microscopy showed that cementite had precipitated along many 
grain boundaries in the form of films (fig. 9) and as films or rows of fine 
particles along lines within the « iron (figs. 10 and 11). Some of these 
form a continuous network, and presumably delineate dislocation walls 
forming sub-boundaries ; other rows of particles seem, however, to be 
completely isolated, and may indicate the intersection of single disloca- 
tions with the surface. 
+The difference in slip behaviour cannot be attributed, as it has been in 
other cases, to difference in surface preparation ; both steels had been electro- 
lytically polished before etching. 
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Both these specimens were subjected to impact. That with normal 
microstructure behaved in the expected way, deforming by fine slip in 
the ferrite. The surface decarburized specimen, however, though showing 
evidence of fine slip, also exhibited a number of twins (figs. 12 and 13). 

The formation of twin lamellae by impact of iron has often been 
reported (no doubt this is the origin of those first observed by Neumann 
in meteorites) but in the present tests they were only seen on this one 
specimen. The microstructural features suggest that the difference in 
behaviour is to be ascribed either to the heterogeneous precipitation of 
cementite within and between the ferrite grains, or to the grain size 
itself. It might be argued, for instance, that precipitation on dislocations 
has so hindered their movement, and grain boundary cementite so impeded 
the spread of plastic flow from grain to grain, that the critical shear stress 
for slip has been raised to a value at which twinning is an alternative 
possibility. The existence of a grain size effect in twinning has also been 
reported, however (e.g. Clark and Craig 1952), so that this factor cannot 
be disregarded. The evidence is rather vague, but does suggest that the 
critical shear stress for twinning may be a decreasing function of grain 
size. No explanation of such an effect in terms of current theories of 
deformation twinning, such as the pole mechanism of Cottrell and Bilby 
(1951), seems in sight at the moment. Further experiments on the subject 
are desirable. 


§3. CONCLUSIONS 


Metallographic examination of compressed medium carbon steel 
specimens has revealed the following points of interest : 


1. In slow deformation, large amounts of slip occur in very narrow 
regions on corrugated surfaces. This behaviour seems to be related to 
the presence in ferrite of an ageing precipitate, probably cementite. 


2. Fine slip in rapidly deformed specimens occasionally forms deforma- 
tion bands, which have not been observed in slowly compressed specimens. 


3. The strain is not homogeneous ; up to the strains of these experi- 
ments, pearlitic areas undergo little deformation. The apparent ‘ grain 
boundary movement ’ previously observed by Campbell and Duby (1956) 
is a localized shear imposed by the requirements of continuity. 


4. Deformation twinning has only been observed in one specimen, 
of which the microstructure was somewhat abnormal. This behaviour 
can be attributed either to raising of the critical shear stress for slip by 
precipitation of cementite on dislocations, or to lowering of the critical 
stress for twinning by an increase in grain size. 
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ABSTRACT 


In the alloys MngjAlygCo) and Mng Zn g9C a) the Bohr magneton number 
per manganese atom for the ferromagnetic state differs by 0-32, being 
greater in the zinc alloy. It was thought that this difference could be 
accounted for by the difference in valencies between Al and Zn, the extra 
valency electron of the aluminium atom over that of zine causing a 
reduction of one Bohr magneton per unit cell, corresponding to 0-33 Bohr 
magneton per manganese atom. If this were the case then a simple 
variation of Bohr magneton number might be expected in the alloy system 
MngpAl Zn 99 C29 

In addition, since the alloy MngyZn Cy) possesses a transformation at 
—42°c (when it becomes an antiferromagnetic substance with a complex 
magnetic structure) whereas the aluminium alloy shows no such transition, 
measurements on the quaternary alloys should illustrate how this transi- 
tion is dependent on the electron concentration of the corner atom in the 
unit cell. 

The results show a complex behaviour which points against any simple 
valency mechanism controlling the ferromagnetic properties. Although 
the lattice parameter of the quaternary alloys varies linearly with Zn—Al 
content, both the ferromagnetic and antiferromagnetic moments are 
reduced when aluminium partially replaces zinc in MnggZng Coo. The 
ferromagnetic—antiferromagnetic transition moves to lower temperatures 
and attains 0°K for the composition Mn, jAl;.;Zn,4.5Co9. For aluminium 
contents greater than 5-5°%, the alloys show only ferromagnetic properties 
and the saturation moment then increases with increasing Al content. 
However, at 10% Al the alloys contain two phases and further magnetic 
measurements are unrewarding. 

Similar measurements over a more restricted range of composition 
were made on alloys of the form MngyGa,Znoy_,Co9 and gave results in 
agreement with those for the Al-Zn alloys. 


+ Communicated by Professor W. Sucksmith, F.R.S. 
t Now at Atomic Energy of Canada Ltd., Chalk River, Ontario, 
§ Now at AB Atomenergi, Stockholm, Sweden, 
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§1. Iyrropuction 


THE alloys of compositiont MnggAly5Cy) and MngpZno Co) have many 
similar features. In particular they are both spontaneously magnetized 
at room temperature and are isomorphous, the structure being highly 
ordered face-centred cubic with the manganese atoms at face-centres, 
carbon atoms at body-centres and the remaining atom type located at 
the cube corner sites. The properties of these alloys have been described 
by Butters and Myers (1955 a,b). In spite of many similarities there is 
one rather striking difference ; whereas the alloy MngjAljoCo9 is ferro- 
magnetic down to 4°K, the lowest temperature it has been studied 
(Brockhouse 1956), MngoZny9Co) has a second order transition at —42°c, 
Below this temperature the structure becomes tetragonal and the magneti- 
zation decreases with decreasing temperature. This transition has been 
studied in detail, Brockhouse and Myers (1956), and it has been shown that 
below the transition temperature the alloy is antiferromagnetic possessing 
a complex magnetic structure. Above —42°c the alloy is a normal ferro- 
magnetic material and extrapolation of the saturation magnetization— 
temperature curve (o—7’ curve) together with neutron diffraction data 
shows the Bohr magneton number per manganese atom to be 1-58 at 0°K. 
The o—T curves for Mn—Al-C alloys show the alloy Mn, Al,,Cy) to have a 
Bohr magneton number per manganese atom of 1-26. Thus the Bohr 
magneton numbers of these two alloys differ by 0-32 per manganese atom. 

It is to be noted that this difference could be accounted for -by a valency 
effect since the aluminium alloy contains one electron per unit cell more 
than the zinc alloy. If this extra electron can be effective magnetically 
then we have a possible reason for the difference in magneton number 
between the two alloys because it would decrease by one the magneton 
value per cell, i.e. 0-33 per manganese atom and this is very close to the 
actual difference observed. If this were the case then we should expect 
to find a simple variation of Bohr magneton number with composition 
in the system Mn, Al, Znoo_ .Cao- 

Such a picture involving a valency conditioned magnetization is 
suited to a collective electron ferromagnetism for the alloys. The rather 
small non-integral magneton values favour this, but on the other hand, the 
highly ordered structure, which is stable at high temperatures and 
believed to persist up to the melting point, suggests that direct bonding 
occurs between certain atoms, in particular the manganese and carbon 
atoms ; this is difficult to reconcile with a collective electron model. 

A further point is the presence of the transition in the zinc alloy and 
the antiferromagnetism which involves at least three different types of 
manganese atom distinguished by different and probably integral magnetic 
moments. Measurements on alloys having the form Mng Al,Zngo_ Coo 
should illustrate how this transition is dependent upon the valency of 
the atom at the cube corner site. 
ee a 

+ All compositions are referred to atomic percentages. 
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During the course of this work it was found that Mn-Ga-C alloys, 
isomorphous with those of the Mn-Al-C and Mn-Zn-C systems, also 
exist, but have only ferromagnetic properties. Alloys of the form 
MngpGa,ZN oC were also studied and their properties are described 
here. 


§2, EXPERIMENTAL PROCEDURE 


The quaternary alloys containing zinc were produced by mixing zinc 
fillings with crushed powder from master ternary alloys containing Mn, 
Aland C or Mn, Ca and C and sintering the mixtures at 600°C in evacuated 
quartz tubes. The product was repeatedly crushed and resintered until 
x-ray powder photographs showed sharp diffraction lines. 

X-ray powder camera photographs were taken as a routine measure. 
Saturation magnetization values were obtained on coarsely crushed 
specimens, a magnetic field of 16 200 oersteds and a Sucksmith ring balance 
being used. It was established that the materials were saturated in this 
field strength. Measurements were made from liquid nitrogen tempera- 
ture to the Curie temperatures of the alloys (usually less than 250°c). 


§3. RESULTS 


_ In both the Al and Ga quaternary alloy systems the highly ordered 
face-centred cubic structure did not occur for all values of x in the 
composition formula; in particular, alloys containing between 10 and 
15% Al, and alloys containing more than 7°% Ga consisted of two phases. 
Magnetic measurements were made only on single phase alloys and of these 
those containing more than 16°, Zn showed the maximum in the o—-T' 
curve that is characteristic of the alloy MngZng Co. When such a 
maximum exists, provided it occurs at a temperature somewhat higher 
than —190°c, two values of Bohr magneton number are calculable. One 
value of magnetization appropriate to 0°K is obtained from the ferro- 
magnetic range being that part of the curve determined above the 
transition temperature (the temperature at which the magnetization has 
its maximum value); this portion of the o—7' curve is extrapolated to 
0°K and a value of magnetization and thereby a Bohr magneton number 
determined. The extrapolation is, on account of the steepness of the 
curve, somewhat uncertain but definite upper and lower limits to the 
magnetization can be set. This value we choose to call ordinary. On 
the other hand at temperatures below that of the transition the magnetiza- 
tion varies more slowly with temperature and a quite accurate extrapola- 
tion to 0°K can be made. Since, by analogy with the behaviour of 
MngoZng9V9, the material is probably in a complex antiferromagnetic 
state in this temperature range we choose to call the magneton number so 
determined the ‘ extraordinary ’ value. 

The quaternary alloys were assayed for all components. Occasionally 
wide departures between assayed and intended compositions were 
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Fig. 1 


Bohr Magneton Value 
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Variation of Bohr magneton numbers with composition in Mn—Al-Zn-C alloys. 
O ordinary value; ( extraordinary value. 


Fig. 2 
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Variation of Bohr magneton numbers with composition in Mn—Ga-Zn-C alloys. 
© ordinary value; (J extraordinary value. 
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found, in such cases the alloys were discarded. Generally the assayed 
compositions were in good agreement with the intended compositions but, 
as might be expected with alloys of this kind, the total Al and Zn, and 
total Ga and Zn were not always exactly 20°, of the whole. In these 


Fig. 3 Fig. 4 
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Variation of transition temperature with Al and Ga content. 


Fig. 5 
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Curie temperature and composition in Mn—Al-Zn-C alloys. 


cases the amounts present were scaled to 20%, and in the graphs generous 
tolerances have been allowed as shown by the horizontal lines through the 
points. On account of the necessity for extrapolation and the difficulties 
involved in controlling the compositions of these quaternary alloys we 
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treat our data as semiquantitative but believe that the trend in properties 
shown by our results is in general correct and to have significance. 

The ordinary and extraordinary Bohr magneton values per manganese 
atom for the two sets of alloys are plotted against composition in figs. 1 
and 2. Both systems show very similar results although the data for the 
Ga series are more restricted than for the Al alloys. The ordinary and 
extraordinary magneton values are reduced, but the former more rapidly 
than the latter when Al or Ga replaces Zn in the alloy Ming Mls, on: 
The transition temperature moves to lower temperatures as the zinc 
content is reduced, figs. 3 and 4. The variation of magneton numbers 


Fig. 6 
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Variation of Bohr magneton numbers with Mn content in Mn—Zn-C alloys 
© ordinary value; A extraordinary value. 


with composition allows us to determine that composition for which the 
transition reaches 0°K for in this case both ordinary and extraordinary 
values are the same. This occurs at about 5-5°% Al or Ga. 

For aluminium contents greater than 5-5° the alloys are normal ferro- 
magnetic materials, the magneton value then increases with increasing 
aluminium content until at about 10° Al the alloys contain two phases. 
At the aluminium rich end of the system single phase alloys are again 
found and the variation of magneton number seems to match very well 
the curve in the region of 10°, Al. The variation of Curie temperature 
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for the Al-Zn alloys is shown in fig. 5. Room temperature lattice para- 
meter measurements on these alloys showed the parameter to vary 
continuously and in a linear fashion. Similar data for the Ga—Zn alloys 
is limited on account of instability in this system but it is in agreement 
with that obtained for the Zn—Al alloys. 


§4, DISCUSSION 


It is clear that the magnetic properties of the ternary alloys Mng Aloo a9 
and MngpZnyoCo9 are very dependent upon the electron concentration of 
the corner atoms in the cell, but this dependence is not due to the simple 
valency mechanism proposed in the earlier section of this paper. 

It has been mentioned that as Al or Ga replaces Zn in Mng)Zng Co both 
ordinary and extraordinary magneton numbers are decreased. This 
result may be compared with the behaviour of alloys of the form 
MngyMn,Zno_,C¥o9 Where manganese replaces zinc, Butters and Myers 
(1955 a). In this case the ordinary magneton value is decreased whereas 
the extraordinary value is increased, fig. 6. The transition temperature 
moves to lower temperatures as in the other alloys. 

These alloys show complex behaviour and we believe that much more 
data will be required for a proper understanding of their properties. In 
particular neutron diffraction measurements should prove valuable but 
the preparation of large homogeneous samples of quaternary alloys 
containing both highly oxidizable and volatile components offers con- 
siderable difficulties. 
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SUMMARY 


Results of experiments in connection with climb of dislocations in 
fluorite crystals are given. Peculiarities of helical dislocations are 
explained. 

Observations concerning rows of closed loops and special features are 
described and discussed. 

Double rows of etchpits were identified with the emergence points of 
the parts of helical dislocations obtained by cleaving through helices. 


§ 1. INTRODUCTION 

ULTRAMICROSCOPIC examination of natural fluorite crystals, which had 
been additively coloured with sodium at 700°c, proved that dislocation 
lines were decorated (W. Bontinck and W. Dekeyser 1956); irregular 
networks and dislocation walls which resulted from slip or climb were 
observed. By treating synthetic CaF, crystals at much higher tempera- 
tures (1200°c) by a method described below, decorated helical dislocations 
have been found as reported in a previous note (W. Bontinck and 
S. Amelinckx 1957). Their relation to whisker growth as well as the 
nature of the colour centres in CaF’, has been discussed in a paper by 
S. Amelinckx et al. (1957). In this paper, additional observations and 
properties of helical dislocations are given. 


§ 2. EXPERIMENTAL PROCEDURE 

The synthetic crystals of CaF, (Harshaw) were first coated with a thin 
silver layer evaporated in a high vacuum, and afterwards heated for 
several hours in a hydrogen atmosphere at respectively 1200° and 1100°c. 
The role of the silver is not clear at present and is still under investigation. 
It turns out that the silver layer prevents a too violent reduction of the 
crystal by the hydrogen; the underlying mechanism is however not 
clearly established. The crystals were cooled very slowly in order to 
enhance colloid formation. 

The decorated patterns were observed by ultramicroscopy, after removal 
of the surface layer which contained large amounts of colloidal specks. 
The observation plane was the (111) plane, which is also the cleavage 


plane. 
+ Communicated by the Author. 
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§3. ResuLTs AND DIscUSSION 
3.1. Properties and Peculiarities of Helical Dislocations 


As shown in a previous paper, dislocations with a mixed character 
can climb into regularly spaced helices with a uniform radius. Many 
such helices were observed which differed by their radius R and their 
pitch h. The table gives some measured values of h and # microns and 
the calculated value of «, when tya—=h/27R. Inspection of these values 
indicates that they vary between wide limits. The value of the angle x 
gives a measure of how densely the helices are wound. 


h (microns) R (microns) x 


95 13° 30’ 
108 
72°8 
Tit 
100 
44 
ie 
166 


1 
2 
3 
4 
5 
6 
7 
8 


The directions of the axes of the helices were found to enclose small 
angles with the <110) directions. Next to helices with the expected 
characteristics, others, very irregular are also rather frequently observed. 

This is probably the result of differences in the concentration of 
vacancies or interstitials in the vicinity of the climbing dislocations. 
Photograph 1 (a) shows a helix where two turns have climbed faster 
than the others. Photograph 1(b) shows a helix with variable pitch 
and radius. We suppose that in this case the helix has been frozen in 
by the decoration in an energetically unfavourable position. Longer 
heat treatment would probably have produced a regular helix. 


3.2. Formation of Rows of Closed Loops and Special Features 


A hairpin shaped dislocation of which the longer sections have mixed 
(nearly screw) character can climb and produce a sequence of prismatic 
dislocation rings situated on the surface of a same cylinder. These 
rings are formed by the interaction of the two helices formed by climb 
of the dislocations AB and CD (fig. 1) (Amelinckx et al. 1957). 

Decorated patterns, where only closed loops are observed, can be 
explained by this mechanism. In most cases, the prismatic dislocations 
were followed by a helical dislocation (photograph 2 (a)). Such patterns 
can be formed by two different mechanisms. 

(a) A dislocation ABCDE consists of parts with respectively more 
Taylor (BC) and more screw character (AB-CE) (fig. 2 (a)). When 
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climb occurs, ABCD will act as a source of prismatic dislocations, DE 
will give rise to a helix, the resulting configuration is shown in fig. 2 (6). 
The axis of the closed loops will not coincide with the axis of the helix. 
The dislocation patterns on photograph 2 (b) could be the result of such a 
mechanism, although it was the only case where such a shift of the axis 
was observed. 


Formation of a helix, followed by a sequence of closed loops. ABCD will 
climb into two helices, which will interact and form a row of closed loops. 


DE will become a helix after climb. 

(b) We have a dislocation with a V-shaped form, where only one part 
had such a character that it could climb into a helix (fig. 3 (a)). During 
climb of the helix the part of the dislocation MN may cut the helical 
dislocation in the points A, B, so that a row of closed loops results (fig. 3 (6)). 
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MN can however have such a form that it cuts the helix in B instead of A. 
A dislocation as shown in fig. 3 (c) results. Those Lissajoux like configura- 
tions were many times observed, photographs 2 (c), (d) show some 
examples. 

It is evident that the same mechanism can occur with a helix AB 
which cuts during climb a neighbouring singular dislocation which has the 
same Burgers vector but of the opposite sign (fig. 4). 


Ze Z 


(b) (ey 


By climb only a part of a dislocation has adopted a helical shape. After 
intersection with the straight part, closed loops (6) and Lissajoux like 
configurations (c) are formed. 


Fig. 4 


oe 
> 
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(a) (b) 


Formation of closed loops by intersection of a helix with a singular dislocation 
line of the same Burgers vector and of suitable sign. 


3.3. Formation of Single Loops 

In many cases single closed loops were observed, they can be formed by 
different mechanisms (photographs 1 (a), 2 (c)). 

Tn our crystals, helices which possess only a single winding after climb 
were observed (photograph 2 (e)). These features were also seen in 
rocksalt (Amelinckx, unpublished results). Such a dislocation can 
interact with another dislocation with the same Burgers vector and of 
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suitable sign, and give rise to a closed loop (fig. 5). A second possibility 
is the precipitation of vacancies in a certain region. The so formed 
prismatic dislocation can climb further by enlarging its radius. 


Fig. 5 


(a) (b) 


Formation of a single closed dislocation line, by intersection of a helical disloca- 
tion with one loop with a singular dislocation with appropriate Burgers 
vector. 


3.4. Features Consisting of Different Helical Dislocations 


When only one, two or all three dislocations forming a node are nearly 
pure screws one or more helices with a common endpoint will result if 
climb is possible. As the Burgers vectors (110) in the (111) plane form 
angles of 120°, the screw dislocations will also form angles of 120° or 
60° and this angle will be preserved between the axis of the helices. 

In photograph 3 two sets of helical dislocations are shown. The 
measured angles are 111° and 78°, which is very acceptable, because 
the dislocations must have screw character without being pure screw. 
This is probably the case for the third dislocation for it is not decorated. 
When the three dislocations are nearly pure screws, three helices will 
result after climb. This has not yet been observed. 


3.5. Relation between the Etch Structure of (111) Cleavage Planes and 
Helical Dislocations 


Fluorite crystals annealed at high temperatures were cleaved and 
etched. Concentrated H,SO, was used as an etchant. Double rows of 
etchpits nearly parallel to (110) were many times observed. From the 
relative positions of the etchpits, we could deduce if it were the emergence 
points of a helical dislocation line or of a row of closed loops (photograph 4). 
The shape of the pits gives already some indication that the dislocation 
line has helical character. The pits are not perfectly hexagonal inverted 
pyramids ; the bottom points of the pits of the two opposite rows are 
displaced towards each other (fig. 6). This proves already that the 
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dislocation lines are oblique to the cleavage plane and are symmetrical 
to the centre line between the two rows. The definite proof could be 
given by continued etching. Photographs 5 (a), (b), (c) show three 
different stages. In (c) nearly all etchpits have disappeared ; the succes- 
sive positions of their centres come closer and closer and finally the 
centres disappear, as one continues the etching. As the instantaneous 
position of an etchpit is always centred on the emergence points of the 
dislocation, this allows the unambiguous conclusion that one has to do 
with half loops. 


SS = 


(6) 


Relation between the shape of the etchpits and the direction of the half loops 
responsible for their formation. 
(a) Cross section of the pits ; (b) top view of the pits. 


Photograph 6 (a) shows an etch pattern of two helical dislocation lines 
which make an angle of 115°. 

Patterns of double etchpits were also observed (photograph 6 (b)) it is 
clear that these are the emergence points of single closed loops as after 
continued etching the pits came closer and closer and disappeared 
altogether. 


3.6. Climb at Lower Temperatures 


At lower temperatures climb of what we think to be Taylor like disloca- 
tions was observed. We searched our crystals for decorated Bardeen— 
Herring sources. Photographs 7 (a), (b) show some patterns which were 
thought to be the first stages of a Bardeen—Herring mechanism. The 
most remarkable fact of lower temperature treatment was the occurrence 
of polygonized features. This indicates that climb at lower temperature 
proceeds in the anisotropic way. Photograph 8 (a), (b) show the most 
commonly observed feature. Photograph 8 (c) shows the first stage of 
a polygonized potential Bardeen—Herring source, and photograph 9 a 
polygonized helix seen on different levels. The bissectrices (fig. 7) of 
these polygonized patterns have a <110) direction. 

These features show that climb at lower temperatures seems to proceed 
along preferential directions. This polygonization can be explained in a 
similar way as the polygonization of growth fronts in crystal growth 
(Amelinckx 1953). The jogs in the climbing dislocation play the role of 
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the kinks on a growing step. We have to differentiate between intrinsic 
jogs which are due to the instantaneous orientation of the line with respect 
to the closed packed directions and those which are generated thermally. 
Densely packed directions in fluorite, which are (110), will contain only 
thermally generated jogs. At a certain temperature the total number of 
jogs will be smallest along these directions. The climb rate will be 
maximum perpendicular to the least densely packed direction. [112] is 
not a closely packed climb front and will climb fastest in a [110] direction. 


Fig. 7 
[10] 


=—-4-— —_ > [112] 


fi10] 


(a) (b) 


Patterns observed in photograph 8; the preferential climb directions in a (111) 
plane are indicated. 


At higher temperatures, the amount of thermally generated jogs will 
increase and the differences of climb velocities for different directions will 
be suppressed. 


3.7. Miscellaneous Features 


Interaction with other dislocations can also influence the shape of a 
helical dislocation. Photograph 10 shows some examples. In (a) two 
turns have come very close, this is normally impossible because of the 
high repulsion between two successive windings. In photograph 10 (6) 
we can see some features probably resulting from interaction. 
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CORRESPONDENCE 


Multipole Order of Lattice Cells of Ionic Crystals 


By F. G. Fumi and M. P. Tost 
Cattedra di Fisica Teorica dell’Université, Palermo, Italy 


[Received October 12, 1956] 


Lartice cells are neutral and thus act as multipoles. Here we determine 
the power of the distance with which the electrostatic potential decreases 
(effective multipole order) for some cubic lattice cells of the most common 
cubic ionic structures, but the method that we use is of general applicability. 
The problem of the multipole nature of lattice cells of ionic crystals had 
already been touched upon by Evjen (1932), but his treatment is not 
complete nor physically transparent, and his results are improperly quoted 
in the standard review article by Born and Goeppert—Mayer (1933). In 
particular, Born and Goeppert—Mayer imply that the effective multipole 
order of a lattice cell of a given structure depends only on the structure, 
while we shall see that it depends also on the cell. We shall also see that 
the results quoted by Born and Goeppert—Mayer are valid for the Evjen 
cells of the NaCl and CsCl structures (Evjen 1932), but not for the 
Evjen cell of the Zincblende structure (Evjen 1932). 

Let us consider a set of charges e, with position vectors r; with reference 
to a point C, which compose a neutral lattice cell of an ionic crystal. 
To analyze the potential that this cell generates in a distant point P 
we use the multipolar expansion 


P| (5) A Pe Cd 


where R is the position vector of point C with ne to point P and 
R,=R-r,. De 2 Deh y; and so on are the components of the multipole 


moments of the cell, which transform as polar tensors of rank 1, 2 and so on. 
To determine the multipole order of the cell, it is sufficient to consider 
its multipole moments with reference to its volume centre. The dipole 
term is not always zero in eqn. (1) (Kittel 1953), as Evjen supposed. 
It is certainly zero with the other terms of eqn. (1) which contain 
even powers of 1/R if one chooses a cell with an inversion centre, since 
then the cell cannot have any multipole moment represented by a polar 
tensor of odd rank. The dipole term vanishes also, with the quadrupole 
term, if the cell has cubic symmetry : indeed in this case the diagonal 
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components of the quadrupole moment are equal and the non-diagonal 
components are zero with the dipole components, while 1/R; satisfies 
Laplace equation. Thus, for cells of cubic symmetry without an inversion 
centre the first term in the multipolar expansion of the potential which 
may be different from zero is the term in 1/R*; for cells of cubic symmetry 
with an inversion centre the first non-vanishing term can be the term in 
1/Rk°. A closer analysis is necessary to decide whether these terms, for a 
given lattice cell, are actually different from zero. 

Let us consider, for instance, the Evjen lattice cell of the NaCl structure 
(Evjen 1932), which has cubic symmetry with an inversion centre (point 
group m3m). The direct-inspection method (Fumi 1952) allows one to 
see at once that for cubic symmetry the only non-vanishing components 
of the fourth-rank polar tensor which enters the term in 1/R® of eqn. (1) 
are w4=yt=z* and ay?=y?2?=222?. For the cell in question one has 
de27i=>de«2Zy2—0, since its faces and edges are neutral.. Thus its 


au v 
electrostatic potential does not go as 1/R°: it goes actually as 1/R” since 
dewey722A0. On the other hand, for the Hdjendahl lattice cell of the 


a 
NaCl structure (H6ojendahl 1938), which does not have an inversion centre 
(point group 437), the electrostatic potential goes as 1/R*, since Ye,7,4,2,40. 


v 
Similarly, the electrostatic potential of the Evjen lattice cells for the 
CsCl (Evjen 1932) and CaF, structures (point group m3m) goes as 1/R® 
since Sex'40, and that of the Evjen lattice cells for the Zincblende 


(Evjen 1 932) and cuprite structures (point group 43m) goes as 1/R* since 
Dery - 

The method that we have used to determine the effective multipole 
order of lattice cells of ionic crystals can obviously be used also to 
determine the effective multipole order of the Wigner—Seitz neutral 
polyhedra which enter Frank’s method of lattice summation (Frank 1950). 
One finds at once that the potential decreases with the fifth power of the 
distance both for the Wigner-Seitz cubes considered by Frank for the 
NaCl structure, and for the Wigner—Seitz truncated octahedra that one 
may consider to extend Frank’s method to the CsCl structure (Fumi and 
Tosi, to appear). 
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Quenching Vacancies in Aluminium 


By F. J. BrapsHaw and S. PEARSON 
Royal Aircraft Establishment, Farnborough 


[Received January 2, 1957] 


Ustne the same techniques as for previous work on platinum and gold 
(Bradshaw and Pearson 1956, 1957) we have quenched ‘ superpure ’ 
(99-995°%) aluminium wires of 0-2mm diameter from temperatures 
between 360°c and 530°c into water at ~2°c. The changes of resistance 
caused by quenching and annealing were measured in liquid nitrogen and 
are interpreted as being due to vacancies. 

In this work, owing to difficulties in welding aluminium to itself, the 
potential leads used to determine the resistance of the specimen were 
of 0:02mm diameter gold and not aluminium. But the annealed 
resistance remained virtually constant throughout the experiments and 
we are satisfied that the contamination of the specimen with gold was 
quite negligible. For temperatures up to 470°c the relation 


Apop=1240 exp (—0-76/k7’) » ohm cm 


could be fitted, where 4p, is the increase in resistivity resulting from 
quenching and 7’ is the absolute temperature (&7' expressed in ev). The 
error in the vacancy formation energy is thought to be about +0-04 ey. 
At temperatures above 470° the quenching was not rapid enough to 
retain all the vacancies. If in the absence of theoretical calculations 
for aluminium we assume that the effect of 1°% of vacancies will be 
comparable to that calculated in Cu, Ag or Au (~1-5u ohm cm), then 
the vacancy concentration in aluminium was 8 exp (—0-76/k7) giving 
~6x 10-4 at the melting point. 

The resistivity was annealed out at temperatures between —10° and 
+40°c (at 20°c the annealing time constant was ~4 min). The annealing 
curve shapes were similar to those found in gold (i.e. approximately 
those of second-order reactions) and there were indications of an increase 
in annealing rate with increasing 4p). The activation energy for annealing 
was found to be 0-444-0-03 ev and the mean value of »B when the 
resistivity increase had been reduced to 4p)/2 was 6x 10-1° where n is 
the probability per jump that a vacancy is lost and B is the entropy 
factor for vacancy migration. Since B is thought to be between 1 and 10, 
the average number of jumps made by a vacancy before it vanished was 
~101°. The grain size of the wire was roughly equal to the wire diameter 
from which it can quickly be shown that the grain boundaries and surfaces 
were not important as sinks for vacancies. 

From oscillograph camera records the initial quenching rates were 
found to be ~30°/m sec for all quenching temperatures and the total 
time to reach 2°c was between 15 and 22 m sec depending on quenching 
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temperature. As in the previous papers we may calculate from this and 
the annealing rate the transition temperature above which some vacancies 
should be lost because of insufficiently fast quenching. This agreed 
within the limits of error with the transition temperature estimated 
from the experimental variation of 4p, with temperature. 

Adding the formation and movement energies we get 1:2 ev for the 
activation energy for self-diffusion. Nowick (1951) has predicted 1-43 ev 
but Seymour (1953), using resonance methods, has measured 0-91 ey. 
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The Influence of Retardation on the London—Van der Waals Force 


By M. R. Aus, E. A. Power anp S. Zrenau 
University College, London 


[Received January 23, 1957] 


In a recent issue J. W. Leech has published a perturbation calculation of 
the retarded interaction between two atoms (Leech 1955). He concluded 
that the interaction energy at large distances varied as R~? in disagree- 
ment with a previous result of Casimir and Polder (1948). Large distance 
in this context means more precisely R>A, where A is the wavelength of 
the light emitted by one of the two identical atoms in its downward 
transition from the lowest excited state. On the other hand Casimir and 
Polder find that the usual London attractive energy proportional to R~6 
for R<A goes over to one proportional to k~’ at RSA which has recently 
found provisional experimental confirmation (Derjaguin ef al. 1956; 
Prosser and Kitchener 1956). 

In this problem the only physical dimensions, apart from # and ce, are 
the distance R, the static polarizability « of the neutral atom and A in 
the simplest model. An energy proportional to R © will therefore vary 
like A~1, one proportional to R~* will only contain A in so far as it enters 
into « and one varying like R~? must be proportional to A *. If this is 
considered in conjunction with the fact that the &* energy law is claimed 
for RSA, one sees that Leech’s result has an unphysical behaviour. 

We have therefore recalculated the retarded force by the conventional 
stationary perturbation method using the interaction Hamiltonian 


Pe é Cit. 
H=q,197(6,;—3h,f;)/Fe+ » {- ne pAt+ Since A) me) 


=1,2 


in dipole approximation, identical to that employed by Leech and by 
Casimir and Polder. We report disagreement with Leech, and obtain by 
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this symmetrical method the result of Casimir—Polder in closed form. 
The leading term in an expansion in powers of A/R is given by 
RTT 
ea) 23aP(g'y*(g?y? 
pr \ Oa =— : 2 
x Lies (k,w;)? — k+w reap ex [i(k+w), hi |= aR he ( ) 
The error in Leech’s work can be traced to his evaluation of the 


expectation value 


Op er 1c) 


which equals 
, fap yt] ee a ee 
— 8124q1)2(q2)2/(A2Kc) 2, [1+ (k,w;)?] kw(k+w) oe 


and although giving rise to an R~? dependence, this term cancels exactly 
with the leading terms of the remaining transverse photon contributions. 
That this method of evaluating the energy shift in lowest order 
~(e*) for two neutral atoms in interaction with the radiation field gives 
the same result as Casimir and Polder’s elegant but unsymmetrical 
procedure is to be expected on general grounds. 
Equation (2) can be obtained in a particularly straightforward manner 
with the aid of the reduced interaction Hamiltonian 
H=— YD qE{l) io Re eee eee) 


1=1,2 


where H;, is the transverse electric field operator. 
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REVIEWS OF BOOKS 


Irradiation of Colours and Luminescence. By K. Przipram. (Pergamon 
Press.) [Pp. 266+48 bibliography.] Price 63s. 


Tuts book gives an account of the coloration occurring in crystals as a result 
of exposure to different types of radiation, in particular those resulting from 
radioactivity. The associated luminescent properties are also dealt with. 
Part I discusses general aspects of coloration, especially in alkali halides, 
with emphasis on the structure sensitive properties, and includes some account 
of current transfer processes, the colour of glasses, and the effects of colloidal 
particles. Luminescence and radio-luminescence are also briefly treated. 

Part I deals with the colours and luminescence encountered in natural 
minerals, with special reference to rock salt and fluorite, and the part played by 
radiation in producing the effects. The impurities present in the crystals, 
and the sources of the radiation leading to colour and luminescence in natural 
minerals are also considered. 
__ The interpretation of the many detailed observations in terms of present 
ideas about colour centres and luminescence avoids the presentation of a mere 
catalogue of information, though here and there the escape from this a very 
narrow one. There is an extensive bibliography, dealing with the literature 
up to 1954. ‘ 

This book will be of general interest to the mineralogist and physicist, but 
is unlikely to compete at a specialist level with books already available dealing 
specifically with alkali halides, luminescent materials and radiation damage. 


Mathematics for Electronics with Applications. By Henry M. NopetMan 
and Freperick W. Smita. (McGraw-Hill Book Company, Inc.) [Pp. 391.] 
Price 52s. 6d. 


Ty planning this book, the authors have made an analysis of the mathematical 
content of a number of journals in the field of electronics to determine the 
incidence and importance of various mathematical processes. The subjects 
covered are based on this survey and include dimensional analysis, determinants, 
matrices, series and the solution of differential equations by both classical 
and Laplace transform methods. In addition, an outline of Boolean algebra 
is given and some practical examples of the use of calculus are included. 

Each subject is presented in two sections. The first is a concise and lucid 
exposition of the mathematical processes with illustration by worked numerical 
or abstract examples and the second is the application of the processes to 
specific electronic problems, again with a number of worked examples and 
related problems. Much of the practical material is taken from journals to 
which reference is made so that the reader may determine the context in which 
a particular problem has arisen. 

The clear way in which the authors have presented their subject and the 
constant recourse to practical problems for illustration will make this book 
of value to both the student who intends to specialise in electronics and the 
more experienced scientist who wishes to have at his disposal a number of 
labour-saving mathematical tools, especially fashioned for the solution of 


electronic problems. 
H. V. Beck. 
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Lectures on Rock Magnetism. By P. M. 8S. Buackerr. (Weizmann Science 
Press of Israel, Jerusalem.) [Pp.131.] Price $5 (Israeli). 


Rock magnetism is a subject touching on many fields of research in geology, 
geophysics, physics, mineralogy and chemistry. It has attracted increasing 
attention in the last ten years but is, as yet, little known and in a stage of rapid 
growth in which a textbook is out of place. Professor Blackett’s book does not 
claim to be a textbook but is a valuable account of the attitude and approach 
of one of the foremost minds of present day physics to a new and enthralling 
study. 

The possibilities of rock magnetism are discussed as well as the results to. 
date and the problems which have arisen in connection with polar wandering, 
continental drift and the reversal of the geomagnetic field. Naturally, in a 
somewhat personal book of this kind, the main emphasis is placed on the work 
of the author’s own laboratory but other contributions from many parts of 
the world are discussed and references quoted. Details are given of apparatus 
designed to meet some of the rather unusual requirements in these investiga- 
tions. 

The book will prove valuable both to workers within the field and, more 
particularly perhaps, in stimulating the interest of those in other fields. This 
stimulation should be of mutual advantage for, while the possible applications 
of rock magnetism are very wide, its study correspondingly requires techniques 
and ideas from many different bordering disciplines. ieee 


Les Moments Dipolaires. By J. Banriou. (Paris: Gouthier-Villars.) [Pp. 186.] 
Price : Unbound 3000 fr. ; Bound 3300 fr. 


THE contents of this excellent book are much wider than the title suggests. 
Besides treating the connection between dipolar moment and structure of 
molecules, it contains an up-to-date survey of the theory of dielectric constant 
and loss, as well as chapters on the dielectric properties of solids, and on 
experimental methods. In a subject in which so much incorrect theoretical 
work is done, the author has succeeded in presenting a correct and well- 
balanced picture. This book can be strongly recommended to physicists, 
chemists and engineers. 


Hoe: 


Digital Calculating Machines. By G. A. Montcomery. (Glasgow: Blackie 
and Son Limited.) [Pp. 262.] Price 30s. 


THE title of this book tells the strict truth, but it might lead some readers to 
expect a treatise devoted exclusively to modern electronic digital computers. 
In fact, most of the book deals’ with desk calculating machines ; punched 
card machines are dealt with in a single chapter, another chapter contains a 
description of the earlier electromechanical automatically sequenced machines, 
and the two final chapters deal respectively with the logical design of, and 
programming for, electronic sequenced machines. 

The author’s method in the earlier chapters is to divide machines into 
groups according to their type, and to give for each group a description of 
typical mechanisms, followed by a detailed discussion of the operating 
procedures appropriate to that group. Some readers may feel that the 
operating procedures are discussed in excessive detail, particularly as the 
author is concerned with operating pure and simple, and not with numerical 
analysis or with the overall planning of large-scale computations. Most 
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computers pick up their technique on the job, learning from others more 
experienced than themselves ; in some ways the present book is of the ‘ teach 
yourself’ type, adapted to the use of those who are denied such opportunities, 
The last chapter on programming for electronic digital computers, brief as it 
is, contains a sound introduction to the subject. 

The author is, undoubtedly, well-qualified in the field he covers, but the 
words “Member of the Association for Computing Machinery’, proudly 
displayed by the publisher on the dust-cover, can hardly be said to prove this, 
since that admirable Association is open to all who are interested in the subject, 
regardless of their qualifications. 


MV. W. 
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Transactions of the Symposium on Partial Differential Equations held at the 
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CORRIGENDUM 


Amendment to “A Method for the Mapping of Vector Potential 
Distributions in Axially Symmetrical Systems’, by the late Dr. G. 
LIEBMANN, published in Phil. Mag. (7), 61, 1143-1151, 1950. By F.C. 
Gatr (A.E.R.E., Aldermaston). 


EqQuatTIon (7) does not follow from (6), which is the wrong transformation 
to use. 
Instead, let 
r={(l—b)p}V-* (kx) 


Pot ee (k=1, ry>=constant) 


but leave z untransformed. 
Then, eqn. (9) becomes : 


For meshes of unit length in both r and z directions, the appropriate 
resistors can be specified by eqns. (10) (41), and the currents to be fed in 
by eqn. (11). 

In both figs. 1 and 2, R, and Rs are labelled the wrong way round. 
In fig. 3, the label on R, has been omitted. 


[The Editors do not hold themselves responsible for the views 
expressed by thew correspondents. | 
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Fatigue specimen repolished to remove approximately 4 microns. x 2500. 
Fig. 2 


A section through a group of damage bands. x 2000. 
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Fig. 4 


Grain boundaries and sub grains after electropolishing and etching for 20 sec 
in Wassermann’s reagent. x 3000. 


Section through approaching fatigue cracks As polished specimen surface after fatigue. x3 
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Fig. 1 


Pattern from 2 A of copper on rocksalt at 160°c [110] azimuth. 


12 A of copper on rocksalt at 160°C. 


Fig. 3 


160 & of copper on rocksalt at 160°c. 
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Pattern from 150 4 of gold on rocksalt at 295°c. 
[110] azimuth. Glancing angle of 0-7°. 


Or 
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As fig. 4, but glancing angle of 1-7°. 
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Light profile of steel slowly compressed 4°. Deformation band in steel rapidly 
< 1700. compressed 4°%. x 1700. 
Fig. 3 Hig. 4 


Coarse slip lines in slowly compressed. ferrite. 
7O00. 


apidly compressed (8%) steel, sectioned. 
electropolished and_ electrolytically 

etched. Oblique illumination. 
L700. 
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Structure of coarse slip step in slowly strained ferrite. — « 750( 


Slip steps following a grain boundary 
in rapidly deformed ferrite. > 7500. 


J. D. CAMPBELL et al. Phil. Mag. Ser. 8, Vol. 2, Pl. 


Fig. 9 


Precipitation of carbide along ferrite grain boundary. 7500. 


Fig. 10 


Carbide precipitation in ferrite grain. x 7500. 
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wins formed by impact in ferrite. x 1700. 
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Photograph 1 — - 


Photograph 3 


(a) A helix with two windings of a larger (a) and (b) Set of two helices forming 
radius. A single closed dislocation angles of +60° and 120°. 
line can be observed. 


(6) Helix with variable radius and pitch. 


Photograph 2 


(a) Helix followed by a row of closed loops. 
(b) Dislocation pattern which may be the result of the formation mechanism 
of fig. 2. - : Can 
(c) and (d) Helices followed by a sequence of closed loops. Lissajoux like 

configurations can be observed. 
(e) Helical dislocation with only one winding. 
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Photograph 4 
Photograph 5 


(a) Etch structure of a helical dislocation line. 
(b) Etch structure of a row of closed loops. 


Photograph 6 


The etch structure of a helical dis- 


(a) Etch pattern of two helical dislocations which location after successive etchings. 
form an angle of 115°. In (c) most of the pits have 
(6) Emergence points of isolated closed loops after disappeared. 


etching. 
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Photograph 7 Photograph 8 


Photograph 7. (a) and (4) First stages of the Bardeen—Herring climb mechanism. 


Photograph 8. (a) and (b) Climbing dislocations. 
(c) A polygonized Bardeen—Herring source. 


Photograph 10 
Photograph 9 


Polygonized helix seen on different depths. (a) Influence of the interaction of 
dislocations on the shape of 
helical dislocations. 

(b) Interaction between two rows 
of closed loops (a and y). 


